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Abstract 

Post-Newtonian relativistic theory of astronomical reference frames based on Ein- 
stein's general theory of relativity was adopted by General Assembly of the Interna- 
tional Astronomical Union in 2000. This theory is extended in the present paper by 
taking into account all relativistic effects caused by the presumable existence of a 
scalar field and parametrized by two parameters, P and 7, of the parametrized post- 
Newtonian (PPN) formalism. We use a general class of the scalar-tensor (Brans- 
Dicke type) theories of gravitation to work out PPN concepts of global and local 
reference frames for an astronomical N-body system. The global reference frame is 
a standard PPN coordinate system. A local reference frame is constructed in the 
vicinity of a weakly self-gravitating body (a sub-system of the bodies) that is a 
member of the astronomical N-body system. Such local inertial frame is required 
for unambiguous derivation of the equations of motion of the body in the field of 
other members of the N-body system and for construction of adequate algorithms 
for data analysis of various gravitational experiments conducted in ground-based 
laboratories and/or on board of spacecrafts in the solar system. 

We assume that the bodies comprising the N-body system have weak gravitational 
field and move slowly. At the same time we do not impose any specific limitations 
on the distribution of density, velocity and the equation of state of the body's mat- 
ter. Scalar-tensor equations of the gravitational field are solved by making use of 
the post-Newtonian approximations so that the metric tensor and the scalar field 
are obtained as functions of the global and local coordinates. A correspondence 
between the local and global coordinate frames is found by making use of asymp- 
totic expansion matching technique. This technique allows us to find a class of the 



post-Newtonian coordinate transformations between the frames as well as equations 
of translational motion of the origin of the local frame along with the law of rel- 
ativistic precession of its spatial axes. These transformations depend on the PPN 
parameters (5 and 7, generalize general relativistic transformations of the lAU 2000 
resolutions, and should be used in the data processing of the solar system gravi- 
tational experiments aimed to detect the presence of the scalar field. These PPN 
transformations are also applicable in the precise time-keeping metrology, celestial 
mechanics, astrometry, geodesy and navigation. 

We consider a multipolar post-Newtonian expansion of the gravitational and 
scalar fields and construct a set of internal and external gravitational multipoles 
depending on the parameters (3 and 7. These PPN multipoles generalize the Thorne- 
Blanchet-Damour multipoles defined in harmonic coordinates of general theory of 
relativity. The PPN multipoles of the scalar-tensor theory of gravity are split in three 
classes - active^ conformal, and sca/ar multipoles. Only two of them are algebraically 
independent and we chose to work with the conformal and active multipoles. We 
derive the laws of conservations of the multipole moments and show that they must 
be formulated in terms of the conformal multipoles. We focus then on the law of 
conservation of body's linear momentum which is defined as a time derivative of 
the conformal dipole moment of the body in the local coordinates. We prove that 
the local force violating the law of conservation of the body's linear momentum 
depends exclusively on the active multipole moments of the body along with a few 
other terms which depend on the internal structure of the body and are responsible 
for the violation of the strong principle of equivalence (the Nordtvedt effect). 

The PPN translational equations of motion of extended bodies in the global 
coordinate frame and with all gravitational multipoles taken into account are derived 
from the law of conservation of the body's linear momentum supplemented by the 
law of motion of the origin of the local frame derived from the matching procedure. 
We use these equations to analyze translational motion of spherically-symmetric and 
rigidly rotating bodies having finite size. Spherical symmetry is defined in the local 
frame of each body through a set of conditions imposed on the shape of the body 
and the distribution of its internal density, pressure and velocity field. We prove 
that our formalism brings about the parametrized post-Newtonian EIH equations 
of motion of the bodies if the finite-size effects are neglected. Analysis of the finite- 
size effects reveal that they are proportional to the parameter P coupled with the 
second and higher-order rotational moments of inertia of the bodies. The finite-size 
effects in the translational equations of motion can be appreciably large at the latest 
stage of coalescence of binary neutron stars and can be important in calculations of 
gravitational waveform templates for the gravitational-wave interferometers. 

The PPN rotational equations of motion for each extended body possessing an 
arbitrary multipolar structure of its gravitational field, have been derived in body's 
local coordinates. Spin of the body is defined phenomenologically in accordance 
with the post-Newtonian law of conservation of angular momentum of an isolated 
system. Torque consists of a general relativistic part and the PPN contribution 
due to the presence of the scalar field. The PPN scalar-field-dependent part is 
proportional to the difference between active and conformal dipole moments of the 
body which disappears in general relativity. Finite-size effects in rotational equations 



of motion can be a matter of interest for calculating gravitational wave radiation 
from coalescing binaries. 
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1 Notations 



1.1 General Conventions 

Greek indices a, /3, 7, ... run from to 3 and mark space-time components of four- dimensional 
objects. Roman indices 2, j, /c, ... run from 1 to 3 and denote components of three-dimensional 
objects (zero component belongs to time). Repeated indices mean the Einstein summation 
rule, for instance, = A^Bq + A^Bi + A^B2 + A^B^ and T| = T\ + T| + T|, etc. 

Minkowski metric is denoted rjap = diag(— 1, +1, +1, +1). Kronecker symbol (the unit ma- 
trix) is denoted 6ij = diag(l, 1, 1). Levi-Civita fully antisymmetric symbol is Sijk such that 
^123 = +1- Kronecker symbol is used to rise and lower Roman indices. Complete metric 
tensor g^p is used to rise and lower the Greek indices in exact tensor equations whereas 
the Minkowski metric riap is employed for rising and lowering indices in the post-Newtonian 
approximations 

Parentheses surrounding a group of Roman indices mean symmetrization, for example, 
A(^ij) = {l/2){Aij + Aji). Brackets around two Roman indices denote antisymmetrization, 
that is y4[jj] = (1/2) (Ajj — Aji). Angle brackets surrounding a group of Roman indices de- 
note the symmetric trace-free (STF) part of the corresponding three-dimensional object, for 
instance, 

^<ij> '^(ij) g^ij'^fcfe ) -^<ijk> ^{ijk) '^^ij-^kpp '^^jk-^ipp "^ikAjpp . 

We also use multi-index notations, for example. 
Sum over multi-indices is understood as 

A,n^ — A. . . r^hh-'-ii p ^ rpbL-1 _ p rpbiii2...ii_i 

Comma denotes a partial derivative, for example, (f)^a = d(j)/dx°', where (p^ = c~^d(j)/dt, (f),i = 
d(j)/dx\ and semicolon T"^ denotes a covariant derivative. L-order partial derivative with 
respect to spatial coordinates is denoted = di^...d^. Other conventions are introduced as 
they appear in the text. We summarize these particular conventions and notations in the 
next section for the convenience of the readers. 



1.2 Particular Conventions and Symbols Used in the Paper 



Symbol 


Description 


Equation(s) 




physical (Jordan-Fierz frame) metric 
tensor 


(3.1.4) 


9fiu 


conformal (Einstein frame) metric ten- 
sor 


(3.1.6) 



Symbol 


Description 


Equation(s) 


g 


the determinant of g^j_y 


(3.1.1) 


g 


the determinant of g^y 


(3.4.2) 




the Minkowski (fiat) metric tensor 


(3.3.4) 


-pa 


the Christoffel symbol 


(3.1.3) 




the Ricci tensor 


(3.1.2) 


R 


the Ricci scalar 


(3.1.1) 




the conformal Ricci tensor 


(3.1.7) 


T 


the energy-momentum tensor of matter 


(3.1.2) 


T 


the trace of the energy-momentum ten- 
sor 


(3.1.2) 




the scalar field 


(3.1.1) 


00 


the background value of the scalar field 




(3.3.1) 


c 


the dimensionless perturbation of the 
scalar field 


(3.3.1) 


^(0) 


the coupling function of the scalar field 


(3.1.1) 




the Laplace-Beltrami operator 


(3.1.3) 


□ 


the D'Alembert operator in the 
Minkowski space-time 


(3.5.5) 


p 


the density of matter in the co-moving 
frame 


(3.2.1) 


p* 


the invariant (Fock) density of matter 


(3.3.18) 


n 


the internal energy of matter in the co- 
moving frame 


(3.2.1) 




the tensor of (anisotropic) stresses of 
matter 


(3.2.1) 




the 4-velocity of matter 


(3.2.1) 




the 3-dimensional velocity of matter in 
the global frame 


(3.3.15) 




the asymptotic value of the coupling 
function 6^(0) 


(3.3.2) 




the asymptotic value of the derivative 
of the coupling function ^(0) 


(3.3.2) 


C 


the ultimate speed of general and spe- 
cial theories of relativity 


(3.1.1) 



Symbol 


Description 


Equation(s) 


e 


a small dimensional parameter, e = 1/c 


(3.3.4) 




the metric tensor perturbation, g^^, — 


(4.1.1) 


(n) 

h fj,u 


the metric tensor perturbation of order 
e m the post-iNewtonian expansion or 
the metric tensor 


(3.3.4) 


N 


(2) 

a shorthand notation for hoo 


(3.3.9) 


L 


(4) 

a shorthand notation for hoo 


(3.3.9) 


Ni 


(1) 

a shorthand notation for hoi 


(3.3.9) 


Li 


(3) 

a shorthand notation for hoi 


(3.3.9) 


Hij 


(2) 

a shorthand notation for hij 


(3.3.9) 


H 


(2) 

a shorthand notation for /iytfc 


(3.3.9) 


N, L 


shorthand notations for perturbations 
of conformal metric cjf^i, 


(5.4.1) 


7 


the 'space-curvature' PPN parameter 


(3.5.1) 




the 'non-linearity' PPN parameter 


(3.5.2) 


7] 


the Nordtvedt parameter, t] = 4/5—7—3 


(5.3.1) 


G 


the observed value of the universal 
gravitational constant 


(3.5.4) 


g 


the bare value of the universal gravita- 
tional constant 


(3.3.3), (3.5.4) 




the global coordinates with x° = ct and 






the local coordinates with = cu and 
= w 




u 


the Newtonian gravitational potential 
in the global frame 


(4.2.1) 


t/(A) 


the Newtonian gravitational potential 
of body A m the global frame 


(4.2.7) 




a vector potential in the global frame 


(4.2.4) 




a vector potential of body A in the 
global frame 


(4.2.7) 




various special gravitational potentials 
in the global frame 


(4.2.2), (4.2.6) 



Symbol 


Description 


Equation(s) 




potentials of the physical metric in the 
global frame 


(5.2.1), 5.2.2) 


(7, Cr* 


the active mass and current-mass den- 
sities in the global frame 


(5.2.3), (5.2.4) 


I<L> 


the active Thorne-Blanchet-Damour 
mass multipole moments in the global 
frame 


(5.2.12) 


S<L> 


the active spin multipole moments in 
the global frame 


(5.2.13) 


V 


potential of the scalar field in the global 
frame 


(5.3.1) 


a 


scalar mass density in the global frame 


(5.3.2) 


I<L> 


scalar mass multipole moments in the 
global frame 


(5.3.6) 


V 


gravitational potential of the conformal 
metric in the global frame 


(5.4.1) 


a 


the conformal mass density in the 
global frame 


(5.4.2) 


I<L> 


the conformal mass multipole moments 
in the global frame 


(5.4.6) 


M 


conserved mass of an isolated system 


(5.5.4) 


pi 


conserved linear momentum of an iso- 
lated system 


(5.5.5) 


§^ 


conserved angular momentum of an iso- 
lated system 


(5.5.6) 




integral of the center of mass of an iso- 
lated system 


(5.5.7) 


A 


symbols with the hat stand for quanti- 
ties in the local frame 




(B) 


sub-index referring to the body and 
standing for the internal solution in the 
local frame 


(6.2.1), (6.2.2) 


(E) 


sub-index referring to the external with 
respect to (B) bodies and standing for 
the external solution in the local frame 


(6.2.1), (6.2.2) 


(C) 


sub-index standing for the coupling 
part of the solution in the local frame 


(6.2.2) 



Symbol 


Description 


Equation(s) 




external STF multipole moments of the 
scalar field 


(6.2.4) 


Ql 


external STF gravitoelectric multipole 
moments of the metric tensor 


(6.2.19) 


c. 


external STF gravitomagnetic multi- 
pole moments of the metric tensor 


(6.2.21) 




other sets of STF multipole moments 
entering the general solution for the 
space-time part of the external local 
metric 


(6.2.21) 


Yl, -Bi, -Dl, -El, -Fl, 


STF multipole moments entering the 
general solution for the space-space 
part of the external local metric 


(6.2.22) 




linear and angular velocities of kine- 
matic motion of the local frame; we put 
them to zero throughout the rest of the 
paper 


(6.2.19), (6.2.20) 




3-dimensional velocity of matter in the 
local frame 


(6.2.10) 




active Thorne-Blanchet-Damour STF 
mass multipole moments of the body 
in the local frame 


(6.3.1) 


Ob 


active mass density of body B in the 
local frame 


(6.3.2) 


X, 


scalar STF mass multipole moments of 
the body in the local frame 


(6.3.3) 




scalar mass density of body B in the 
local frame 


(6.3.4) 


Xl 


conformal STF mass multipole mo- 
ments of the body in the local frame 


(6.3.5) 




conformal mass density of body B in 
the local frame 


(6.3.6) 




current mass density of body B in the 
local frame 


(6.3.7) 




spin STF multipole moments of the 
body in the local frame 


(6.3.8) 




Relativistic corrections in the post- 
Newtonian transformation of time and 
space coordinates 


(7.2.1), (7.2.2) 



Symbol 


Description 


Equation(s) 


'T'^ 71^ 
■^Bl ^B1 


position, velocity and acceleration of 
the body's center of mass with respect 
to the global frame 


(7.2.2), (7.2.6), (7.2.9) 


Fh„ 

B 


X* — x'ft), i.e. the spacial coordinates 
taken with respect to the center of mass 
of body B in the global frame 


(7.2.2) 


^<L> 


functions appearing in the relativistic 
transformation of time 


(7.2.6), (7.2.13) 


T^<L>i ^ <L>1 £<L> 


functions appearing in the relativistic 
transformation of spacial coordinates 


(7.2.14) 




matrix of transformation between local 
and global coordinate bases 


(7.3.1) 


03, 2), m\ ^\ 


PN corrections in the matrix of trans- 
formation A^^ 


(7.3.3)-(7.3.6) 


U, U\ etc. 


external gravitational potentials 


(8.2.1)-(8.2.4) 


Ul{Xb). U\{Xb) 


/-th spatial derivative of an external po- 
tential taken at the center of mass of 
body B 


(8.4.1) 

V / 




PN correction in the formula of match- 
ing of the local Newtonian potential 


(8.3.10) 




the matrix of relativistic precession of 
local coordinates with respect to global 
coordinates 


(8.5.33) 




Newtonian-type mass, center of mass, 
and linear momentum of the body in 
the local frame 


(9.1.1)-(9.1.3) 


M 


general relativistic PN mass of the 
body in the local frame 


(9.3.3) 




active mass of the body in the local 
frame 


(9.3.2) 




conformal mass of the body in the local 
frame 


(9.3.1) 


J{2) 


rotational moment of inertia of the 
body in the local frame 


(9.3.4) 




a set or S i r multipole moments 


(9.3.8) 


•pi 


PN linear momentum of the body in 
the local frame 


(9.3.9) 




scalar-tensor PN correction to 


(9.4.2) 



Symbol 


Description 


Equation(s) 




conformal anisotropic mass of the body 
in the local frame 


(9.4.4) 


Ft, A¥%, ¥i^, A¥i^ 


gravitational forces in the expression 
for Qi 


(9.4.5)-(9.4.7) 

V / V / 




the bare post-Newtonian definition of 
the angular momentum (spin) of a 
body 


(10.1.3) 




the post-Newtonian torque in equa- 
tions of rotational motion 


(10.2.2) 


AT' 


the post-Newtonian correction to the 
torque 


(10.2.3) 


AS' 


the post-Newtonian correction to the 
bare spin iS* 


(10.2.4) 




velocity-dependent multipole moments 


(10.2.5) 


sx 


the (measured) post-Newtonian spin of 
the body 


(10.2.6) 


r 


radial space coordinate in the body's 
local frame, r = it> 


(11.1.1) 




angular velocity of rigid rotation of the 
body B referred to its local frame 


(11.1.3) 


j{2l) 


l-th rotational moment of inertia of the 
body B 


(11.1.10) 




multipole moments of the multipolar 
expansion of the Newtonian potential 
in the global coordinates 


(11.2.1) 


Rb 


-Rb , where Rg = x — Xg 


(11.3.1) 


r)i 

^BC 


•^C -^B 


(11.4.3) 


pi pi pi pi Xpi 

AT' EIH^ ^ S'> IGfl' " XGR 


forces from the equation of motion of 
spherically-symmetric bodies 


(11.4.11)-(11.4.15) 


^B 


Nordtvedt's gravitational mass of the 
body B 


(11.4.16) 



2 Introduction 



2. 1 General Outline of the Paper 

This paper consists of 11 sections and 3 appendices. In this section we give a brief introduction 
to the problem of relativistic reference frames and describe our motivations for doing this 
work. Section 3 outhnes the statement of the problem, field equations and the principles of 
the post-Newtonian approximations. Section 4 is devoted to the construction of the global 
(barycentric) reference frame which is based on the solution of the field equations in entire 
space. We make a multipolar expansion of the gravitational field in the global coordinates 
and discuss the post-Newtonian conservation laws in section 5. Section 6 is devoted to the 
construction of local coordinates in the vicinity of each body being a member of N-body 
system. A general structure of the post-Newtonian coordinate transformations between the 
global and local coordinate frames is discussed in section 7. This structure is specified in 
section 8 where the matching procedure between the global and local coordinates is employed 
on a systematic ground. We use results of the matching procedure in section 9 in order 
to derive PPN translational equations of motion of the extended bodies in the N-body 
system. PPN equations of rotational motion of each body are derived in section 10. These 
general equations are applied to the case of motion of spherically-symmetric bodies which 
is considered in section 11. Appendix A gives solution of the Laplace equation in terms 
of scalar, vector and tensor harmonics. Appendix B provides with explicit expressions for 
calculation of the Christoffel symbols and Riemann tensor in terms of the post-Newtonian 
perturbation of the metric tensor. Appendix C compares our results with those obtained by 
Klioner and Soffel [1] by making use of a different approach. 

2.2 Motivations and Historical Background 

General theory of relativity is the most powerful theoretical tool for experimental gravita- 
tional physics both in the solar system and outside of its boundaries. It passed all tests with 
unparallel degree of accuracy [2,3,4]. However, alternative theoretical models are required 
for deeper understanding of the nature of space-time gravitational physics and for studying 
possible violations of general relativistic relationships which may be observed in near-future 
gravitational experiments designed for testing the principle of equivalence [5], mapping as- 
trometric positions of stars in our galaxy with micro-arcsecond precision [6] and searching 
for extra-solar planets [7], testing near-zone relativistic effects associated with finite speed 
of propagation of gravitational fields [8,9,10], detection of freely propagating gravitational 
waves [11] -[12], etc. 

Recently, International Astronomical Union (lAU) has adopted new resolutions [13] -[15] 
which lay down a self-consistent general relativistic foundation for further applications in 
modern geodesy, fundamental astrometry, and celestial mechanics in the solar system. These 
resolutions combined two independent approaches to the theory of relativistic reference 
frames in the solar system developed in a series of publications of various authors ^ . The 

^ These approaches are called Brumberg-Kopeikin (BK) and Damour-Soffel-Xu (DSX) formalisms. 
The reader is invited to review [15] for full list of bibliographic references. 



goal of the present paper is to incorporate the parametrized post-Newtonian (PPN) formal- 
ism [16] -[20] to the lAU theory of general relativistic reference frames in the solar system. 
This will extend domain of applicability of the resolutions to more general class of gravity 
theories. Furthermore, it will make the lAU resolutions fully compatible with the JPL equa- 
tions of motion used for calculation of ephemerides of major planets, Sun and Moon. These 
equations depend on two PPN parameters (3 and 7 [21] and they are presently compatible 
with the lAU resolutions only in the case of /5 = 7 = 1. 

PPN parameters (3 and 7 are characteristics of a scalar field which perturbs the metric tensor 
and makes it different from general relativity. Scalar fields has not yet been detected but they 
already play significant role in modern physics. This is because scalar fields help us to explain 
the origin of masses of elementary particles [22], to solve various cosmological problems [23] 
-[25], to disclose the nature of dark energy in the universe [26], to develop a gauge-invariant 
theory of cosmological perturbations [27,28] joining in a very natural way the ideas contained 
in the original gauge-invariant formulation proposed by Bardeen [29] ^ with a coordinate- 
based approach of Lifshitz [31,32]. In the present paper we employ a general class of the 
scalar-tensor theories of gravity initiated in the pioneering works by Jordan [33,34], Fierz 
[35] and, especially. Brans and Dicke [36] -[38] ^ . This class of theories is based on the 
metric tensor gap representing gravitational field and a scalar field that couples with the 
metric tensor through the coupling function 9{(j)) which we keep arbitrary. We assume that 
and 9{(j)) are analytic functions which can be expanded about their cosmological background 
values and 9. Existence of the scalar field brings about dependence of the universal 
gravitational constant G on the background value of the field which can be considered as 
constant on the time scale much shorter than the Hubble cosmological time. 

Our purpose is to develop a theory of relativistic reference frames in an N-body problem 
(solar system) with two parameters (3 and 7 of the PPN formalism. There is a principal 
difficulty in developing such a theory associated with the problem of construction of a local 
reference frame in the vicinity of each self-gravitating body (Sun, Earth, planet) comprising 
the N-body system. Standard textbook on the PPN formalism [20] does not contain solution 
of this problem in the post-Newtonian approximation because the original PPN formalism 
was constructed in a single, asymptotically flat, global coordinate chart (PPN coordinates) 
covering the entire space-time and having the origin at the barycenter of the solar system. 
PPN formalism admits existence of several fields which are responsible for gravity - scalar, 
vector, tensor, etc. After imposing boundary conditions on all these fields at infinity the 
standard PPN metric tensor combines their contributions all together in a single expression 
so that they get absorbed to the Newtonian and other general relativistic potentials and their 
contributions are strongly mixed up. It becomes technically impossible to disentangle the 
fields in order to find out relativistic space-time transformation between local frame of a self- 
gravitating body (Earth) and the global PPN coordinates which would be consistent with the 
law of transformation of the fields imposed by each specific theory of gravitation. Rapidly 
growing precision of optical and radio astronomical observations as well as calculation of 
relativistic equations of motion in gravitational wave astronomy urgently demands to work 
out a PPN theory of such relativistic transformations between the local and global frames. 



^ See [30] for review of more recent results related to the development of Bardeen's theory of 
cosmological perturbations. 

^ For a well-written introduction to this theory and other relevant references can be found in [20] 
and [39]. 



It is quite straightforward to construct the post-Newtonian Fermi coordinates along a world 
line of a massless particle [40]. Such approach can be directly applied in the PPN formal- 
ism to construct the Fermi reference frame around a world line of, for example, an artificial 
satellite. However, account for gravitational self-field of the particle (extended body) changes 
physics of the problem and introduces new mathematical aspects to the existing procedure of 
construction of the Fermi frames as well as to the PPN formalism. To the best of our knowl- 
edge only two papers [1,41] have been published so far by other researchers where possible 
approaches aimed to derive the relativistic transformations between the local (geocentric, 
planetocentric) and the PPN global coordinate frame were discussed in the framework of the 
PPN formalism. The approach proposed in [41] is based on the formalism that was originally 
worked out by Ashby and Bertotti [42,43] in order to construct a local inertial frame in the 
vicinity of a self-gravitating body that is a member of an N-body system . In the Ashby- 
Bertotti formalism the PPN metric tensor is taken in its standard form [20] and it treats 
all massive bodies as point-like monopole massive particles without rotation. Construction 
of a local inertial frame in the vicinity of such massive particle requires to impose some 
specific restrictions on the world line of the particle. Namely, the particle is assumed to be 
moving along a geodesic defined on the "effective" space-time manifold which is obtained 
by elimination of the body under consideration from the expression for the standard PPN 
metric tensor. This way of introduction of the "effective" manifold is not defined uniquely 
bringing about an ambiguity in the construction of the "effective" manifold [45]. Moreover, 
the assumption that bodies are point-like and non-rotating is insufficient for modern geodesy 
and relativistic celestial mechanics. For example, planets in the solar system and stars in 
binary systems have appreciable rotational speeds and noticeable higher-order multipole mo- 
ments. Gravitational interaction of the multipole moments of a celestial body with external 
tidal field does not allow the body to move along the geodesic line [45]. Deviation of the 
body's center-of-mass world line from the geodesic motion can be significant and important 
in numerical calculations of planetary ephemerides (see, e.g., [46] and discussion on page 
307 in [47]) and must be taken into account when one constructs a theory of the relativistic 
reference frames in the N-body system. 

Different approach to the problem of construction of a local (geocentric) reference frame in 
the PPN formalism was proposed in the paper by Klioner and Soffel [1]. These authors have 
used a phenomenological approach which does not assume that the PPN metric tensor in 
local coordinates must be a solution of the field equations of a specific theory of gravity. 
The intention was to make the Klioner-Soffel formalism as general as possible. To this end 
these authors assumed that the structure of the metric tensor written down in the local 
(geocentric) reference frame must have the following properties: 

A. gravitational field of external bodies (Sun, Moon, planets) is represented in the vicinity of 
the Earth in the form of tidal potentials which should reduce in the Newtonian limit to 
the Newtonian tidal potential, 

B. switching off the tidal potentials must reduce the metric tensor of the local coordinate 
system to its standard PPN form. 

Direct calculations revealed that under assumptions made in [1] the properties (A) and (B) 
can not be satisfied simultaneously. This is a direct consequence of the matching procedure 
applied in [1] in order to transform the local geocentric coordinates to the global barycentric 

^ Fukushima (see [44] and references therein) developed similar ideas independently by making use 
of a slightly different mathematical technique. 



ones. More specifically, at each step of the matching procedure four kinds of different terms 
in the metric tensors have been singling out and equating separately in the corresponding 
matching equations for the metric tensor (for more details see page 024019-10 in [1]): 

• the terms depending on internal potentials of the body under consideration (Earth); 

• the terms which are functions of time only; 

• the terms which are linear functions of the local spatial coordinates; 

• the terms which are quadratic and higher-order polynomials of the local coordinates. 

These matching conditions are implemented in order to solve the matching equations. It is 
implicitly assumed in [1] that their application will not give rise to contradiction with other 
principles of the parametrized gravitational theory in the curved space-time. 

We draw attention of the reader to the problem of choosing the right number of the matching 
equations. In general theory of relativity the only gravitational field variable is the metric 
tensor. Therefore, it is necessary and sufficient to write down the matching equations for 
the metric tensor only. However, alternative theories of gravity have additional fields (scalar, 
vector, tensor) which contribute to the gravitational field as well. Hence, in these theories one 
has to work out matching equations not only for the metric tensor but also for the additional 
fields. This problem has not been discussed in [1] which assumed that it will be sufficient 
to solve the matching equations merely for the metric tensor in order to obtain complete 
information about the structure of the parametrized post-Newtonian transformation from the 
local to global frames. This might probably work for some (yet unknown) alternative theory 
of gravity but the result of matching would be rather formal whereas the physical content of 
such matching and the degree of applicability of such post-Newtonian transformations will 
have remained unclear. In the present paper we rely upon quite general class of the scalar- 
tensor theories of gravity and consistently use the matching equation for the metric tensor 
along with that for the scalar field which are direct consequences of the field equations. We 
have found that our results diverge pretty strongly from the results of Klioner-Soffel's paper 
[1]. This divergence is an indication that the phenomenological (no-gravity- field equations) 
Klioner - Soffel approach to the PPN formalism with local frames taken into account has 
too many degrees of freedom so that the method of construction of the parametrized metric 
tensor in the local coordinates along with the PPN coordinate transformations proposed in 
[1] can not fix them uniquely. Phenomenological restriction of this freedom can be done in 
many different ways ad liberum, thus leading to additional (researcher-dependent) ambiguity 
in the interpretation of relativistic effects in the local (geocentric) reference frame. 

We have already commented that in Klioner-Soffel approach [1] the metric tensor in the local 
coordinates is not determined from the field equations ^ but is supposed to be found from 
the four matching conditions indicated above. However, the first of the matching conditions 
requires that all internal potentials generated by the body's (Earth's) matter can be fully 
segregated from the other terms in the metric tensor. This can be done, for example, in 
general relativity and in the scalar-tensor theory of gravity as we shall show later in the 
present paper. However, complete separation of the internal potentials describing gravita- 
tional field of a body under consideration from the other terms in matching equations may 
not work out in arbitrary alternative theory of gravity. Thus, the class of gravity theories to 
which the first of the matching conditions can be applied remains unclear and yet has to be 
identified. Further discussion of the results obtained by Klioner and Soffel is rather technical 



^ Observe the presence of a free function ^' in Eq. (3.33) of the paper [1]. 



and deferred to appendix C. 



Our point of view is that in order to eliminate any inconsistency and undesirable ambigui- 
ties in the construction of the PPN metric tensor in the local reference frame of the body 
under consideration and to apply mathematically rigorous procedure for derivation of the 
relativistic coordinate transformations from the local to global coordinates, a specific theory 
of gravity must be used. The field equations in such known and the number of 

functions entering the PPN metric tensor in the local coordinates is exactly equal to the 
number of matching equations. Hence, all of them can be determined unambiguously. Thus, 
we propose to build a parametrized theory of relativistic reference frames in an N-body 
system by making use of the following procedure: 

(1) Chose a class of gravitational theories with a well-defined system of field equations. 

(2) Impose a specific gauge condition on the metric tensor and other fields to single out 
a class of global and local coordinate systems and to reduce the field equations to a 
solvable form. 

(3) Solve the reduced field equations in the global coordinate system x°' = {x^,x'^) by 
imposing fall-off boundary conditions at infinity. 

(4) Solve the reduced field equations in the local coordinate system = defined 
in the vicinity of a world line of the center-of-mass of a body. This will give N local 
coordinate systems. 

(5) Make use of the residual gauge freedom to eliminate nonphysical degrees of freedom 
and to find out the most general structure of the space-time coordinate transformation 
between the global and local coordinates. 

(6) Transform the metric tensor and the other fields from the local coordinates to the global 
ones by making use of the general form of the coordinate transformations found at the 
previous step. 

(7) Derive from this transformation a set of matching (first-order differential and/or alge- 
braic) equations for all functions entering the metric tensor and the coordinate trans- 
formations. 

(8) Solve the matching equations and determine all functions entering the matching equa- 
tions explicitly. 

This procedure works perfectly in the case of general relativity [15] and is valid also in the 
class of the scalar-tensor theories of gravity as we shall show in the present paper. We do not 
elaborate on this procedure in the case of vector-tensor and tensor-tensor theories of gravity. 
This problem is supposed to be solved somewhere else. 

The scalar-tensor theory of gravity employed in this paper operates with one tensor, Qa/B, 
and one scalar, 0, fields. The tensor field Qap is the metric tensor of the Riemannian space- 
time manifold. The scalar field is not fully independent and is generated by matter of 
the gravitating bodies comprising an N-body system. We assume that the N-body system 
(solar system, binary star) consists of extended bodies which gravitational fields are weak 
everywhere and characteristic velocity of motion is slow. These assumptions allow us to use 
the post-Newtonian approximation (PNA) scheme developed earlier by various researchers 
® in order to find solutions of the scalar-tensor field equations with non-singular distribution 

^ PNA solves the gravity field equations by making use of expansions with respect to the weak-field 
and slow-motion parameters. The reader is referred to the cornerstone works [48] -[55] which reflect 
different aspects of the post-Newtonian approximations. 



of matter in space. The method, we work out in the present paper, is a significant extension 
and further improvement of the general relativistic calculations performed in our previous 
papers [45,47,56,57,58,59,60,61]. It takes into account the post-Newtonian definition of mul- 
tipole moments of an isolated self-gravitating body (or a system of bodies) introduced by Kip 
Thorne [62] which has been mathematically elucidated and further developed by Blanchet 
and Damour [63] (see also [64] and references therein). We do not specify the internal struc- 
ture of the bodies so that our consideration is not restricted with the case of a perfect fluid 
as it is usually done in the PPN formalism. 



3 Statement of the Problem 

3. 1 Field Equations in the Scalar- Tensor Theory of Gravity 

The purpose of this paper is to develop a relativistic theory of reference frames for N-body 
problem in the PPN formalism which contains 10 parameters [20]. Michelson-Morley and 
Hughes-Drever type experiments strongly restricted possible violations of the local isotropy 
of space whereas Eotvos-Dicke-Braginsky type experiments verified a weak equivalence prin- 
ciple with very high precision [20]. These remarkable experimental achievements and modern 
theoretical attempts to unify gravity with other fundamental fields strongly restrict class of 
viable alternative theories of gravity and very likely reduce the number of parameters of the 
standard PPN formalism [20] to two - /3 and 7 . These parameters appear naturally in 
the class of alternative theories of gravity with one or several scalar fields [20,65] which can 
be taken as a basis for making generalization of the lAU resolutions on relativistic refer- 
ence frames. For this reason, we shall work in this paper only with the class of scalar-tensor 
theories of gravity assuming that additional vector and/or tensor fields do not exist. For sim- 
plicity we focus on the case with one real-valued scalar field (f) loosely coupled with gravity 
by means of a coupling function 6'(</)). 

Field equations in such scalar-tensor theory are derived from the action [20] 



where the first, second and third terms in the right side of Eq. (3.1.1) are the Lagrangian 
densities of gravitational field, scalar field and matter respectively, g = det [5^0/3] < is 
the determinant of the metric tensor gap, R is the Ricci scalar, indicates dependence 
of the matter Lagrangian C on matter fields, and ^(0) is the coupling function which is 
kept arbitrary. This makes the class of the theories we are working with to be sufficiently 
large. For the sake of simplicity we postulate that the self-coupling potential of the scalar 
field is identically zero so that the scalar field does not interact with itself. This is because 
we do not expect that this potential can lead to measurable relativistic effects within the 
boundaries of the solar system. However, this potential can be important in the case of a 



^ Experimental testing of the Lorentz-invariance of the gravity field equations (that is Einstein's 
principle of relativity for gravitational field) requires introducing more parameters [10,20]. We 
assume in this paper that the Lorentz-invariance is not violated. 




(3.1.1) 



strong gravitational field and its inclusion to the theory can lead to interesting physical 
consequences [65]. 

Equations of gravitational field are obtained by variation of the action (3.1.1) with respect 
to Qaf} and it's spatial derivatives. It yields 



R,u = 1^ (t,. - \g,,T^ + e{<P)^ + ^ (0;,. + lg,.n,<P ) , (3.1.2) 



where 



is the scalar Laplace-Beltrami operator and T^,^ is the stress-energy-momentum tensor of 
matter comprising the N-body system. It is defined by equation [66] 

The field equation for the scalar field is obtained by variation of the action (3.1.1) with 
respect to and it's spatial derivatives. After making use of use of the contracted form of 
Eq. (3.1.2) it yields 

1 /Stt „ , , de^ 



In what follows, we shall also utilize another version of the Einstein equations (3.1.2) which 
is obtained after conformal transformation of the metric tensor 

9,u = ^g,u , ^" = ^9'". (3.1.6) 

Here (po denotes the background value of the scalar field which will be introduced in (3.3.1). 
It is worth noting that the determinant g of the conformal metric tensor relates to the 
determinant g of the metric g^^i, a.s g = (0/0o)^5'. Conformal transformation of the metric 
tensor leads to the conformal transformation of the Christoffel symbols and the Ricci tensor. 
Denoting the conformal Ricci tensor by R^i, one can reduce the field equations (3.1.2) to a 
simpler form 

^ V ~ 2^^" J — 202 — '^■'^ "^'^ ■ ^^■^■^> 

The metric tensor g^i, is called the physical (Jordan-Fierz-frame) metric [20] because it is 
used in real measurements of time intervals and space distances. The conformal metric g^p is 
called the Einstein-frame metric. Its main advantage is that this metric is in many technical 
aspects more convenient for doing calculations than the Jordan-Fierz frame metric. Indeed, if 
the last (quadratic with respect to the scalar field) term in Eq. (3.1.7) was omitted, it would 
make them look similar to the Einstein equations of general relativity. Nevertheless, we 
prefer to construct the parametrized post-Newtonian theory of reference frames for N-body 
problem in terms of the Jordan-Fierz-frame metric in order to avoid unnecessary conformal 
transformation to convert results of our calculations to physically meaningful form. 



3.2 The Tensor of Energy- Momentum 



In order to find the gravitational field and determine the motion of the bodies comprising 
the N-body system one needs: 

(1) to specify a model of matter composing of the N-body system, 

(2) to specify the gauge condition on the metric tensor ga/s, 

(3) to simplify (reduce) the field equations by making use of the chosen gauge, 

(4) to solve the reduced field equations, 

(5) to derive equations of motion of the bodies by making use of the solutions of the field 
equations. 

This program will be completed in the present paper for the case of an isolated system 
of N bodies moving slowly and having weak gravitational field. In principle, the formalism 
which will be developed in the present paper allows us to treat N-body systems consisting 
of black holes, neutron stars, or other compact relativistic bodies if the strong field zones 
are excluded and the appropriate matching of the strong-field and weak-field zones is done 
[67]. This problem will be considered somewhere else. The most important example of the 
weak-field and slow-motion N-body system represents our solar system and one can keep 
this example in mind for future practical applications of the PPN formalism developed in 
the present paper. 

We assume that the N-body system is isolated which means that we neglect the tidal influence 
of other matter in our galaxy on this system. Thus, the space-time very far away outside of 
the system is considered as asymptotically-fiat so that the barycenter of the N-body system 
is either at rest or moves with respect to the asymptotically fiat space along a straight line 
with a constant velocity. We assume that the matter comprising the bodies of the N-body 
system is described by the energy-momentum tensor with some equation of state which we 
do not specify. Following Fock [48] ^ we define the energy-momentum tensor as 

c^T>"' = p[c^ + n)u>'u'' + 71'"' , (3.2.1) 

where p and 11 are the density and the specific internal energy of matter in the co-moving 
frame, = dx°'/cdT is the dimensionless 4- velocity of the matter with r being the proper 
time along the world lines of matter, and the anisotropic tensor of stresses defined in 

such a way that it is orthogonal to the 4-velocity 

u^Tiap = . (3.2.2) 

Original PPN formalism treats the matter of the N-body system as a perfect fluid for which 

where p is an isotropic pressure [20]. Perfect-fluid approximation is not sufficient in the 
Newtonian theory of motion of the solar system bodies where tidal phenomena and dissi- 
pative forces play essential role [68]. It is also inappropriate for consideration of last stages 
of coalescing binary systems for which a full relativistic theory of tidal deformations must 



^ See also [49] which develops similar ideas. 



be worked out. For this reason we abandon the perfect-fluid approximation and incorpo- 
rate the anisotropic stresses to the PPN formahsm. General relativistic consideration of the 
anisotropic stresses has been done in papers [69,70,71,72]. 

Conservation of the energy- momentum tensor T'^'^.^, = leads to the equation of continuity 

and the second law of thermodynamics that is expressed as a differential relationship between 
the specific internal energy and the stress tensor 

pu"n,^ + 7r"^u„.^ = . (3.2.5) 

These equations define the structure of the tensor of energy-momentum and will be employed 
later for solving the field equations and derivation of the equations of motion of the bodies. 

3.3 Basic Principles of the Post-Newtonian Approximations 

Field equations (3.1.2) and (3.1.5) all together represent a system of eleventh non-linear 
differential equations in partial derivatives and one has to find their solutions for the case of 
an N-body system. This problem is complicated and can be solved only by making use of ap- 
proximation methods. Two basic methods are known as the post-Minkowskian (see [67,73,83] 
and references therein) and the post-Newtonian (see [67] and references therein) approxima- 
tion schemes. The post-Newtonian approximation (PNA) scheme deals with slowly moving 
bodies having weak gravitational field which makes it very appropriate for constructing the 
theory of the relativistic reference frames in the solar system than the post-Minkowskian ap- 
proximation (PMA) scheme. This is because PMA does not use the slow-motion assumption 
and solves the gravity field equations in terms of retarded gravitational potentials which are 
not very convenient for description of relativistic celestial mechanics of isolated systems. For 
this reason, we shall mostly use the PNA scheme in this paper though some elements of the 
post-Minkowskian approximation (PMA) scheme will be used for definition of the multipole 
moments of the gravitational field. 

Small parameters in the PNA scheme are ei ~ f /c and e2 ~ U/c^, where w is a characteristic 
velocity of motion of matter, c is the ultimate speed (which is numerically equal to the speed 
of light in vacuum), and U is the Newtonian gravitational potential. Due to validity of the 
virial theorem for self-gravitating isolated systems one has e2 ~ and, hence, only one small 
parameter ei can be used. For the sake of simplicity we introduce parameter e = 1/c and 
consider it formally as a primary parameter of the PNA scheme so, for example, ei = ev, 
€2 = e^U, etc. 

One assumes that the scalar field can be expanded in power series around its background 
value 00, that is 

= 0o(l + C) , (3.3.1) 
where ( is dimensionless perturbation of the scalar field around its background value. The 
background value 0o of the scalar field can depend on time due to cosmological evolution of 
the universe but, according to Damour and Nordtvedt [74], such time-dependence is expected 
to be rather insignificant due to the presumably rapid decay of the scalar field in the course 



of cosmological evolution following immediately after the Big Bang. According to theoretical 
expectations [74] and experimental data [3], [4], [20] the variable part ( of the scalar field 
must have a very small magnitude so that we can expand all quantities depending on the 
scalar field in Taylor series using ^ as a small parameter. In particular, decomposition of the 
coupling function 6[(f)) can be written as 

9{(P)=u + u;' C + 0{C^) , (3.3.2) 

where u = 9{(l)o), u' = {d9/d()^^^^^, and we assume that ( approaches zero as the distance 
from the N-body system grows to infinity. 

Accounting for the decomposition of the scalar field and Eq. (3.1.5) the gravity field equations 
(3.1.2) assume the following form 



C)C2 



' UJ + l uj'C 



+ V {lu + l){2iu + 3) 



(3.3.3) 



2^^" 2u; + 3 (1 + C)' l + C 



where ^ = l/0o is the bare value of the universal gravitational constant and we have taken 
into account only linear and quadratic terms of the scalar field which is sufficient for devel- 
oping the post-Newtonian parametrized theory of the reference frames in the solar system. 

We look for solutions of the field equations in the form of a Taylor expansion of the metric 
tensor and the scalar field with respect to the parameter e such that 



(1) 2 S ^/ 4\ 

fi'a/3 = '7Q/3 + e /la/3 +e /ia/3+e hap+0[e), (3.3.4) 



or, more explicitly. 



o (2) . (4) ^ 

goo = -1 + hoo hoo +0{e^) , (3.3.5) 

(1) n (3) K 

9oi = e hoi+e"^ hoi+0{e^) , (3.3.6) 

(2) (4) 

gij = 6ij + hij hij +0{e^) , (3.3.7) 

(2) (4) 

C = e2 C C +0(e6), (3.3.8) 

(n) (") 

where hap and C, denote terms of order (n = 1,2,3...). It has been established that 
the post-Newtonian expansion of the metric tensor in general theory of relativity is non- 
analytic [67] . However, the non-analytic terms emerge in the approximations of higher post- 
Newtonian order and does not affect our results since we restrict ourselves only with the first 
post-Newtonian approximation. The first post-Newtonian approximation involves explicitly 

(2) (4) (1) (3) (2) (2) 

only terms hoo, hoo, hoi, hoi, hij and Q . In what follows we shall use simplified notations 
for the metric tensor and scalar field perturbations: 

_(2) _(4) _(1) _(3) _(2) _(2) 

N =hoo , L=hoo, Ni=hoi, =hoi , =hij , H =hkk , (3.3.9) 



and 



(2) 



(3.3.10) 



The post-Newtonian expansion of the metric tensor and scalar field introduces a correspond- 
ing expansion of the energy-momentum tensor 



Too = Too +e' Too +0(e^) , (3.3.11) 

(1) o (3) r- 

To. = eTo^+e'To. +0(e^), (3.3.12) 

O (2) A , p., , , 

T,, = Tij T^, +0(e^) , (3.3.13) 

(n) 

where again Tap (n = 1,2,3...) denote terms of order e". In the first post-Newtonian 

(0) (2) (1) (2) 

approximation we need only Too, Too, Toi and Tij which are given by the following equations 



(0) 

Too = P*, (3.3.14) 



To^ = - P* {v' + N') , (3.3.15) 

Tij = p* (v' + N') (v^ + N^) + tt'^ , (3.3.16) 

Too = P* (^y - v'Nk - ^N'^Nk + H - iV - y ^ . (3.3.17) 
Here we have used the invariant density [48] 

p* = ^u^p = p + e^p Qi/ + ^v^ + ^NkN' + v'^Nk^ , (3.3.18) 

that replaces density p and is more convenient in calculations because it satisfies the exact 
Newtonian-like equation of continuity (3.2.4) which can be recast to [20,48] 

cp*o+(p*^'),i = 0, (3.3.19) 
where v = (w*) is the 3-dimensional velocity of matter such that = cu^ /vP. 



3.4 The Gauge Conditions and the Residual Gauge Freedom 



The gauge conditions imposed on the components of the metric tensor had been proposed 
by Nutku and are chosen as follows [51,52] 




(3.4.1) 



By making use of the conformal metric tensor one can recast Eq. (3.4.1) to the same form 
as the de Bonder (or harmonic) gauge conditions in general relativity [48,49] 



(3.4.2) 



In what follows, we shall use a more convenient form of Eq. (3.4.1) written as 



(^'^T" = Inf , (3.4.3) 



^0, 



so the Laplace-Beltrami operator (3.1.3) assumes the form 

(3.4.4) 




\dx^^dx'" 

Dependence of this operator on the scalar field is a property of the adopted gauge condition. 

Any function F(x") satisfying the homogeneous Laplace-Beltrami equation, \Z\g F{x°') = 0, 
is called harmonic. Notice that DgO;" = —(In </>)'" 7^ 0, so the coordinates defined by the 
gauge conditions (3.4.3) are not harmonic functions. Therefore, we shall call the coordinate 
systems singled out by the Nutku conditions (3.4.1) as quasi-harmonic. They have many 
properties similar to the harmonic coordinates in general relativity. The choice of the quasi- 
harmonic coordinates for constructing theory of the relativistic reference frames in the scalar- 
tensor theory of gravity is justified by the following three factors: (1) the quasi-harmonic 
coordinates become harmonic when the scalar field (p 4'o, (2) the harmonic coordinates are 
used in the resolutions of the lAU 2000 [15] on relativistic reference frames, (3) the condition 
(3.4.1) significantly simplifies the field equations and makes it easier to find their solutions. 
One could use, of course, the harmonic coordinates too as it has been done, for example, by 
Klioner and Soffel [1]. They are defined by the condition = but as we found the field 

equations and the space-time transformations in these coordinates look more complicated in 
contrast to the quasi-harmonic coordinates defined by the Nutku conditions (3.4.1). 

Post-Newtonian expansion of the gauge conditions (3.4.3) yields 

Nk,k = 0, (3.4.5) 

I r , AT \ U AT ATK \ _ ( , ■p.T , Tj at atA: 



^ ^ . N + H-NkN') =- i^— + N + H-NkN''] (3.4.6) 
2 \uj -\- Z / z \uj + Z 



1 f 2ip 



N + H- NkN" ) = AT,, + Hik,k - cN,,o + {N,,k - 2Nk,^) . (3.4.7) 



2 + 2 

It is worth noting that in the first PNA the gauge- condition Eqs. (3.4.5) - (3.4.7) do not 

(4) 

restrict the metric tensor component hoo= L. 

Gauge equations (3.4.5) - (3.4.7) do not fix the coordinate system uniquely. Indeed, if one 
changes coordinates 

— > u;" = w'^ (x") , (3.4.8) 

the gauge condition (3.4.3) demands only that the new coordinates must satisfy the 
homogeneous wave equation 

which have an infinite set of non-trivial solutions. 



Eq. (3.4.9) describe the residual gauge freedom existing in the class of the quasi-harmonic 
coordinate systems restricted by the Nutku gauge conditions (3.4.3). This residual gauge 
freedom in the scalar-tensor theory is described by the same equation (3.4.9) as in the case 
of the harmonic coordinates in general relativity. We shall discuss this gauge freedom and 
its applicability to the theory of astronomical reference frames in more detail in section 7. 



3.5 The Reduced Field Equations 



Reduced field equations for the scalar field and the metric tensor are obtained in the first 
post-Newtonian approximation from Eqs. (3.1.5) and (3.3.3) after making use of the post- 
Newtonian expansions, given by Eqs. (3.3.5) - (3.3.13). Taking into account the gauge con- 
ditions (3.4.5) - (3.4.7) significantly simplifies the field equations. 

The scalar-tensor theory of gravity with variable coupling function 9{(j)) has two additional 
(constant) parameters uj and tu' with respect to general relativity. They are related to the 
standard PPN parameters 7 and (3 as follows [20] 



7 = 7H = ^lj, 3.5.1 

We draw attention of the reader that in the book [20] (equation (5.36)), parameter A = /? — 1 
is introduced as A = ct;'(2co'+3)~^(2co'+4)~^. The difference with our definition (3.5.2) given in 
the present paper arises due to different definitions of the derivative of the coupling function 
6 with respect to the scalar field, that is (ci;')wiii = 0o^(i^')this paper where 0o is the asymptotic 
value of the scalar field ^ . All other parameters of the standard PPN formalism describing 
possible deviations from general relativity are identically equal to zero [20]. General relativity 
is obtained as a limiting case of the scalar-tensor theory when parameters 7 = /5 = 1. In 
order to obtain this limit parameter uj must go to infinity with u' growing slower than u^. If 
this was not the case one could get lim^^oo 7 = 1 but lim^^oo /5 7^ 1 which is not a general 
relativistic limit. 

One can note also that the scalar field perturbation (3.3.10) is expressed in terms of 7 as 

C = (1 - 1)V ■ (3.5.3) 

As it was established by previous researchers (see, for instance, [20]) the background scalar 
field 00 and the parameter of coupling uj determine the observed numerical value of the 
universal gravitational constant 

'-^) 

where Q = l/4>o. Had the background value of the scalar field driven by cosmological evolu- 
tion, the measured value of the universal gravitational constant would depend on time and 
one could hope to detect it experimentally. The best upper limit on time variability of G is 
imposed by lunar laser ranging (LLR) as \G/G\ < 0.5 x 10"^^ yr~^ [3]. 



^ We thank CM. Will for pointing out this difference to us. 



After making use of the definition of the tensor of energy-momentum, Eqs. (3.3.14) - (3.3.17), 
and that of the PPN parameters, Eqs. (3.5.1) - (3.5.4), one obtains the final form of the 
reduced field equations: 



□ = -AnGp* 



/V2 

L + — + 2{(3-l)^' 



(3.5.5) 
(3.5.6) 



-SvrGp* + - (iV,,fe - iVfc,,) (iV,,fc - iV,,,) + e^l H^,,^N^,^ 



-SttGp* 
□ iVi = 0, 



( 7 + 1 ) 2 + n + 7 ^ - 1^ - ( 2/3 - 7 - 1 ) <^ 
/ p* b 



nLi = 87rGp* {l + -fy + N' 



i,Ok 



(3.5.7) 
(3.5.8) 
(3.5.9) 



where □ = 'q'^'^d^du is the D'Alembert (wave) operator of the Minkowski space-time, and 
H^ij^ = Hij—6ijH/3 is the symmetric trace-free (STF) part of the spatial components of the 
metric tensor. In these field equations we keep the terms quadratic on Ni, but cubic terms 
and ones proportional to the products of and perturbations of the metric are omitted. 

Equations (3.5.5) - (3.5.9) are valid in any coordinate system which is admitted by the 
residual gauge freedom defined by the gauge conditions (3.4.1). We shall study this residual 
gauge freedom in full details when constructing the global coordinates for the entire N-body 
system and the local coordinates for each of the bodies. Global coordinates in the solar system 
are identified with the barycentric reference frame and the local coordinates are associated 
with planets. The most interesting case of practical applications is the geocentric coordinate 
frame attached to Earth. 



4 Global PPN Coordinate System 



4-1 Dynamic and Kinematic Properties of the Global Coordinates 



We assume that the gravitational and scalar fields are brought about by the only one system 
comprising of N extended bodies which matter occupies a finite domain in space. Such 
an astronomical system is called isolated [48,49,75] and the solar system consisting of Sun, 
Earth, Moon, and other planets is its particular example. Astronomical systems like a galaxy, 
a globular cluster, a binary star, etc. typify other specimens of the isolated systems. A number 
of bodies in the N-body system which must be taken into account depends on the accuracy 
of astronomical observations and is determined mathematically by the magnitude of residual 
terms which one must retain in calculations to construct relativistic theory of reference frames 
being compatible with the accuracy of the observations. Since we ignore other gravitating 
bodies residing outside of the N-body system the space-time can be considered on the global 



scale as asymptotically-flat so the metric tensor gai3 at inflnity is the Minkowski metric 
r]a[3 = diag(-l, +1, +1, +1). 

In the simplest case the N-body system can be comprised of several solitary bodies, like it 
is shown in Fig 1, but in the most general case it has more complicated hierarchic structure 
which consists of a sequence of sub-systems each being comprised of Mp bodies where p is 
a serial number of the sub-system (see Fig 2). In its own turn each of the sub-systems can 
contain several sub-sub-systems, and so on. In order to describe dynamical behavior of the 
entire N-body system one needs to introduce a global 4-dimensional coordinate system. We 
denote such global coordinates = {x^, x*), where = ct is time coordinate and x^ = x are 
spatial coordinates. Adequate description of dynamical behavior of the sub-systems of bodies 
and /or solitary celestial bodies requires introducing a set of local coordinates attached to 
each of the sub-systems (or a body) under consideration. Hence, a hierarchic structure of the 
coordinate charts in the N-body system repeats that of the N-body system itself and is fully 
compatible with mathematical notion of differentiable manifold [76,77,78]. We shall discuss 
local coordinates later in section 6. 

Let us deflne the metric tensor perturbation with respect to the Minkowski metric (c.f. Eq. 
(3.3.4)) 

Kp{t,x) = ga(3{t,x) - r]ap . (4.1.1) 
We demand that quantities rhap and r'^ha^^^ are bounded, and 

Jim Kp{t,x)=^, (4.1.2) 

t+r/c=const. 

where r = \x\. Additional boundary condition must be imposed on the derivatives of the 
metric tensor to prevent appearance of non-physical radiative solutions associated with the 
advanced wave potentials [48]. It is written as 

hm [{rhap) ,r + {rhap) ,o] = . (4.1.3) 

i+r/c=const. 

Eq. (4.1.3) is known as a " no- incoming-radiation" boundary condition [48,79]. In the case of 
an isolated astronomical system this condition singles out a causal solution of the D'Alembert 
wave equation depending on the retarded time t — r/c only. Similar boundary conditions are 
imposed on the perturbation of the scalar fleld deflned in Eq. (3.3.2) 

Jim C(t,=«) = 0, (4.1.4) 

t+r/c=const. 

Jim [(rC),. + (rO,o] = 0. (4.1.5) 

f+r/c=const. 

In principle, the boundary conditions (4.1.3) and (4.1.5) are not explicitly required in the 
flrst post-Newtonian approximation for solving equations (3.5.5) - (3.5.9) because the grav- 
itational potentials in this approximation are time-symmetric. However, they are convenient 
for doing calculations and are physically motivated. Therefore, we shall use the (radiative) 
boundary conditions (4.1.3) and (4.1.5) later on for giving precise deflnitions of the multipole 
moments of the gravitational fleld of the isolated astronomical system. 

The global coordinates cover the entire space-time and they set up a primary basis for 
construction of the theory of relativistic reference frames in the N-body system [45]. In what 



follows, we shall assume that the origin of the global coordinates coincides with the barycenter 
of the N-body system at any instant of time. This condition can be satisfied after choosing 
a suitable definition of the post-Newtonian dipole moment D* of the N-body system and 
equating its numerical value to zero along with its first time derivative (see section 5.5). This 
can be always done in general relativity in low orders of the post-Newtonian approximation 
scheme if one neglects the octuple and higher-order multipole gravitational radiation [80]. In 
the scalar-tensor theory of gravity one has to take into account gravitational wave emission in 
the form of scalar modes [65] but it does not affect the first post-Newtonian approximation 
which is our main concern in the present paper. There are alternative theories of gravity 
which violate the third Newton's law so that the dipole moment D* of an N-body system is 
not conserved even in the first post-Newtonian approximation [20] but we do not consider 
such extreme cases. 

We shall also assume that spatial axes of the global coordinates do not rotate in space either 
kinematically or dynamically [59]. Spatial axes of a coordinate system are called kinemati- 
cally non-rotating if their orientation is kept fixed with respect to a Minkowski coordinate 
system defined at the infinite past and at the infinite distance from the solar system . Such 
kinematically non-rotating coordinate system can be built on the stellar sky by making use 
of quasars as reference objects with accuracy better than 100 /iarcsec (see [81] and refer- 
ences therein). Quasars are uniformly distributed all over the sky and have negligibly small 
parallaxes and proper motions . Thus, kinematically non-rotating coordinate system can 
be determined only through the experimental analysis of global properties of the space-time 
manifold including its global topology. This consideration reveals that the theory of reference 
frames in N-body system based on the assumption that the space-time is asymptotically-flat 
may be corrupted by the influence of some cosmological effects. Hence, a more appropriate 
approach to the reference frames taking into account that the background space-time is the 
Friedmann-Robertson- Walker universe is to be developed. We have done a constructive work 
in this direction in papers [27,28] but the results of these papers are still to be matched with 
the post-Newtonian approximations. 

Dynamically non-rotating coordinate system is defined by the condition that equations of 
motion of test particles moving with respect to these coordinates do not have any terms that 
might be interpreted as the Coriolis or centripetal forces [59]. This definition operates only 
with local properties of the space-time and does not require observations of distant celestial 
objects like stars or quasars. Dynamical definition of spatially non-rotating coordinates is 
used in construction of modern ephemerides of the solar system bodies which are based 
primarily on radar and laser ranging measurements to planets and Moon (see [21,81,82] and 
references therein). Because of the assumption that the N-body system under consideration 
is isolated, we can postulate that the global coordinates does not rotate at all in any sense. 



10 Angular velocities of dynamic and kinematic rotations of a reference frame in classic celestial 
mechanics are equal. However, they have different values already in the first post-Newtonian ap- 
proximation due to the presence of the relativistic geodetic precession caused by the orbital motion 
of the body. 

At the relativistic language the domain of the asymptotically-flat space-time is located both at 
the infinite distance and infinite past. This boundary consisting of null rays is called past null 
infinity [76]. 

Proper motion of an astronomical object in the sky is defined as its transverse motion in the 
plane of the sky being orthogonal to the line of sight of observer located at the barycenter of the 
solar system. 



4-2 The Metric Tensor and the Scalar Field in the Global Coordinates 



The metric tensor g^pit^x) is obtained by solving the field equations (3.5.5) - (3.5.9) after 
imposing the boundary conditions (4.1.2) - (4.1.4). We chose solution of the homogeneous 
equation (3.5.7) as Ni = 0. This is because TVj describes rotation of spatial axes of the 
coordinate system but we assumed in the previous section that the global coordinates are 
not rotating. It yields solution of the other field equations in the following form 



ip{t, X 

N{t,x 
L{t, X 
Li(t, X 

(ty X 



■■U{t,x), 
■■2U{t,x), 

■■ 2$(t, x) - 2/3U\t, x) - c\^oo{t, x) 
:-2(l + 7)[/,(t,a;), 
■.2-fSijU{t,x), 



where 



$(t, a;) = (7 + -)<l>i(t, x) + {l- 2p)<l>2{t, x) + <l>3(t, x) + j<^^{t, x) 



(4.2.1) 
(4.2.2) 
(4.2.3) 
(4.2.4) 
(4.2.5) 



(4.2.6) 



and the gravitational potentials U, U\ x, and (^ = 1, ■ ■ , 4) can be represented as linear 
combinations of the gravitational potentials of each body, that is 



AAA 



A 



(4.2.7) 



Herein, the gravitational potentials of body A are defined as integrals taken over the volume 
of this body 



U^^\t,x) 


= GlL1{p*} , 


(4.2.8) 


Urit,x) 


= GlLl{pV} , 


(4.2.9) 


x'^'Kt^x) 


= -Gir{p*}, 


(4.2.10) 


^i^\t,x) 


= GlLl{pV}, 


(4.2.11) 


^^f{t,x) 


= GlL1{p*?7} 


(4.2.12) 




= GVl'l{p*U} , 


(4.2.13) 






(4.2.14) 



where notation H^^if} {n = 1, 2, 3...) is used to define the volume integral 

ll'\f}{t,x) = J fit,x')\x-xTd'x' , 



(4.2.15) 



Va 



with n being an integer, and Va - the volume of integration. Potential x is determined as a 
particular solution of the inhomogeneous equation 

V\ = -2U, (4.2.16) 

with the right side defined in a whole space. Nevertheless, it proves out that its solution (see 
Eq. (4.2.10)) is spread out over volumes of the bodies only. It is worthwhile to emphasize 
that all integrals defining the metric tensor in the global coordinates are taken over the 
hypersurface of (constant) coordinate time t. Space-time transformations can change the 
time hypersurface, hence transforming the corresponding integrals. This important issue will 
be discussed in section 7. 



5 Multipolar Decomposition of the Metric Tensor and the Scalar Field in the 
Global Coordinates 

5. 1 General Description of Multipole Moments 

In what follows a set of certain parameters describing properties of gravitational and scalar 
fields and depending on integral characteristics of the N-body system will be indispensable. 
These parameters are called multipole moments. In the Newtonian approximation they are 
uniquely defined as coefficients in Taylor expansion of the Newtonian gravitational potential 
in powers of 1/i? where R = \x\ is the radial distance from the origin of a coordinate system 
to a field point. All Newtonian multipole moments can be functions of time in the most 
general astronomical situations. However, very often one assumes that mass is conserved 
and the center of mass of the system is located at the origin of the coordinate system 
under consideration. Provided that these assumptions are satisfied the monopole and dipole 
multipole moments must be constant. 

General relativistic multipolar expansion of gravitational field is in many aspects similar to 
the Newtonian multipolar decomposition. However, due to the non-linearity and tensorial 
character of gravitational interaction proper definition of relativistic multipole moments is 
much more complicated in contrast to the Newtonian theory. Furthermore, the gauge freedom 
existing in the general theory of relativity clearly indicates that any multipolar decompo- 
sition of gravitational field will be coordinate-dependent. Hence, a great care is required 
for unambiguous physical interpretation of various relativistic effects associated with certain 
multipoles . It was shown by many researchers that in general relativity the multipo- 
lar expansion of the gravitational field of an isolated gravitating system is characterized by 
only two independent sets - mass-type and current-type multipole moments. In particular, 
Thorne [62] had systematized and significantly perfected works of previous researchers 



See, for example, section 11 where we have shown how an appropriate choice of coordinate 
system allows us to eliminate a number of coordinate-dependent terms in equations of motion 
of spherically-symmetric bodies depending on the "quadrupoles" defined in the global coordinate 
system. 

For a comprehensive historical review see papers by [62] , [83] , [84] and references therein 
Some of the most important of these works are [85,86,87,88,89]. 



and defined two sets of the post-Newtonian multipole moments as follows (see Eqs. (5.32a) 
and (5.32b) from [62]) 



where numerical coefficients 



/(/_!)(/ + 9) 6/(/-l) 



2(/ + l)(2/ + 3) ' (/ + l)(2/ + 3)' ^ ' 



rio- 2/(/-l) ,o_ (/-!)(/ + 4) 

(/ + l)(2/ + 3) ' 2(/ + 2)(2/ + 3) ' ^ ^ ■ ^ 

and the multipolar integer-valued index / runs from to infinity. In these expressions r"^ is 
the effective stress-energy tensor evaluated at post-Newtonian order in the post-Newtonian 
harmonic gauge [62] 

r°/3 = fl + Ae^u) T"^ + 7-^77 A^"^ , (5.1.6) 
^ ^ 167rG 

where 
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(5.1.7) 
(5.1.8) 
(5.1.9) 



and [/, W are gravitational potentials of the isolated astronomical system defined in Eqs. 
(4.2.7). Thorne [62] systematically neglected all surface terms in solution of the boundary- 
value problem of gravitational field equations. However, the effective stress-energy tensor 
falls off as distance from the isolated system grows as For this reason, the multipole 

moments defined in Eqs. (5.1.1), (5.1.2) are to be formally divergent. This divergency can 
be completely eliminated if one makes use of more rigorous mathematical technique devel- 
oped by Blanchet and Damour [63] for the mass multipole moments and used later on by 
Damour and Iyer [91] to define the spin multipoles. This technique is based on the theory of 
distributions [90] and consists in the replacement in Eqs. (5.1.1), (5.1.2) of the stress-energy 
pseudo-tensor r"^ defined in the entire space with the effective source r"^ which has a com- 
pact support inside the region occupied by matter of the isolated system [63,91]. Blanchet 
and Damour proved [63] that formal integration by parts of the integrands of Thome's mul- 
tipole moments (5.1.1), (5.1.2) with subsequent discarding of all surface terms recovers the 
multipole moments derived by Blanchet and Damour by making use of the compact-support 
effective source r"^. It effectively demonstrates that Thome's post-Newtonian multipole 
moments are physically (and computationally) meaningful provided that one takes care and 
operates only with compact-support terms in the integrands of Eqs. (5.1.1), (5.1.2) after their 



rearrangement with the proper use of integration by parts of the non-hnear source of gravi- 
tational field iV"^ given by Eqs. (5.1.7)-(5.1.9). This transformation was done by Blanchet 
and Damour [63] who extracted the non-divergent core of Thome's multipole moments. We 
shall use their results in this paper. 

In the scalar-tensor theory of gravity the multipolar series gets more involved because of 
the presence of the scalar field. This brings about an additional set of multipole moments 
which are intimately related with the multipolar decomposition of the scalar field outside 
of the gravitating system. We emphasize that definition of the multipole moments in the 
scalar-tensor theory of gravity depends not only on the choice of the gauge conditions but 
also on the freedom of conformal transformation of the metric tensor as was pointed out 
by Damour and Esposito-Farese [65] who also derived (in global coordinates) the set of 
multipole moments for an isolated astronomical system in two-parametric class of scalar- 
tensor theories of gravity. In this and subsequent sections we shall study the problem of 
the multipolar decomposition of gravitational and scalar fields both of the whole N-body 
system and of each body comprising the system in the framework of the scalar-tensor theory 
of gravity under discussion. In this endeavor we shall follow the line of study outlined and 
elucidated in [20,62,63,65]. The multipole moments under discussion will include the sets 
of active, conformal and scalar multipole moments. These three sets are constrained by 
one identity (see Eq. (5.4.7)). Hence, only two of the sets are algebraically (and physically) 
independent. The multipole moments we shall work with will be defined in different reference 
frames associated both with an isolated astronomical system and with a single body (or sub- 
system of the bodies) comprising the isolated system. We call all these post-Newtonian 
moments as Thorne-Blanchet-Damour multipoles after the names of the researchers who 
strongly stimulated and structured this field by putting it on firm physical and rigorous 
mathematical bases. Let us now consider the multipole moments of the scalar-tensor theory 
of gravity in more detail. 



5.2 Thorne-Blanchet-Damour Active Multipole Moments 



Let us introduce the metric tensor potentials 
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L + — + 2{(3-l)^' 



(5.2.1) 
(5.2.2) 



which enter goo{t, x) and goi{t, x) components of the metric tensor respectively. Furthermore, 
throughout this chapter we shall put A^j = and assume that the spatial metric component 
Hij is isotropic, that is H^ij^ = 0. Then, the field equations for these potentials follow from 
Eqs. (3.5.6), (3.5.8) and read 



(5.2.3) 
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where we have introduced the active mass density 



In 9 ^ 

7 + Tjh' + n + 7 — 
/ p* 



kk 



H 
6" 



(2/3-7- IV 



(5.2.5) 



and the active current mass density 



cr' = v" 



(5.2.6) 



It is worthwhile to observe that in the global coordinates one has H = Q^U {t, x) and 
(p[t,x) = U(t,x). Hence, the expression (5.2.5) for the active mass density in these coordi- 
nates is simplified and reduced to 
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Solutions of Eqs. (5.2.3) and (5.2.4) are retarded wave potentials [66] determined up to the 
solution of a homogeneous wave equation and satisfying the boundary conditions (4.1.2) - 
(4.1.3). Taking into account that potentials V and Vi are in fact components of the metric 
tensor, solutions of Eqs. (5.2.3) and (5.2.4) can be written down as 
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(5.2.8) 
(5.2.9) 



where V designates a domain of integration going over entire space, and the gauge functions 
and are solutions of the homogeneous wave equation. We notice that because the 
densities a and cr* vanish outside the bodies the integration in Eqs. (5.2.8) and (5.2.9) is 
performed only over the volume occupied by matter of the bodies. 



We take a special choice of the gauge functions as proposed in [63] (the only difference is 
the factor 2(1 + 7) instead of 4 in [63], coming from the field equation for goi component), 
namely 
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Such gauge transformation preserves the gauge conditions (3.4.1) and also does not change 
the post-Newtonian form of the scalar multipole moments which will be discussed in the 
next section. Then one can show that potentials V and can be expanded outside of the 
N-body system in a multipolar series as follows [63] 
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where overdot denotes differentiation with respect to time t. Eqs. (5.2.12) and (5.2.13) define 
the active Thorne-Blanchet-Damour mass multipole moments, J^, and the spin moments, 
Sl, which can be expressed in the first PNA in terms of integrals over the N-body system's 
matter as follows 
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(5.2.16) 



As one can see the mass and current multipole moments of the scalar-tensor theory define 
gravitational field of the metric tensor outside of the N-body system as well as in general 
relativity [62,63]. When /3 = 7 = 1 these multipole moments coincide with their general 
relativistic expressions [63]. However, in order to complete the multipole decomposition of 
gravitational field in the scalar-tensor theory one needs to obtain a multipolar expansion of 
the scalar field as well. 



5.3 Thorne-Blanchet-Damour Scalar Multipole Moments 



In order to find out the post-Newtonian definitions of the multipole moments of the scalar 
field we again shall use the same technique as in [62,63]. We take Eq. (3.1.5) and write it 
down with the post-Newtonian accuracy by making use of a new (scalar) potential 



r/-(7-l)(7-2) 



(5.3.1) 



Then, Eq. (3.1.5) assumes the form 



UV = -AirGa , 



(5.3.2) 



where the conventional notation 77 = 4/5 — 7 — 3 for the Nordtvedt parameter [20] has been 
used and the scalar mass density a is defined as 
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We can easily check out that in the global coordinates, where H = 6'-fU{t, x) and ip{t, x) = 
U{t,x), the scalar mass density is simplified and is given by 



Solution of Eq. (5.3.2) is the retarded scalar potential 
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Multipolar decomposition of the potential (5.3.5) has the same form as in Eq. (5.2.12) with 
the scalar mass multipole moments defined as integrals over a volume of matter of the N-body 
system 
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We conclude that in the scalar-tensor theory of gravity the multipolar decomposition of 
gravitational field requires introduction of three sets of multipole moments - the active mass 
moments II, the scalar mass moments 1^, and the spin moments Sl- Neither the active 
nor the scalar mass multipole moments alone lead to the laws of conservation of energy, 
linear momentum, etc. of an isolated system; only their linear combination does. This linear 
combination of the multipole moments can be derived after making conformal transformation 
of the metric tensor, solving the Einstein equations for the conformal metric, and finding its 
multipolar decomposition in the similar way as it was done in section 5.2. 



5.4 Thome- Blanchet-D amour Conformal Multipole Moments 



Let us now define the conformal metric potential 
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The conformal field equations (3.1.7) in the quasi-harmonic gauge of Nutku (3.4.2) yield 



□ 1/=-47rG'a, 
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where we have introduced a conformal mass density 
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which has been calculated directly from Eq. (3.1.7) by making use of the definition of the 
conformal metric (3.1.6) and the post-Newtonian expansions of corresponding quantities 
described in section 3.3. Remembering that in the global coordinates H = 6'yU{t,x) and 
ip{t, x) = U{t, x) one can simplify expression for the conformal mass density which assumes 
the form 
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This equation coincides precisely with the post-Newtonian mass density as it is defined in 
general relativity (see [20,48] and [65] for more detail). The con/orma/ current density a* is 
defined in the approximation under consideration by the same equation as Eq. (5.2.6), that 
is 5"* = (T*. The field equation for the conformal vector potential V'^ has the form (5.2.4), 
therefore = V\ 

Solution of Eq. (5.4.2) gives the retarded conformal potential 

V{t, x) = G I ^3^, _ ^5 4 5) 

Jv \x — x'\ 

Multipolar expansion of conformal potentials V and is done in the same way as it was 
done previously in section 5.2. It turns out that the conformal spin moments coincide with 
the active spin moments (5.2.16), and the expansion of the potential V{t,x) acquires the 
same form as that given in Eq. (5.2.12) but with all active multipole moments replaced with 
the conformal multipoles, /<l>, defined as follows 



i<L>{t) = /^a(t,x')x'<^>rfV + ^^^f^ [^—^jj{t,x')x'<'^>x'^d\' (5.4.6) 

_g(2[+_l)^ I a\t,x')x'<'^>d'x'\ . 
l + l dtJv ^ ' J 

These conformal mass multipole moments coincide exactly with those introduced by Blanchet 
and Damour [63] who also proved (see appendix A in [63]) that their definition coincides pre- 
cisely (after formal discarding of all surface integrals which have no physical meaning) with 
the mass multipole moments introduced originally in the first post-Newtonian approximation 
in general relativity by Thorne [62]. 

There is a simple algebraic relationship between the three mass multipole moments, J^, II 
and II in the global frame. Specifically, one has 

I<L> = ^^^<L> + ll<L> . (5.4.7) 

We shall show later in section 6.3 that relationship (5.4.7) between the multipole moments 
obtained in the global coordinates for the case of an isolated astronomical N-body system 
preserves its form in the local coordinates for each gravitating body (a sub-system of the 
bodies) as well. 



5.5 Post-Newtonian Conservation Laws 



It is crucial for the following analysis to discuss the laws of conservation for an isolated 
astronomical system in the framework of the scalar-tensor theory of gravity. These laws will 
allow us to formulate the post-Newtonian definitions of mass, the center of mass, the linear 
and the angular momenta for the isolated system which are used in derivation of equations 
of motion of the bodies comprising the system. In order to derive the laws of conservation 
we shall employ a general relativistic approach developed in [66] and extended to the Brans- 
Dicke theory by Nutku [52]. 



To this end it is convenient to recast the field equations (3.1.2) to the form 
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where t^'^ is an analog of the Landau-Lifshitz pseudo-tensor of the gravitational field in the 
scalar-tensor theory of gravity. This pseudotensor is defined by the equation 
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(5.5.2) 



where r|^^ is the (standard) Landau-Lifshitz pseudotensor [66] expressed in terms of the 
conformal metric ^o/3 and its derivatives. 

The conservation laws are now obtained from Eq. (5.5.1) 
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(5.5.3) 



They are a direct consequence of anti-symmetry of the right side of Eq. (5.5.1) with respect 
to the upper indices u and a. In what follows, we concentrate on the laws of conservation 
in the first post-Newtonian approximation only. Hence, we neglect the energy, linear and 
angular momenta taken away from the system by gravitational waves (see [65] where this 
problem has been tackled). For this reason, the conserved mass M, the linear momentum P', 
and spin S* of the isolated gravitating N-body system are defined as 



-e'l^e^Wx, (5.5.4) 
Q°'d^x, (5.5.5) 
u;^e°'= d^x . (5.5.6) 



V 

In these definitions integration is performed over the whole space. Let us remark that the 
integrals are finite since in the first PNA 0"° and are of 0{r~^) for large r. Moreover, 
in this approximation the domain of integration can be reduced to the volume of the bodies 
comprising the system - observe that in (5.5.9) - (5.5.11) the functions under the integrals 
are compactly supported. Taking into account the asymptotic behavior of one can prove 
that the linear momentum P* can be represented as the time derivative of the function 

W = e^ I e^Vrf^x , (5.5.7) 
Jv 

which is interpreted as the integral of the center of mass. Hence, 

W{t)=rt + M.\ (5.5.8) 

where K* is a constant vector defining displacement of the barycenter of the N-body system 
from the origin of the global coordinate frame. One can chose K* = and P* = 0. In such 
case W = 0, and the center of mass of the N-body system will always coincide with the 



origin of the global reference frame. Such global reference frame is called barycentric. It is 
used in description of ephemerides of the solar system bodies, navigation of spacecrafts in 
deep space and reduction of astronomical observations of various types. 



Direct calculations of the pseudotensor (5.5.2) with subsequent comparison with the con- 
formal multipole moments (5.4.6) reveal that for the isolated system the post-Newtonian 
conserved quantities are 
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where by definition 



W\t,x) = G 
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and the integration is performed over the hypersurface of constant global coordinate time t. 
It is evident from Eqs. (5.5.11) and (5.5.12) that it is the con/orma/ moments, / and which 
define the conserved mass M and linear momentum W of the N-body system. The active 
monopole and dipole moments defined by Eq. (5.2.14) for / = 0, 1 are not consistent with 
the laws of conservation and, hence, can not serve to define the conserved quantities. We fix 
position of the center of mass (barycenter) of the N-body system in the global coordinates 
by equating conformal dipole moment of the system to zero, that is /* = 0. 

Now we are prepared to begin construction of a local coordinate system in the vicinity of a 
gravitating body or a sub-system of bodies which are members of the entire N-body system. 
For concreteness and for the sake of simplicity we shall focus on the construction of the local 
coordinate system around one body (Earth, planet, etc.). 



6 Local PPN Coordinate System 



6. 1 Dynamic and Kinematic Properties of the Local Coordinates 



Local coordinate system (local coordinates) is constructed in the vicinity of each body com- 
prising the N-body system . Thus, in principle, N local coordinate systems must be 
introduced in addition to one global coordinate system (see Fig. 1). In the case of the 
N-body system which is divided on sub-systems of bodies the number of the local coordi- 
nates increases in accordance with the underlying hierarchic structure of the N-body system 

Precise definition of body's center of mass will be given in subsequent sections along with deriva- 
tion of its equations of motion. 



(see Fig. 2). The principles of construction of the local coordinates are the same for any 
weakly gravitating body (a sub-system of bodies). For this reason, it is sufficient to work out 
description of only one local coordinate system, to" = (cm, w), as the other local coordinate 
charts must have a similar structure [58]. For practical applications in the solar system the 
most important local coordinates are associated with the Earth and they are called geocen- 
tric coordinates. Local coordinates are not asymptotically Minkowskian far away from the 
body because the gravitational field of the body under consideration must smoothly match 
with the gravitational field of external bodies. 

We assume that each body consists of matter which admits continuous distribution of mass 
density, anisotropic stresses and internal velocity field. If one had "turned off" gravitational 
field of all external bodies (Moon, Sun, planets) but the gravitational field of the body under 
consideration (Earth), it would be described by a set of the (internal) multipole moments 
defined by equations given in previous section. However, we can not neglect gravitational field 
of the external bodies if one wants to take into account classic [76,92] and relativistic effects 
associated with tides [72,93]. The tidal deformation of the body will be comprehensively large, 
for example, at the latest stage of coalescence of neutron stars in binary systems emitting 
gravitational radiation and is to be taken into account in calculations of the templates 
of gravitational waves being emitted by such systems. We also know that already in the 
Newtonian limit this external gravitational field reveals itself in the vicinity of the Earth 
as a classical tidal force [76]. Gravitational potential of the tidal force is represented as a 
Taylor series with respect to the local geocentric coordinates with time-dependent coefficients 
which are called external (tidal) multipole moments [94]. This series usually starts in the 
Newtonian approximation from the second order (quadratic) term because the monopole 
and dipole external multipole moments are not physically associated with the tidal force. In 
general relativity this monopole- dipole effacing property of the external gravitational field 
is retained in the post-Newtonian approximation as a consequence of the Einstein principle 
of equivalence (EEP) [20,66,76]. In particular, EEP suggests that it is possible to chose 
such (local) coordinates that all first derivatives of the metric tensor (i.e., the Christoffel 
symbols) vanish along a geodesic world line of a freely falling particle [40] . This is equivalent 
to making a suitable coordinate transformation on the space-time manifold from the global 
to local frame [77,78]. In general relativity this property of EEP is also valid for a self- 
gravitating body moving in external gravitational field of other bodies. The original proof 
was given in [94,95,96] and elaborated on later in a series of papers by other researchers 
[42,43,45,61,69,70]. 

As contrasted with general relativity the scalar-tensor theory of gravity has a scalar (helicity- 
0) component of the gravitational field which can not be eliminated by a coordinate trans- 
formation to the local frame of the body being in a free fall. This is because the scalar field 
does not change its numerical value under pointwise coordinate transformations and can not 
be eliminated if it has a non-zero value on space-time manifold. It means that scalar fields do 
not obey the principle of equivalence and the gravitational field in the scalar-tensor theory 
can not be reduced in the local coordinate system to the tidal field only. In particular, this 
was the reason why Einstein had rejected a theory of gravity based exceptionally on a scalar 
field (for more detail see [20,76]). 

This argument makes it clear that in order to incorporate a local coordinate system to the 
standard PPN formalism [16] -[20] one needs to know the nature of the fundamental fields 
(scalar, vector, tensor, spinor, etc.) present in the theory because these fields have different 



behavior under coordinate transformations. To construct a local coordinate system, solution 
of the field equations for all fundamental fields must be found directly in the local coordinate 
system. Then, this solution must be matched to the solution of the same equations in the 
global coordinates and the transformation laws of the additional fields must be used along 
with the transformation law of the metric tensor in order to find relativistic space-time 
transformation between the global and local coordinate systems. Nonetheless, one has to 
keep in mind that the scalar field is not observed directly but is organically incorporated to 
the metric tensor which obeys to EEP. It means that the scalar field and its first derivative 
at each point of the manifold can be absorbed in the metric tensor and its first derivatives. 
Thus, the metric tensor in the origin of the local coordinates can be reduced to the Minkowski 
metric as long as the body's gravitational field is not to be a matter of concern. 

We demand that the origin of the local coordinates coincides with the body's center of mass 
at any instant of time. This requires a precise definition of the center of mass of each body 
with respect to its local coordinates. But when one takes into account the post-Newtonian 
corrections, the notion of the body's center of mass becomes ambiguous because it can be 
chosen in several different ways depending on what kind of definition of the internal dipole 
moment of the body in the multipolar expansion of the local metric tensor is chosen. We have 
proven by straightforward calculations that it is the con/orma/ dipole moment (5.5.10) which 
gives a physically correct definition of the body's center of mass because only this moment 
allows to derive equations of translational motion of the body which does not contain self- 
accelerated terms violating the Newton's third law of action-counteraction. This property 
of the conformal moment is closely related to its conservation for an isolated system of N 
bodies as demonstrated in section 5.5. 

In general, the body (Earth) as a part of the N-body system is not isolated and interacts 
gravitationally with other bodies (Moon, Sun, etc.). For this reason, the second and higher 
order time derivatives of the conformal dipole moment of the body are not equal to zero by 
themselves. It means that there is a local force exerted on the body by external gravitational 
field which prevents its linear momentum (the first time derivative of the body's dipole 
moment) to conserve. Nevertheless, it is possible to prove that all time derivatives of the 
body's dipole moment can be kept equal to zero if one chose the origin of the local coordinates 
to move not along a geodesic world line. This goal is achieved by making use of a specific 
choice of external dipole moment in the multipolar expansion of the homogeneous solution of 
the gravitational field equations (see section 9.4). The correct choice of the body's center of 
mass allows us to eliminate the ill-behaved coordinate-dependent terms in the equations of 
motion of the body and facilitates discussion of the strong equivalence principle's violation 
(the Nordtvedt effect) for extended bodies. 

We admit that the local coordinates can be, in general, dynamically rotating. It means that 
translational equations of motion of a test particle written down in the local coordinates can 
include the Coriolis and centrifugal forces. If one excludes the dynamical rotation of the local 
coordinates, their spatial axes will slowly rotate in the kinematic sense with respect to the 
spatial axes of the global coordinates [13,14]. This effect is called a geodetic precession and it 
obeys the law of parallel transport of vectors on curved space-time manifold [76]. Nowadays, 
the lAU recommends to use a kinematically-nonrotating geocentric system which spatial 
axes are anchored to distant quasars used as reference points of the international celestial 
reference system (ICRS) (see [97] for more detail). The metric tensor of the kinematically- 
nonrotating geocentric coordinates has an external dipole moment in gQi{u,w) component 



of the geocentric metric tensor describing the dynamical rotation of the spatial axes of the 
geocentric coordinates. This term would be zero if the geocentric coordinates were chosen 
to be dynamically nonrotating. The angular velocity of the dynamical rotation is equal to 
that of the geodetic precession and is fixed by the corresponding lAU resolution. At this step 
of development of our formalism we shall not specify the angular velocity of the dynamical 
rotation in order to keep the formalism as general as possible. 

6.2 The Metric Tensor and the Scalar Field in the Local Coordinates 

We denote the local (for example, geocentric) coordinates by = {w^, w^) = {cu, w'') where 
u stands for the local coordinate time. All quantities related to the (central) body around 
which the local coordinate frame is constructed will be labelled by subindex B standing for 
"body". We are looking for the solution of the field equations (3.5.5) - (3.5.6) inside a world 
tube containing the world line of the body's center of mass and spreading up to the nearest 
external body, so that the only source of matter inside the region covered by the local frame 
is the matter of the central body. Thus, the right side of equations (3.5.5) - (3.5.9) contains 
the energy-momentum tensor of the body's matter only. Spatial domain of applicability of 
the local coordinates can be extended after finding the space-time transformation from local 
to global coordinates [61]. 

Solution of the differential equations (3.5.5) - (3.5.6) is a linear combination of general solu- 
tion of the homogeneous equation and a particular solution of the inhomogeneous equation. 
For example, solution for a scalar field in the local coordinates is written as 

if{u, w) = (^'""(u, w) + (f^'^^u, w) , (6.2.1) 

whereas the metric tensor, g^y{u, w) = rj^y + h^y{u, w), is given in the form 

V(n, w) = h';j (n, w) + h%\u, w) + /i|f,'(^, ^) , (6-2.2) 

where terms with sub-index B refer to the central body (Earth) and describe the (internal) 
solution of the inhomogeneous equations, terms with sub-index E refer to the external bodies 
(Moon, Sun, etc.) and describe the (external) solution of the homogeneous equations, and 
terms with sub-index C (which stands for coupling) arise because of the non-linearity of the 
gravity field equations for the metric tensor. One notices that in the first post-Newtonian 
approximation the coupling terms appear only in goo{u, w) component of the metric tensor. 
We do not impose any other specific limitations on the structure of the metric tensor in local 
coordinates. All information about its structure can be obtained from the solution of the 
field equations (3.5.5) - (3.5.9). We draw attention of the reader that we put a hat over all 
quantities referred to the local coordinates w". This is because functional dependence of one 
and the same quantity looks different in different coordinates. For example, for any scalar 
function F{x) and coordinate transformation x = x{w) one has F{x) = F[x{w)] = F{w) 
while F{w) differs from F{x) [77,78]. 

6.2.1 The Scalar Field: Internal and External Solutions 

Eq. (3.5.5) gives internal, (^'^'(n, it>), and external, 0^'^\u,w), solutions for the scalar field 
in the following form 



1=0 



(6.2.3) 
(6.2.4) 



Here Pl = Pl{u) are external STF multipole moments in the multipolar decomposition 
of the scalar field generated by the bodies which are external with respect to the central 
body. These external moments are functions of the local time u only. The internal solution 
(fi^'^^u, w) describes the scalar field which is generated by the central body only. 



6.2.2 The Metric Tensor: Internal Solution 



The boundary conditions imposed on the internal solution for the metric tensor are identical 
with those given in Eqs. (4.1.2) - (4.1.3). For this reason the internal solution for the metric 
tensor has a form which is similar with that obtained in the global coordinates where all 
quantities must be referred now only to the central body. We obtain 



iV'"'(M,io)=2f/'"'(M,io) , 



Lr{u,w) = -2{l + ^)Ur{u,w 



HlJ\u,w)=2^6i,U^^\u,w) 
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(6.2.5) 
(6.2.6) 
(6.2.7) 
(6.2.8) 



where all gravitational potentials of the central body are taken over the volume of the body's 
matter defined as a cross-section of the body's world tube with the hypersurface of constant 
local coordinate time u. Specifically, one has 



l>'^'(n, lo) = (7 + ^)<l>r(M, w) + {l- 2/5)$^''^ (m, w) 
+$f^(n,iu) +7<lf' (n,tu), 
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(6.2.16) 



the symbol z/* = dw'^/du is the velocity of the body's matter with respect to the origin of the 
local coordinates, and we have introduced a special notation 



ijfH/}(w,^)= / f{u,w')\w-wrd'w' 

Vb 



(6.2.17) 



for integrals over the body's volume. We emphasize once again that the integrand of I^^'{/}(n, w) 
is a function which is taken over the hypersurface of constant time u. 

The local metric given by Eqs. (6.2.3), (6.2.5) - (6.2.8) must obey the gauge condition (3.4.1) 
which yields 

This is the only gauge condition which can be imposed on the local metric in the first post- 
Newtonian approximation. We note that Eq. (6.2.18) is satisfied due to the validity of the 
equation of continuity (3.3.19). 



6.2.3 The Metric Tensor: External Solution 



Solution of the homogeneous field equations for the metric tensor given in this section is 
based on and extends the multipolar formalism for description of vacuum gravitational fields 
developed in [94,98]. Brief introduction to this formalism is given in Appendix A. Boundary 
conditions imposed on the external solution must ensure its convergency on the world line of 
the origin of the local coordinates where w = 0. However, the external solution for the metric 
tensor diverges as the radial distance r = |it>| from the origin of the local coordinates grows. 
This is because of the gravitational field of external bodies which does not asymptotically 
vanish in the local coordinates for large r [94,98]. 

Explicit form of the external solution for the linearized metric tensor perturbation in local 
coordinates is given by 



oo -| 

N^^\u, w) = 2Y^ j^Qlw^ + ^^w^ - MWu;" , (6.2.19) 
1=0 ^■ 



iVi^'(u, w) = V^ + eip.Qpw'' , (6.2.20) 
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where f2j is the angular velocity of kinematic rotation of the local frame with respect to the 
global coordinates, V* is the velocity of the local frame with respect to the local frame moving 
along geodesic world line (see below), and either symbol "Sym(ij)" or the round brackets 
around indices denote symmetry with respect to the indices, for instance, [T^x]^^™'-*-''' = 
T(ij)i = {l/2)[TijL + Tjii]. In (6.2.19) - (6.2.22) we keep only the terms of 0(9) and {Vl?) 
which are relevant for the discussion of Newtonian geodesies. It is worth to notice that 
external solutions for the metric tensor in local coordinates contain monopole terms Q and 
Y . Term Q defines the unit of measurement of the coordinate time u at the origin of the 
local frame and Y defines the unit of measurement of spatial distances with respect to the 
international system of units. Both these terms could be equated to zero from the very 
beginning but we prefer to keep them in our equations for generality. This is because the 
lAU resolutions [14] explicitly introduce the non-zero values of Q and Y and this is why we 
are interested in the impact of these functions on the PPN theory of reference frames . 

In order to understand physical meaning of various components of the external solution for 
the metric tensor in the local coordinates it is instructive to write down the Newtonian 
equation of motion of a test particle falling freely in the field defined by the external metric. 
This equation is a geodesic world line so that after calculation of the Christoffel symbols one 
has particle's acceleration 

j2 i 

^ = - Vi - 2£,,fcfiV'= - Ei.u^'w'' + {^f5ij - (6.2.23) 

oo I 
1=2 ■ 

where = dw'^/du and we have neglected the post-Newtonian corrections. First two terms 
in the right side of this equation, Qi — Vi, describe kinematic acceleration of the particle with 



Functions Q = Y = Lc = 1.48082686741 x 10"^ ± (2 x 10^^'^ [99,100] (see also [101,102] for 
detailed theoretical review on the relativistic time scales in the solar system. 



respect to the coordinate system moving along geodesic. The third term, 2ejjfcf2V^, in the 
right side of Eq. (6.2.23) is the famous Coriohs acceleration [103] caused by motion of the 
particle and rotation of spatial axes of the local frame with angular velocity fi*. The forth 
term, eijkCl^w^, in the right side of Eq. (6.2.23) is acceleration due to the non-uniform rotation 
of the local frame. The fifth term, {Q'^6ij — QiQj) , describes a centrifugal acceleration of 
the particle. The sixth term, Q<ij>w^ , is a quadrupole tidal acceleration due to the presence 
of external gravitational field from other bodies besides the central one. Last term in the 
right side of Eq. (6.2.23) is the tidal acceleration due to the higher order multipoles of the 
external gravitational field of other bodies. It is interesting to note that the centrifugal and 
the quadrupole tidal accelerations have similar structure. The difference, however, is that the 
matrix of the centrifugal acceleration, Q'^Sij — QiQj, is not trace- free in contrast to the tidal 
matrix, Qij. However, the trace- free part of Q^6ij — QiQj can be singled out and absorbed 
to the definition of Qij. 

It is convenient to construct the external part of the metric tensor in such a way that makes 
it Minkowskian (orthogonal) at the origin of the local coordinates. This can be achieved if 
one chooses function V* = 0. This condition also allows us to give a unique interpretation 
of the dipole term Qi as equal to the inertial force per unit mass exerted on the free falling 
particle due to the accelerated motion of the local frame under consideration with respect to 
the geodesic world line. In other words, the metric tensor with V* = and Qi ^ specifies 
a local coordinate system such that its origin moves with acceleration Qi with respect to a 
geodesic world-line defined on the background space-time which is determined exclusively by 
the external part of the metric tensor. We also notice that the dipole term Zj in Eq. (6.2.21) 
is just a post-Newtonian correction to V* and would also destroy orthogonality of spatial 
axes of the local frame at its origin. Thus, in addition to the condition, V* = 0, we also 
demand, = 0. Therefore, any post-Newtonian corrections to the equations of motion of 
the origin of the local coordinates are hidden in the inertial acceleration Qi in Eq. (6.2.19). 

A set of eleven external STF multipole moments Pl, Ql, Cl, Zl, Sl, Yl, Bl, Dl, El, Gl 
(we omit < > for simplicity, i.e. Pl = P<l>, etc.) is defined on the world line of the origin 
of the local coordinates so that these multipoles are functions of the local coordinate time 
u only. Furthermore, the external multipole moments are symmetric and trace-free (STF) 
objects with respect to any of two indices, and they are transformed as tensors with respect 
to linear coordinate transformations. In what follows, we shall assume that the angular ve- 
locity of rotation of the local frame, Q\ is so small that the metric tensor component eN-'^^ 
is comparable with that e^L['^\ For this reason, we shall neglect all terms which are either 
quadratic with respect to N-"^^ or are products of N-'^^ with one of each of the other compo- 
nents, iV'^', Z/-^' or HiP. Only linear with respect to N-'^^ terms and their first derivatives 
will be retained in our calculations. 

Imposing the gauge conditions (3.4.1) on the metric tensor given by Eqs. (6.2.19) - (6.2.22) 
reveals that only 7 from 11 external multipole moments are algebraically independent. More 
specifically, the gauge condition (3.4.7) leads to the following relationship between the mo- 
ments 
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Three other relationships are obtained after accounting for Eq. (6.2.24) in the gauge condition 
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Eqs. (6.2.24) - (6.2.28) allow us to eliminate the external multipole moments Bl, El, Dl, Sl 
from the local metric so that the space-time and space-space components of the external 
metric tensor assume the form 
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Remaining multipole moments Pl, Ql, Zl, Cl, Yl, Fl, Gl and the angular velocity of 
rotation, Qi, can not be constrained by imposing the gauge conditions. However, a residual 
gauge freedom described by differential Eq. (3.4.9) allows us to find further limitations on 
the remaining 7 sets of the multipole moments which are explicitly shown in the right side 
of Eqs. (6.2.29) and (6.2.30). Exclusion of the residual degrees of freedom makes it clear 
which multipole moments are in fact physically relevant, that is can not be excluded by 
infinitesimal coordinate transformations ty" = -|- ^°'{x^). 

In order to exclude the external multipoles, which do not carry out information about gravita- 
tional degrees of freedom, we shall make use of a well-known property of the gauge-invariance 
of the (linearized) Riemann tensor under infinitesimally small coordinate transformations 
(see Box 18.2 in [76]) 

Raf3-y5{u, W) = Ra(^^s{u, w) . (6.2.31) 

Eq. (6.2.31) must be understood as invariance of the functional form of the Riemann tensor 
after making an infinitesimally small gauge transformation. Computing all components of 
the Riemann tensor, which are functions of the external metric tensor (6.2.19), (6.2.29) and 



(6.2.30) only, one finds (see Appendix B.2) that the external part of the Riemann tensor 
depends only on four sets of the external multipole moments Pl, Ql, Cl and G^. However, 
one can notice that the multipoles Cl and Gl enter Eq. (B.2) only in the form of a linear 
combination which can not be split algebraically. This means that only three sets of the 
external moments have a real physical meaning. In what follows, we shall choose Pl, Ql, 
and Cl as the primary external multipoles. Other multipole moments Yl, Zl, Fq and Cl 
can be chosen arbitrary which reflects the presence of 4 residual gauge degrees of freedom 
generated by the coordinate transformation confined by Eq. (3.4.9). 

Hereafter we assume that the angular velocity of rotation of the local frame 

a = , (6.2.32) 

which gives N-'^\u, w) = 0. This assumption greatly simplifies subsequent calculations with- 
out missing any significant physics. We have to notice, however, that rotating local coordi- 
nate systems have a great practical value for satellite geodesy and global positioning system 
(GPS) [104]. 

Various authors used the residual gauge freedom differently. We shall follow the convention 
accepted in papers [15,45,69,94] and postulate ^® that the space-space components gij{u, w) 
of the local metric tensor must form a diagonal matrix proportional to the Kronecker symbol 
6ij. It allows to simplify the non-linear term H^ij^N^ij in Eq. (3.5.6) for the time-time 
component of the metric tensor by making use of the Laplace equation for function N{t, w) 
which converts the non-linear term to that having a compact support, that is H^ij^Nij ~ 
p*H. In order to diagonalize gij{u, w) one chooses the external multipoles Fl, Cl as follows 

Fl = , (6.2.33) 

Cl = , (6.2.34) 

for alH > 0. Furthermore, we chose 

YL = QL + h- 1)Pl , (6.2.35) 

for all Yl with / > 1 but the monopole moment Y which is left arbitrary . We preserve 
some gauge freedom and do not fix the external multipoles Zl (/ > 2). 

Finally, the external metric tensor assumes the following simple form 



We draw attention of the reader that in the standard PPN formalism [20] and in the global 
coordinates the space-space components of the metric tensor are always diagonal in standard gauge. 
However, here we discuss the local coordinates which admit more general degrees of freedom than 
the global ones. Hence, our postulate is not redundant. 

Restrictions (6.2.33) ~ (6.2.35) on the external multipoles result of the requirement that the 
metric tensor must be diagonal in the entire domain of validity of the local coordinates but not in 
a single point. 
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where we still keep the (time-dependent) monopole term Q ^ 0. 



Now we can compute /iQQ'('u,iy) component of the external metric tensor up to the post- 
Newtonian order by making use of Eq. (3.5.6). Its the most general solution is determined 
up to that of a homogeneous wave equation which we shall assume to be incorporated to the 
post-Newtonian correction to the multipole moments Ql defined by Eq. (6.2.36). Hence, we 
obtain 
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It is interesting to note that summation in the first two terms in the right side of Eq. (6.2.39) 
was originally started from / = (see Eqs. (6.2.4) and (6.2.19)). However, product of a har- 
monic polynomial with the monopoles Q and P represent a homogeneous solution of the 
Laplace equation and, for this reason, can be absorbed by the Newtonian-like polynomial 
QlW^ in Eq. (6.2.36) by means of re-definition of the multipoles Ql- This is possible be- 
cause the mathematical structure of the multipoles Ql as functions of the parameters of the 
external bodies has not yet been specified. Our remark helps to realize that the influence of 
the scalar field on the external solution of the metric tensor in the local coordinates starts 
from quadratic, with respect to coordinates w, terms only. External scalar field can not be 
eliminated by pointwise coordinate transformation but it enters to the external metric tensor 
in such a way that it can be absorbed to the multipole moments Ql of the metric tensor 
gravitational field. Hence, the external multipoles Pl do not contribute linearly to the equa- 
tions of translational motion of test particles and extended bodies - only their non-linear 
combination is observable (see section 9). 



6.2.4 The Metric Tensor: The Coupling Terms 

The coupling terms in the metric tensor in local coordinates are given as a particular solution 
of the inhomogeneous equation (3.5.6) with the right side taken as a product of the internal 
and external solutions found on previous step of the approximation procedure. Solving Eq. 
(3.5.6) yields the coupling terms of the metric tensor in the local coordinates 



L(^'(n, w) = -2f/(^'(n, w) x (6.2.40) 
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This completes derivation of the metric tensor in the local coordinates. 

6.3 Multipolar Decomposition of the Body's Gravitational Field in the Local Coordinates 



The local coordinates are introduced in the vicinity of each of the gravitating body comprising 
the N-body system. We consider one of them and call it the " central" body which is indexed 
by the letter 'B'. This body, for example, can be the Earth and the local coordinates in such 
case are called the geocentric coordinates [14,15]. Gravitational field of the central body taken 
alone, that is when all other (external) bodies are ignored, is described in the local coordinates 
in terms of the metric tensor and scalar field which depend on the internal field potentials 
U'^^\ UIb)^ etc., defined in Eqs. (6.2.9) - (6.2.15). Multipolar decomposition of the 

internal metric tensor of the central body is totally equivalent to the procedure of multipolar 
decomposition of the gravitational field of the N-body in global coordinates described in 
section 11.3. However, from the point of view of precise theory the central body is not 
gravitationally isolated from the other bodies of the N-body system because it interacts with 
them gravitationally. This interaction brings about the coupling terms to the metric tensor in 
the local coordinates which can contribute to the numerical values of the body's multipole 
moments in the multipolar decomposition of the local metric tensor. The presence of the 
coupling terms introduces a post-Newtonian ambiguity of the multipolar decomposition of 
gravitational field in the local coordinates and puts a question about what specific definition 
of the multipole moments of the central body must be used in deriving translational and 
rotational equations of motion of the body in the gravitational field of N-body system. This 
problem was seemingly pointed out for the first time by Thorne and Hartle [94]. 

Solution of this problem can be found only by doing complete calculations of the equations 
of motion of the central body with taking into account all its multipoles. One has two 
possibilities: either to include or to exclude the contribution of the coupling terms to the 
multipole moments of the body and we have explored both of them. It turns out the final 
form of the equations of motion can be significantly simplified if the coupling terms are 
included to the definition of the multipole moments. In fact, if one excludes the contribution 
of the coupling terms to the body's multipoles it produces a lot of additional terms in the 
equations of motion but all these terms can be eliminated after a suitable re-definition of the 
multipole moments. By inspection one can check that the final form of such renormalized 
equations of motion coincides with that which would be obtained if one included contribution 
of the coupling terms in the local metric tensor to the definition of the multipole moments of 
the central body from the very beginning. The significantly simple form of the renormalized 
equations of motion is a direct indication that the coupling terms must be included to the 
definition of the multipole moments of the body in the local coordinates. This resolves the 
Thorne-Hartle ambiguity [94] in the definition of the multipole moments. 



Thus, the formal procedure of the multipolar decomposition of the gravitational field in 



the local coordinates is based on the same field equations (5.2.3), (5.3.2) and (5.4.2) for 
active, scalar and conformal potentials whose right sides depends on the active, scalar and 
conformal mass densities defined by Eqs. (5.2.5), (5.3.3) and (5.4.3) respectively. All these 
densities depend on the trace of the space-space component of the metric tensor, H, and the 
scalar field, if. In accordance with our procedure of definition of the multipole moments in the 
local coordinates, these functions must include the contribution of the external gravitational 
and scalar fields. In other words, computation of the mass densities in Eqs. (5.2.5), (5.3.3) 
and (5.4.3) in the local coordinates must be relied upon the trace of the metric tensor, H, 
defined by the sum of Eqs. (6.2.8) and (6.2.38), and the scalar field, ip, defined by the sum of 
Eqs. (6.2.3) and (6.2.5). Solving Eqs. (5.2.3), (5.3.2), (5.4.2) with the mass densities defined 
herewith and expanding the metric potentials in the multipolar series yields the multipole 
moments of the central body in the local coordinates. 

In section 4 we have constructed three sets of the mass multipole moments - active, scalar, 
and conformal. The same type of multipoles presents in the local coordinates as well. The 
active STF mass multipole moments of the central body are 
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where Vb denotes the volume of the central body under consideration and the active mass 
density in the body's interior is defined as 
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where f/^^' is the gravitational potential of the body given by Eq. (6.2.9). 
The scalar STF mass multipole moments of the body are defined as 
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where the scalar mass density of the body's matter is defined by 
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The conformal STF mass multipole moments of the body are 
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with the conformal mass density of the body's matter defined as 
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The conformal density does not depend on the PPN parameters (3 and 7. 
The current density is defined in the local coordinates by 



and the spin multipole moments of the body are determined by the formula 
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It is important to emphasize that the algebraic relationship (5.4.7) preserves its form for the 
set of the mass multipole moments taken for each body separately, that is 

1 + 1 ~ 1 - 

1l = —^Tl + -^l- (6.3.9) 

Validity of this relationship can be checked out by a straightforward calculation. We also 
draw attention of the reader that the hypersurface of the integration in Eqs. (6.3.1), (6.3.3), 
(6.3.5) is that of the constant local coordinate time u which does not coincide with the 
hypersurface of the constant time t in the global coordinates (see Fig. 3). This remark is of 
a great importance in the matching procedure of the local and global coordinates discussed 
in next section. It will be also important in section 11 for correct derivation of equations of 
motion of the body which requires comparison of the definition of the multipole moments of 
the body in the local coordinates used in this paper and that in the global coordinates used 
by some other authors . 



We discuss the post-Newtonian definition of body's spin in section 10. 
See section 11 for more detail. 



7 Parametrized Post-Newtonian Transformation from Local to Global Coordi- 
nates 

1.1 Preliminary Remarks 

General relativistic post-Newtonian coordinate transformation from local, w", to global, x", 
coordinates are used in standard algorithms of data processing of various astronomical obser- 
vations to reduce the observable quantities to the barycentric coordinates of the solar system 
[21,82] where they are "stored" (catalogued) at a certain astrometric epoch recommended 
for the international usage by the lAU. The most commonly used epochs are J2000 and 
B1950 (i.e. instantaneous orientation of the barycentric coordinates fixed by the position of 
the vernal equinox as it was located on the sky in 2000 AD or 1950AD). Post-Newtonian 
coordinate transformation from the local to global frame is an inalienable part in the pro- 
cedure of derivation of the general relativistic equations of motion of test and/or extended 
bodies in the solar system [1,41,42,43,45,56,57,58,59,60,61,69,70,105,106] as well as com- 
pact relativistic stars in binary systems emitting gravitational waves [67,79,95,96,107,108]. 
It is widely accepted that scalar fields must exist since they provide elegant explanation for 
various effects and phenomena encountered by modern theoretical physics of fundamental 
particle interactions, gravitation and cosmology [109]. Therefore, it is natural to generalize 
the existing form of the general relativistic post-Newtonian transformation from the local to 
global coordinates to make them fully compatible with the scalar-tensor theory of gravity. 
Solution of this problem in terms of the PPN parameters (3 and 7 is displayed in the present 
section. 

The PPN coordinate transformations from the local (geocentric) coordinates to global (barycen- 
tric) coordinates can be found by making use of the asymptotic expansion matching technique 
[114] as proposed in [95,96]. The metric tensor and the scalar field are solutions of the field 
equations which are expressed in the form of two different post-Newtonian expansions in 
the global and local coordinates. These two expansions of the metric tensor and the scalar 
field must match smoothly in the spatial domain where both coordinate charts overlap. This 
matching domain is originally defined as a region in which the post-Newtonian expansion of 
the metric tensor and scalar field is not divergent. In the case of a weak gravitational field 
this domain extends from the origin of the local coordinates associated with a central body 
(Earth) up to the next closest celestial body (Moon) . However, after the structure of the 
PPN coordinate transformation is found the domain of applicability of the local coordinates 
can be stretched out up to much larger distances (see text below after Eq. (7.3.12) and 
papers [47,61,105] for more detail). 

Special relativistic transformation from local to global coordinates is linear and takes into 
account only kinematic aspects of the transformation, i.e. that the local frame moves with 
respect to the global one with constant velocity. The special relativistic (Poincare-Lorentz) 
transformation is described at each instant of time by 10 parameters characterizing the 
intrinsic properties of fiat space-time which is tangent to curved space-time manifold at any 



In the N-body system comprised of neutron stars and/or black holes the local coordinates would 
overlap with the global ones in the so-called buffer region in which the gravitational field of both 
the neutron star /black hole and the external bodies is weak. More details on this subject can be 
found in papers [94]-[96]. 



point on the world line of the origin of the local coordinates [48,66,76]. These parameters 
include 4 space and time translations, 3 spatial rotations, and 3 Lorentz boosts depending on 
the velocity of the origin of the local frame with respect to the global coordinates. General 
relativity generalizes the Poincare-Lorentz transformation by accounting for the presence of 
dynamic effects (acceleration) in the motion of the local frame as well as for the effects of the 
external gravitational field of the other bodies and the background curvature of the space- 
time. It was shown (see, e.g., [43,45,95,96] and references therein) that in general relativity 
the post-Newtonian coordinate transformation from local to global coordinates is non-linear, 
has more than 10 parameters, and all parameters depend on time. Scalar-tensor theory of 
gravity brings about additional complications to the relativistic theory of reference frames 
caused by the dependence of the post-Newtonian transformation on the scalar field which is 
described by two more parameters, 7 and j3. However, scalar-tensor theory does not change 
general structure of the parametrized post-Newtonian (PPN) transformation which remains 
the same as that used in general theory of relativity. 

The PPN coordinate transformation between the global and local coordinates belongs to 
the class of transformations which must comply with the gauge condition (3.4.1). Therefore, 
one begins with finding of the most general structure of such coordinate transformations. As 
soon as this structure is known it can be further specialized by reducing the number of the 
gauge degrees of freedom by making use of the residual gauge transformations and match- 
ing technique applied for establishing a one-to-one correspondence between two asymptotic 
expansions of the metric tensor and scalar field written down in the global and local coordi- 
nates. This procedure allows us to arrive to a single, ingenious form of the PPN coordinate 
transformation. We shall show that the PPN transformation from local to global coordinates 
can be represented as a power series expansion with respect to two small parameters, e = 1/c 
and r //?, where r and R are distances from the origin of the local coordinates (the central 
body) to matching point in space and to the other (external) gravitating body respectively. 
Coefficients of the power series expansion are symmetric and trace-free (STF) multi-index 
functions of time which are determined simultaneously with the external multipole moments, 
-Pl, Ql^ Cl, etc., of the local metric in the course of the matching procedure. The STF co- 
efficients of the PPN coordinate transformation are functions of the local coordinate time u 
and are "pinned down" to the origin of the local coordinates. The matching procedure shows 
that the STF coefficients of the most general form of the PPN coordinate transformation 
couples linearly with the external STF multipole moments of the local metric tensor entering 
Eqs. (6.2.19) - (6.2.22) so that the matching allows to derive a set of equations defining only 
their algebraic sum. This reflects existence of the residual gauge freedom which can be used 
in order to simplify the structure of the PPN coordinate transformation and/or that of the 
metric tensor in local coordinates. We have used this gauge freedom already in section 6.2.3 
to suppress the number of the external multipole moments which have no physical signifi- 
cance. Elimination of the non-physical multipole moments from the metric tensor leads to 
more simple structure of the PPN transformations as well. This issue is discussed in the next 
sections in more detail. 



7.2 General Structure of the Coordinate Transformation 

The most general structure of the coordinate transformation from the local, = {cu,w^), 
to global, = (ct, x*), coordinates in the weak-field and slow-motion post-Newtonian ap- 



proximation is given by two equations: 



u = t + e^^\t,x) , (7.2.1) 
w' = K + e^C{t,x) , (7.2.2) 



where and ^* are the post-Newtonian corrections to the Gahlean transformation, i?^ 
X* — Xg(t), and x%{t) is the position of the origin of the local frame at time t with respect 
to the origin of the global coordinates. We shall prove later on that the origin of the local 
coordinates can be always chosen at any instant of time at the center of mass of the body 
(Earth, Sun, planet) around which the local coordinate chart has been constructed. In what 
follows, we shall denote velocity and acceleration of the origin of the local coordinates as 
v'g = x'g and = respectively, where here and everywhere else the overdot must be 
understood as a time derivative with respect to time t. 

Pointwise matching equations for the scalar field, the metric tensor, and the Christoffel 
symbols are given by their general laws of a coordinate transformation [76,77,78] 



(p(t,x) =ip{u,w) , (7.2.3) 
9^.u{t,x)=g^p{u,w)^^ , (7.2.4) 



, ^ , ^ dx^ dwP dw"' dx^ d'^w" 
rut, X) = C(n, w)—-—-—^ + . (7.2.5) 
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One recalls that goi component of the metric tensor does not contain terms of order 0(e) 
because we have assumed that both the global and the local frames are not rotating. This 
fact, being used in Eq. (7.2.4), implies that function ^°(t, x) from Eq. (7.2.1) must be subject 
to the following restriction: = — + O(e^). This is a partial differential equation which 
can be integrated so that function can be represented as 

i\t, x) = -A{t) - vlRl + e^K{t, x) + O(e^) , (7.2.6) 

where A{t) and n{t, x) are analytic, but in all other aspects yet unspecified, functions . 

Let us now use the gauge conditions (3.4.3) in order to impose further restrictions of the 
PPN functions and ^* from Eqs. (7.2.1) and (7.2.2). The gauge conditions can be written 
in arbitrary coordinates as an exact equality 



a'-'Kp = • (7.2.7) 



The law of transformation of the Christoffel symbols, Eq. (7.2.5), being substituted to Eq. 
(7.2.7) yields a partial differential equation of the second order 



Notice that function A{t) depends only on time t. 



which describes any possible freedom of the PPN transformation in the quasi-harmonic gauge 
of the scalar-tensor theory of gravity. Let us now substitute functions = cu and from 
Eqs. (7.2.1) and (7.2.2), and ^° from Eq. (7.2.6) to Eq. (7.2.8). One obtains 



W'^Kit, x) = 'iv^al -A- dlRl + 0{e^) , (7.2.9) 

V'C{t,x) = -a', + 0{e^) . (7.2.10) 

General solution of these elliptic-type equations can be written in the form of the power series 
expansion in terms of the scalar and vector spherical harmonics (see appendix A for more 
detail). Furthermore, solution for functions n{t,x) and ^^{t,x) in Eqs. (7.2.9) and (7.2.10) 
consist of two parts - general solution of the homogeneous Laplace equation and a particular 
solution of the inhomogeneous elliptic equation. We shall omit that part of the solution of the 
homogeneous equation which is singular at the origin of the local coordinates, that is at the 
point = 0. In general, such ill-behaved terms can be present in other alternative theories 
of gravity which violate the law of conservation of linear momentum. But in the scalar-tensor 
theory of gravity such ill-behaved functions do not appear in the geocentric metric and we 
do not need singular terms in the coordinate transformations (7.2.1) and (7.2.2) to match 
the scalar field and the metric tensor in two coordinate systems. Integration of equations 
(7.2.9) and (7.2.10) resuhs in 



^= - Ia)rI - ^a'.RlRl + E{t,x) , (7.2.11) 

C = -la],Rl + E%t,x), (7.2.12) 



where functions S and S* are solutions of the homogeneous Laplace equation. These solutions 
can be written in the form of harmonic polynomials 

Eit,x) = J2-B^L>R'^, (7.2.13) 
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where the coefficients B<^l>, ^<l>, ^<l>, and <£^<l> of the polynomials are different multi- 
index STF functions. These functions are defined on the world line of the origin of the 
local coordinates and depend only on the global time t. Explicit form of these functions 
will be obtained in the process of matching of the global and local metric tensors as well 
as the scalar field. Notice that we have denoted all, yet unknown STF functions in Eqs. 
(7.2.13) and (7.2.14) by capital calligraphic letters, while the STF multipole moments in 
the expressions for the local metric, Eqs. (6.2.36) - (6.2.38), and the scalar field (6.2.4) 
have been denoted by capital roman letters. This is supposed to help to distinguish the STF 
functions having different origin. We also emphasize that the specific form of functions in the 
coordinate transformations (7.2.11) and (7.2.12) is the most general one which preserves the 
gauge conditions (7.2.7) and allows us to analyze the residual gauge freedom in construction 
of the local coordinates by operating with the same type of functions which have been used 
in the external solution for the metric tensor in the local coordinates. 



1.3 Transformation oj the Coordinate Bases 



Derivation of the PPN coordinate transformation originates from the matching equation 
(7.2.4) for the metric tensor apphed in the joint domain of vahdity of the local and global 
coordinates. This equation contains the matrix of transformation A^^ = dw^/dx°' between 
the two coordinate bases, = d/dw°' and = d/dx°', in the local, w"", and global, x", 
coordinates respectively. Transformation between them reads 



In that region of space-time where the matrix of transformation A^^, is non-singular it can 
be inverted so that one gets the inverse matrix (A~^)^ defined by a standard rule [76] 



The inverse matrix is required to get the inverse transformation between the two bases. The 
matrix A"^ can be expanded in the post-Newtonian series, which is a consequence of the 
post-Newtonian expansion of the coordinate transformation described in a previous section. 
The post-Newtonian expansion of the matrix of transformation is as follows 
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where the coefficients of the expansion are the following functions of the global coordinates: 
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Elements of the inverse matrix (A"^)*^ can be deduced from Eqs. (7.3.7) - (7.3.11) by ap- 
plying relationship (7.3.2). 

Formulas (7.3.7) - (7.3.11) allow us to evaluate the range of applicability of the local coordi- 
nates. Radius of this range is determined by the condition that determinant of the coordinate 
transformation matrix A*^ is zero. Calculating the determinant of the matrix one obtains 



det (A",) 



-A + 3S 



(^B ~ ) Rl 



(7.3.12) 



- ([+l)(2[+3) 



1=2 



(2/ + 1)/! 



S^T.->R^ 



O e 



Radius of convergence of the polynomial in the right side of Eq. (7.3.12) crucially depends 
on the choice of functions £l. In what follows , we shall prove that it is possible to make 
function £i = and all other functions Sl = for any / > 2. Then if one also takes into 
account (8.5.2) and (8.5.8) (putting Q = Y = there) the determinant (7.3.12) turns to 



See Eqs. (8.5.9), (8.5.10) 



zero when the distance case of the local geocentric frame attached to 

the Earth and moving around the Sun with acceleration ^ 0.6cm/s^ this distance Rg is 
about 10^^ cm. Hence, the local geocentric frame covers a region which includes the entire 
solar system. In case of a binary pulsar with a characteristic size of the orbit ~ 10^'^ cm 
the local coordinate system attached to the pulsar is spread out of the binary system at 
the distance about 10^^ cm which also significantly exceeds the distance between the pulsar 
and its companion. This remark can be important for researchers doing analysis of physical 
processes going on in pulsar's magnetosphere [117]. 

This consideration suggests that the metric tensor defined originally in the local coordinates 
only in the domain restricted by the distance to the nearest external gravitating body can be 
extrapolated far beyond this boundary. Such extrapolation can be accomplished by choos- 
ing another form of the solution of the homogeneous field equations describing background 
gravitational field of external bodies. Research in this direction has been pursued in [105] . 



8 Matching the Post-Newtonian Expansions of the Metric Tensor and a Scalar 
Field 

8.1 Historical Background 

The method of matched asymptotic expansions has been developed for finding solutions of 
ordinary and/or partial differential equations and is well known in mathematics for a long 
time (see, for instance, [114], [115] and references therein). The idea of implying this method 
in general relativity goes back to earlier work of Einstein and Rosen [118] where the authors 
discussed the problem of motion of gravitationally interacting particles by treating them as 
topological structures ('bridges') in the space-time manifold being a regular solution of the 
Einstein field equations. Fock [48] applied the matching technique to join the metric tensor 
expansions of the 'near' and 'far-radiative' zones of an isolated astronomical system emitting 
quadrupole gravitational radiation (see also [84] where this procedure is discussed in more 
detail). Manasse [119] studied radial fall of a small black hole onto a massive gravitating 
body and calculated distortion of the shape of the black hole's horizon by making use of the 
matching technique. Thorne [120] and Burke [121] suggested to use the matching technique 
for imposing an outgoing-wave radiation condition on the post-Newtonian metric tensor for 
an isolated system emitting gravitational waves. This method helps to chose a causal solu- 
tion of the homogeneous Einstein equations in the post-Newtonian approximation scheme 
and to postpone appearance of ill-defined (divergent) integrals, at least, up to the fourth 
PNA [122,123,124]. Demiansky & Grishchuk [125] used the matching technique to show that 
a black hole orbits its companion of a comparable mass in accordance with the Newtonian 
equations of motion. At about the same time, D'Eath [95,96] explored the idea proposed 
by Hawking and worked out a detailed analysis of the problem of motion of two Kerr black 
holes comprising a binary system by making use of matching of the internal (local coordi- 
nates) and external (global coordinates) solutions of the Einstein equations. D'Eath derived 
general relativistic equations of motion of the black holes in the first post-Newtonian (1 
PN) approximation. Kates [110,111] extended his analysis and obtained the gravitational 
radiation-reaction force (2.5 PNA) for the black holes in a binary system. He has also elabo- 
rated on a rigorous mathematical treatment of the matched asymptotic expansions technique 



for various applications in general relativity [112]. Damour [79] used the asymptotic match- 
ing to solve the problem of motion of two spherically-symmetric and non-rotating compact 
bodies with the gravitational radiation reaction force taken into account. He proved that 
mass of each of the bodies, which appears in the external solution of the two-body problem 
as a constant parameter, is the same as that characterizing the Schwarzschild metric of a 
non-rotating black hole. Thorne and Hartle [94] applied the matching technique to study 
the problem of translational motion and precession of compact bodies having quadrupole 
gravitational fields. Their method combined with the mathematical technique of D'Eath [95] 
was employed in [45] to derive the post-Newtonian equations of motion of extended bodies 
making up an N-body system in the weak-field and slow-motion approximation. The paper 
[45] also demonstrates for the case of N-body problem how to construct a local coordinate 
system with the origin moving exactly along the world line of the center of mass of an ex- 
tended body which has arbitrary shape and is rotating. The matching technique used in 
[45] led to the development of the Brumberg-Kopeikin (BK) formalism in the theory of as- 
tronomical reference frames and have been improved in a series of subsequent publications 
[58,59,60,61,101,105]. Similar matching technique and Damour-Soffel-Xu (DSX) formalism 
were formulated in [69,70] to describe the post-Newtonian celestial mechanics of isolated 
astronomical systems. Both BK and DSX formalisms were used later as a basis of the res- 
olutions on the relativistic reference frames and time scales in the solar system adopted 
by the 20-th General Assembly of the lAU [14]. In the present paper we extend the gen- 
eral relativistic theory of reference frames and apply the matched asymptotic expansions 
technique to the case of the scalar-tensor theory of gravity. This will allow us to incorpo- 
rate unambiguously the adopted lAU resolutions on reference frames in the solar system to 
the parametrized post-Newtonian (PPN) formalism [20], thus, making more close the link 
between the experimental gravity and modern observational astronomy. 

For the benefit of readers we emphasize here the differences between the results obtained in 
BK formalism [45,58,59,61] and in this paper and those contained in DSX formalism [69,70] 
and in [1,65]. These differences follows: 

• Original BK formalism deals with the Newtonian definitions of the multipole moments 
while DSX formalism operates with the post-Newtonian multipoles. This paper extends 
the BK formalism and incorporates the post-Newtonian multipoles into the matching 
technique and equations of motion. 

• BK and DSX formalisms have been developed in the framework of the general theory of 
relativity only. The present paper extends the BK formalism to the class of scalar-tensor 
theories of gravity. 

• Damour and Esposito-Farese [65] generalized definitions of the post-Newtonian multipole 
moments for an isolated N-body system for the case of the scalar-tensor theory of gravity 
and concentrated on discussion of experimental tests of this theory for binary pulsars 
and gravitational wave astronomy. They did not work out any matching procedure for 
construction of local frames of reference and derivation of equations of motion of each 
body in the local coordinates. The present paper develops the matching procedure in the 
scalar-tensor theory of gravity and constructs a set of global and local coordinate frames 
for description of both global and local dynamics of the N-body system as a whole and 
each body from the system separately. In the present paper we also construct the post- 
Newtonian definition of the multipole moments for each body being a member of the 
N-body system. 

• Klioner and Soffel [1] used DSX formalism to supplement two-parametric PPN formalism 



by the DSX matching technique and construct a set of global and local coordinate frames 
for description of the dynamics of the N-body system. They also used PPN-parametrized 
definitions of the post-Newtonian multipoles given in [65]. Klioner and Soffel did not 
rely upon a particular class of alternative gravitational theory and abandoned the use of 
gravitational field equations. The present paper makes use of the field equations of the 
scalar-tensor theory of gravity and applies the matching technique of the BK formalism. 
We show that Klioner- Soffel results [1] do not match with the scalar-tensor theory of 
gravity which makes domain of applicability of the phenomenological theory of reference 
frames constructed in [1] uncertain. Further discussion of the Klioner- Soffel approach and 
comparison of results of their paper [1] with those obtained in the present review article 
are deferred to appendix C. 

8.2 Method of the Matched Asymptotic Expansions in the PPN Formalism 

Method of the matched asymptotic expansions is one of the powerful mathematical tools 
for solving differential equations with a small parameter present for which a regular pertur- 
bation series method is inadequate. It occurs as often as a solitary solution can not match all 
the boundary conditions in a differential equation [114], [115]. If a regular perturbation series 
can not be applied, there may be one or more regions where the solution can be represented 
in the form of asymptotic expansion with respect to one or more small parameters which 
satisfies at least one of the boundary conditions. Matching the asymptotic expansions in the 
buffer region (boundary layer) where at least two of them are valid and convergent, allows us 
to find the law of transformation from one expansion to another and to retrieve coefficients 
of these expansions. 

In the present paper the main asymptotic post-Newtonian expansions, which are used in the 
matching procedure, are solutions of the gravity field equations for the metric tensor and 
scalar field found in the global and local coordinates that are subject to different bound- 
ary conditions imposed respectively at infinite distance and/or at the origin of the local 
coordinates. These solutions are shown in Eqs. (4.2.1) - (4.2.3) and (6.2.3), (6.2.4), (6.2.36) 
- (6.2.38). The solution for the metric tensor and scalar field in the global coordinates is 
valid everywhere inside and outside of the N-body system up to infinity. This is because 
we assumed that the gravity field is weak everywhere and there are no singularities in the 
space-time manifold . 

Because we do not deal with space-time singularities the reader may think that one could 
use the global coordinates alone to describe relativistic celestial dynamics of the bodies from 
the N-body system. However, this idea does not work out for two reasons. First, the local co- 
ordinates are still required to give physically meaningful definition of the multipole moments 
of each body and one must know how these definitions correlate with the definitions of these 
multipoles in the global coordinates. This relationship between the two definitions of the 
multipole moments is a key element in the procedure of derivation of equations of motion of 

It is also called the boundary layer method [113]. 

Had one worked with the tensor of energy-momentum of the point-like massive particles one would 
have singularities on the particle's world lines [116]. We argue that in such a case the method of 
the matched asymptotic expansions is the only way to derive equations of motion of these particles 
without ambiguities present in the higher-order post-Newtonian approximations [126]. 



extended bodies having finite size. We discuss importance of tliis issue in section 11 in more 
detail. Second, the global coordinate frame is an inappropriate reference for analysis of gravi- 
tational experiments conducted in the vicinity of the Earth. This is because the Earth is both 
moving and embedded to the gravitational field of other bodies of the solar system. Simple 
translation of the origin of the global coordinates to the geocenter (the Galilean-Newtonian 
transformation) frequently used in early publications does not take into account relativis- 
tic aspects of the coordinate transformations on the space-time manifold and, hence, can not 
eliminate a large number of coordinate-dependent (and for this reason nonphysical) effects 
which will be present in such treatment of astronomical observations [45,59,106,130]. The 
number of the coordinate-dependent effects is much smaller if one uses correct relativistic 
procedure to transform the global to local coordinates. Such post-Newtonian transformation 
simplifies drastically analysis of astronomical observations and description of the dynamics 
of lunar motion and/or that of artificial satellites [61,131,132]. 

The internal solution for the metric tensor and scalar field in the local coordinates contain 
the external multipole moments which can not be found as explicit functions of time without 
matching of the local solution of the gravity field equations to the global one. The matching 
allows us to express the external multipole moments in terms of the gravitational potentials 
(4.2.7) characterizing the global metric tensor and scalar field. At the same time the matching 
procedure determines the structure of the PPN coordinate transformation between the global 
and local coordinates. 

Matching of the local and global solutions of the metric tensor and scalar field is based on 
Eqs. (7.2.3) and (7.2.4), and consists of the following steps : 

Step 1. One re-writes the local metric tensor and scalar field in the right side of Eqs. (7.2.3) 
and (7.2.4) in terms of the global coordinates {t, x). This is achieved by making use of 
a Taylor expansion of (p{u, w) and gai3{u, w) around the point = (ct, x). 

Step 2. One calculates the partial derivatives of the local coordinates with respect to the global 
ones, that is the matrix of transformation of the coordinate bases given in section (7.3). 

Step 3. One separates the 'global' gravitational potentials in the left side of Eqs. (7.2.3) and 
(7.2.4) relating to the central body (Earth) and to the external bodies (Moon, Sun, 
etc.) respectively: 

U{t,x) = m^\t,x) + U{t,x) , (8.2.1) 
U'{t, x) = Ul^\t, x) + U\t, x) , (8.2.2) 

X{t,x) = x'^\t,x)+x{t,x), (8.2.3) 
<l>fc(t, x) = $r(t, x) + <lfc(t, x), (A; = 1, ... , 4) , (8.2.4) 



See, for example, [20,48,127,128,129] and references therein. One may notice, however, that Will 
understood this problem fairly clear long time ago as follows from the discussion in section 6.2 of 
his book [20]. 

Solution of some problems (for instance, in cosmology [133]) requires to match not only the metric 
tensor but also its first derivatives. In this paper this requirement is redundant as the metric tensor 
and scalar field are smooth differentiable functions in the matching domain and their derivatives 
of any order have no jumps. Therefore, Eq. (7.2.5) is a consequence of Eq. (7.2.4). 



where functions with index (B) are given by integrals (4.2.8)-(4.2.15) taken over the 
volume of the central body only, and the bar over other functions indicates here and 
hereafter that the corresponding sum in the definitions (4.2.7) of these functions excludes 
the central body with the index (B), that is the sum takes into account only external 
bodies 

U=Y: U^^^ , f/. = E Ur , ^k=T.^i"'^ X = E . (8.2.5) 

Aj^B Aj^B A^B A^B 

Step 4. One expands the gravitational potentials of the external masses (that is functions with 
bars in Eqs. (8.2.1) - (8.2.5)) in Taylor's series in powers of i?^ = — in the vicinity 
of the origin of the local coordinates, that is the point = x^. 

Step 5. One equates similar terms of these Taylor expansions from the left side of the matching 
equations (7.2.3) and (7.2.4) with the corresponding Taylor expansions entering their 
right side. 

Step 6. One separates symmetric and anti-symmetric tensor parts in the matching equations 
and determine all coefficients in the local metric tensor and scalar field, which remained 
undetermined so far, as well as coefficients in the coordinate transformations. This fixes 
the residual gauge freedom and gives equations of translational and rotational motion 
of the local reference frame. 

Let us now explain each step of the matching procedure in more detail. 



8.3 Transformation of Gravitational Potentials from the Local to Global Coordinates 



8. 3. 1 Transformation of the Internal Potentials 

At the first step of the matching procedure one has to transform the metric tensor and the 
scalar field in the right side of matching equations (7.2.3) and (7.2.4) from the local, = 
(cu, lu), to global, x" = (ct,a;), coordinates . We remind that the internal gravitational 
potentials associated with the scalar field, Eq. (6.2.3), and the metric tensor, Eqs. (6.2.9) - 
(6.2.16), are defined in the local coordinates w° = (cm, w) as integrals over hypersurface of a 
constant coordinate time u. On the other hand, the corresponding gravitational potentials, 
Eqs. (4.2.8) - (4.2.14) are defined in the global coordinates x" = (t, a;) as integrals over 
hypersurface of a constant coordinate time t. These two hypersurfaces do not coincide and 
can intersect only at the points that form a 2-dimensional sub-hypersurface (see Fig. 3). For 
this reason, in order to transform the internal potentials defined in the local coordinates, 
= {cu,w), to those defined in the global coordinates, x° = {ct,x), one needs to make 
a pointwise transformation given by Eqs. (7.2.1) - (7.2.2) along with a Lie transform of 
integrands of the integrals which moves the integrands from the hypersurface of constant 
time u to that of constant time t (see Fig 4). This procedure was worked out in [58] and 
is described below. 

It is also conceivable to make a reciprocal transform of all functions in the left side of equations 
(7.2.3) and (7.2.4) to the local coordinates w" = {cu, w). However, it is more convenient and simpler 
to transform the metric tensor and scalar field from the right side of Eqs. (7.2.3) and (7.2.4) to 
the global coordinates = {ct,x) in accordance to the transformations (7.2.1), (7.2.2), (7.2.11) - 
(7.2.14) which are already displayed in terms of the global coordinates. 
30 See also [61,134]. 



Let us assume that the field point P, at which matching of the internal and external solutions 
of the metric tensor and the scalar field is done, has global coordinates x"(P) = {ct, x) and 
local coordinates w°'{P) = {cu, w) (see Fig. 4). These coordinates are related by the pointwise 
transformation given by Eqs. (7.2.1) - (7.2.2). By definition, the matching point P belongs 
simultaneously both to the hypersurface of constant time u and to that of the constant 
time t. Let us consider a matter element of the central body B located at the point N and 
lying on the same hypersurface of the constant time u inside of the body's world tube. We 
assume that the point N has local coordinates, w"(N) = (cm, w'{u)), and global coordinates, 
a;"(N) = [ct', x'{t')). We emphasize that the time coordinate u of the points P and N is the 
same because they are located on the same hypersurface of the constant time u. However, 
the time coordinate t' of the point is different from the time coordinate t of the point 
P, that is f 7^ t, because the hypersurfaces of the constant global time passing through the 
points P and N are different. Let us consider a world line of the element of the body's matter 
passing through the point N and intersecting a hypersurface of the constant time t at the 
point M. This world line allows us to map coordinates of the element of the body's matter 
from the hypersurface t' to the hypersurface t. By the construction, the point M must have 
the global coordinates x"(M) = {ct,x'{t)). 

One can expand the spatial coordinates of the world line of the chosen element of the body's 
matter in a Taylor series with respect to time 

x'\t') = x'\t) + v'\t){t' -t)+0 [At^) , (8.3.1) 

where f is the spatial velocity of the body's element at the point M and At = t' — t. One 
denotes by letters and Q position of the origin of the local coordinates on two hypersurfaces, 
t and t', respectively (see Fig. 4). Global coordinates of the origin of the local coordinates, 
taken on two different hypersurfaces, are related by equation 

x'^(t') = xi{t) + vliW -t) + [At^) . (8.3.2) 



Now we shall find the time interval At = t' — t separating the two hypersurfaces of the 
constant global coordinate time, t and t' , under condition that the matching point P is fixed. 
At the point N the relationship between the local time u and the global time t', according 
to Eq. (7.2.1), is 

u = t' + e'^\t',x') + 0(e^) . (8.3.3) 



Subtracting Eq. (7.2.1) from Eq. (8.3.3) and accounting for the fact that the space-time 
interval between points N and P is small we obtain 



At = t' -t = e^ 



(8.3.4) 



^<^(t,x)-eit,x')\+0(e^' 

where the point x' = x'{t), that is, it belongs to the hypersurface t. 

Local coordinates of the point N are transformed to the global coordinates as follows 

w" = x'\t') - Xs{t') + e^C {f, x'{t')) . (8.3.5) 

Expanding functions in the right side of Eq. (8.3.5) in the vicinity of the time instant, t, and 
taking into account Eqs. (8.3.1) - (8.3.4) yields 



w 



(8.3.6) 



where = ^(t, x'), ^'^ = x'), and R'^ = x''{t) - xi{t)- 



Transformation (8.3.6) is used for deriving relationship between the absolute values of func- 
tions \w' — w\ and | a;' — a3 1 which enter denominators of all integrands in the integrals defining 
the internal gravitational potentials. Subtracting Eq. (7.2.2) from Eq. (8.3.6) and taking the 
absolute value of the difference gives 



\w — w\ 



\x — x\ + e 



O e 



(8.3.7) 



where = ^*(t, x'), ^* = ^*(t, x), and n* = (x'* — x'^)/\x' — x\. 



We must also perform a Lie transform to find a relationship between the volume elements 
d^w' and dP'x' taken at the points N and M respectively. We note that the invariant density 
p* introduced in Eq. (3.3.18) has one more remarkable property in addition to Eq. (3.3.19) 
[58]. Specifically, the Lie derivative of the product of the invariant density and the volume 
element is zero, which means that 



p* (t, x'it)) c/V(M) = p* it', x'it')) rfV(N) = p* (m, w') (fw'iU) . 



(8.3.81 



Locally measurable velocity z/'* of the body's element at the point N is defined with respect 
to the origin of the local coordinate system. It relates to velocity f of the same element of 
the body taken at point M in global coordinates by equation: 

y'\u) = v'\t) - v'^{t) + O (e^) , (8.3.9) 
which can be derived from Eq. (8.3.6) by direct differentiation with respect to time. 

By making use of Eqs. (8.3.7) - (8.3.9) we finally obtain transformations of the internal 
gravitational potentials from the local to global coordinates. These transformations are given 
by the following equations 



f/(S), 



M, W 



ur{u,w 



(n, w 



^Ij"" (m, w 



U^''\t,x) + e^U'^^\t,x) + 0{e^) , 
Ul^\t,x)-vlU'-\t,x)+0{^), 
x'-\t,x) + 0{e'), 

$r (t, + vlU'^\t, x) - 2vlUi,^{t, X) + 0{e^) 
¥,'^\t,x)-GV^l{p*U{t,x)} + 0{e^) , 
:$r(t,a;) + 0(a. 



¥i'\u,w) = ^f\t,x) + 0{e^) , 
where the post-Newtonian correction, U^'^^t, x), to the Newtonian potential, U^'^^lt, x 

U'^^ {t, x) = Gl'^^ {p*n'' [v"" - vl) - e) - P*n'' " } • 



8.3.10) 
8.3.11) 
8.3.12) 
8.3.13) 
8.3.14) 
8.3.15) 

8.3.16) 

reads 

8.3.17) 



This correction is the result of the post-Newtonian coordinate transformation (7.2.1), (7.2.2) 
and the Lie-transfer of the integrand of the Newtonian gravitational potential from local to 
global coordinates. Transformation of all other internal gravitational potentials from local 
to global coordinates does not require to take into account their relativistic corrections as it 
would exceed the accuracy of the first post-Newtonian approximation which is redundant. 

The matching procedure also requires to derive transformation of the second time derivative 
of the potential x iii explicit form. This transformation can be directly obtained from the 
definition of the potential given by Eq. (4.2.10) and the mapping equation (8.3.12). After 
straightforward calculation one gets 

^^'2^' = c^x!oo'(i, ^) + alxfit, x) + 2cvlxfuit. ^) + ^>ix3'(t, x) + 0{e') , (8.3.18) 

where x\o^ ^ind x^T denote partial derivatives of x'^' with respect to the global time coordi- 
nate = ct and the spatial coordinate respectively. 



8.3.2 Transformation of the External Potentials 

External potentials in the internal solution for the metric tensor in local coordinates u;" = 
(cti, It)) depend on the external multipole moments Ql = Ql{u), Cl = Cl{u), Pl = Pl{u) 
depend only on time u and are 'pinned down' to the origin of the local coordinates located 
at the point B where the hypersurface of the constant time u intersects with the world line 
of the origin (see Fig 4). On the other hand, all functions entering the left side of matching 
equations (7.2.3), (7.2.4) are defined on the hypersurface of constant global coordinate time 
t. Hence, before doing the pointwise coordinate transformation of the external potentials 
we must perform a Lie transfer of the potentials along the world line of the origin of the 
local coordinates from point B to point A located on the hypersurface of the constant global 
coordinate time t. Time shift At along this world line is determined by Eq. (8.3.4) where 
one has to associate the point x' with the origin of the local coordinates, that is x' = x^, 
under condition that the matching point P is taken the same as defined in previous section. 
Keeping in mind that the external potentials are scalars with respect to the Lie transport, 
one obtains 

Ql(B) = QLiA) + gL(A)At + O [At'') , (8.3.19) 

where overdot means differentiation with respect to time t. After making use of Eqs. (7.2.1) 
and (8.3.4), and accounting that local coordinates of the point B are, w°(B) = {cu, 0), while 
the global coordinates of the point A are, x"(B) = {ct,Xg(t)), one gets 

Ql{u) = QL{t) - e'QLit) [e{t, xg) - e°(t, x)] + O (e^) . (8.3.20) 

Similar formulas are valid for the multipole moments Cl and Pl as well. 

Now we can do the pointwise coordinate transformation of the space coordinates given by 
Eq. (7.2.2). For the STF product of the local coordinates one has: 

^<hi2-k> ^ j^<hi2...H> ^ le^Rf^'^-'^-^^n> + O (e^) , (8.3.21) 



where we have used Eq. (7.2.2). After combining Eqs. (8.3.20) - (8.3.21) all together the post- 
Newtonian transformation of the Newtonian part of the external potential of the internal 
solution of the metric tensor assumes the following form 



oo -| oo 
1=0 1=0 



oo -| 



1=0 



(8.3.22) 



+e^y 

This is the most comphcated transformation of the external potential that we need. It takes 
into account the post-Newtonian nature of the PPN coordinate transformation and sup- 
plements Eq. (8.3.10) for the internal Newtonian gravitational potential. All other external 
potentials present in the local metric (6.2.36) - (6.2.38) are transformed without taking into 
account the post-Newtonian corrections by making use of only the very first term in the 
right side of Eq. (8.3.21). 



8.4 Matching of a Scalar Field 



Scalar field is used in the post-Newtonian terms only. For this reason matching of its asymp- 
totic expansions derived in local and global coordinates is quite straightforward. We operate 
with external and internal solutions of the scalar field given by Eqs. (4.2.1) and (6.2.3), 
(6.2.4) respectively. Matching equation (7.2.3) reveals that the internal potentials referred 
to the central body B cancel out in its left and right sides due to Eq. (8.3.10) while the 
potentials Pl match to the Newtonian potential and its derivatives referred to the external 
bodies. More specifically, for any number / > the matching yields 

PL = UA^B) + Oie^), (8.4.1) 

where the external Newtonian potential tj is defined in Eq. (8.2.5) and taken at the origin 
of the local coordinates, that is the point = x%(t), at the instant of time t. Thus, the STF 
multipole moment Pl of the external scalar field is fully determined as l-th spatial derivative 
of the Newtonian gravitational potential U. We remind that the scalar field was normalized 
to the factor 7 — 1 , so that physically observed scalar field C = (1 — 7)^9. It vanishes in 
the limit of general relativity where the multipole moments Pl play no role. 



8. 5 Matching of the Metric Tensor 



8.5.1 Matching ofgQQ{t,x) and gai3{u,w) in the Newtonian approximation 

Newtonian approximation of the matching equation (7.2.4) for goo{t,x) component of the 
metric tensor in its left side yields 

N{u, w) = N{t, x) + 2«B(t, x)-vl + 0{e^) . (8.5.1) 

Function ^(t, x) is taken from Eq. (7.3.7) while the global and local metric tensors are taken 
from Eqs. (4.2.2) and (6.2.36) respectively. One finds that after making use of Eq. (8.3.10) 
the internal gravitational potentials U^^\u, w) and U^^\t, x) drop out of the left and right 
sides of Eq. (8.5.1). Expanding the external gravitational potential U(t, x) in a Taylor series 



See Eq. (3.3.10) where a; + 2 = 1/(1 — 7) in accordance with definition of the PPN parameter 7. 



around the origin of the local coordinates, x^, and equating similar terms with the same 
power of specifies the matching conditions 



Q + A=-vl + Uix^) + Oie'), (8.5.2) 
Q, = f/,i(a;s)-at + 0(e2), (8.5.3) 
QL = U,L{^B)+0{e^), (/>2). (8.5.4) 



Equation (8.5.2) makes it evident that two functions A and Q can not be determined from 
the matching procedure separately, only their linear combination Q + A can be determined. 
Hence, one of these functions can be chosen arbitrary. The most preferable choice is to 
take Q = as it was done, for example, in [45,69,94]. This choice is also consistent with 
recommendations of the lAU [14,15] and it makes goQ{u,w) component of the local metric 
tensor equal to —1 at the origin of the local coordinates if gravitational field of the central 
body (Earth) is neglected. However, if one chooses Q = the rate of the coordinate time u 
can be different from that of the coordinate time t because the average value of functions 
and U is not zero for elliptic orbits. Hence, the choice oi Q = can be inconvenient 
for astronomical data reduction programs in the solar system. Therefore, two time scales, 
TDB = kst and TDT = ksu, have been introduced in such a way that their rate at the 
origin of the local coordinate system is the same . This makes function 



Q{t) = - <A>=a + bt + ct'^ + ... , (8.5.5) 

to be a polynomial of time which numerical coefficients are calculated by means of numerical 
integration of Eq. (8.5.2) over sufficiently long interval of time [99,101]. Time-rate adjustment 
coefficients and kE relate to each other as [99,101] 

^ = l + c-2Q(r), (8.5.6) 
kE 

where t* is a certain astronomical epoch chosen by convention. 

In accordance with the interpretation given in [76,94], function Qi from (8.5.3) must be 
understood as acceleration measured by accelerometer being rigidly fixed at rest in the 
origin of the local coordinates under condition that the internal gravitational field of the 
central body B is neglected. The choice, Qi = 0, (see, e.g., [42,94]) leads to construction of 
a freely falling local coordinate frame which origin moves along a geodesic world line in the 
background space-time manifold defined by the gravitational potentials of all the celestial 
bodies of the N-body system but the central body B that is located near the origin of the 
local coordinates. Such choice of Qi, however, disengage the world line of the center of mass 
of the central body B from that of the origin of the local coordinates (see Fig. 5). Indeed, 
the center of mass of the central body (Earth) does not move along the geodesic world line 
due to the interaction of its internal quadrupole and higher-order multipole moments with 
the tidal gravitational field of external bodies [45,47,58]. For this reason, a special procedure 
must be applied for finding function Qi which will make the origin of the local coordinates 

32 See [14,15,99,101] for more detail. 



be following the same world line as the body's center of mass. This procedure is described 
later in section 9.4 in full details. 

STF external multipoles Ql are defined in the Newtonian approximation by Eq. (8.5.4). They 
describe gravitoelectric [106,135] tidal gravitational field of external bodies at the origin of the 
local coordinates. Post-Newtonian corrections to the Newtonian value of the multipoles Ql 
can be important for lunar laser ranging and other experimental tests of relativistic gravity 
theory in the near-Earth space. The post-Newtonian corrections to the external multipoles 
can be also important in construction of the relativistic theory of accretion disc around a 
star in a close binary system [136]. These corrections are derived later and shown explicitly 
in Eq. (8.5.34). 

8.5.2 Matching ofgij{t,x) and gai3{u,w) to the order of 0{e^) 

We substitute gij{t,x) component of the metric tensor in the global coordinates to the left 
side of the matching Eq. (7.2.4) and write this equation by taking into account all post- 
Newtonian terms of order O(e^). We find that in this approximation the matching equations 
reads 

Hij{u, w) = Hij{t, x) - ^)it, x) - ^K{it, x) + vivi + 0(e2) , (8.5.7) 

where the matrix of transformation 9^*^(t,a;) from the local, w"', to global, x°', coordinates 
is given in Eq. (7.3.11) whereas components of the metric tensor are defined by Eq. (4.2.5) 
in the global coordinates and by Eq. (6.2.38) in the local coordinates. Comparing similar 
terms we find that the internal gravitational potentials of the central body B, tl^'^^u^w) 
and f/'^'(t, a;), entering the left and right sides of the matching equation, cancel each other 
and drop out of the matching equation completely. The remaining terms belong to the 
background gravitational field of external bodies and their matching gives the following set 



of equations 

Y + £ = \vl + 7t/(x,) + 0{e^) , (8.5.8) 
6 

8i = a\ + 0{e^) , (8.5.9) 

£L = 0{e^) (/>2), (8.5.10) 

V,, = ]^vfvi> + 0{e'), (8.5.11) 

VL = 0{e^) (/>3), (8.5.12) 

TL = 0{e^) il>2), (8.5.13) 



which define the structure of the PPN transformation between spatial coordinates of the 
global and local coordinates. 

Function Y = Y{t) defines the unit of measurement of spatial distances in the local coordi- 
nates. It would be convenient to chose F = as it was done, for example, in [45,94]. However, 
introduction of TDB and TDT time scales in ephemeris astronomy must be accompanied 



by the choice of y 7^ to compensate effects of function Q 7^ from Eq. (8.5.2) on spatial 
barycentric (global) and geocentric (local) coordinates [101,102]. In view of this complexity, 
related to a pure practical application of the PPN theory of reference frames in the solar 
system, we do not impose any restriction on function Y at this step of the matching proce- 
dure. It will allow us to trace how functions Q and Y affect equations of motion of extended 
bodies. 

8.5.3 Matching ofgQi{t,x) and gai3{u,w) to the order of 0{e^) 
Matching equation (7.2.4) for gQi{t,x) taken in its left side, reads 



Li{u, w) = Li{t, x) + Q3i(t, x) - %{t, x) + v\ 



Hij{t,x)-m{{t,x) 



^{t,x) + N{t,x) 



(8.5.14) 



where we have employed Eqs. (8.5.1) and (8.5.7) for doing simplification. Subsequent substi- 
tution of the metric tensor given in the global, Eqs. (4.2.2), (4.2.4) and (4.2.5), and the local, 
Eqs. (6.2.7) and (6.2.37), coordinates, to Eq. (8.5.14) along with Eqs. (8.2.2) and (8.3.11) for 
vector-potentials Ui{t, x) and ul'^\u, w) respectively, reveals that all gravitational potentials 
depending on the internal structure of the central body B (Earth) cancel each other. Remain- 
ing terms in the matching equation (8.5.14) depend only on the background values of the 
gravitational potentials of external bodies. They yield a number of relationships which allow 
us to express the external multipole moments from the metric tensor in local coordinates 
and functions present in the PPN coordinate transformations (7.2.11)-(7.2.14) in terms of 
the gravitational potentials of external bodies. These relationships are as follows 



= 2(1 + 7)^'(^b) - (1 + 2^)viU{x,) - ^vivl - viQ + 0{e') 
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(8.5.15) 
(8.5.16) 
(8.5.17) 



and 



e^pk {C, + j-p) = -2(1 + 7)f/''''k^B) + (1 + 2^)vlU^^\x,) + vIQ^^ + 0{e^) 
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(8.5.18) 



(8.5.19) 



As one can see the matching equation (8.5.14) gives rise to two sets of equations, Eqs. 
(8.5.15) - (8.5.17), and Eqs. (8.5.18 - (8.5.19) which have different properties of symmetry. 



More specifically, Eqs. (8.5.15) - (8.5.17) depend only on objects which are symmetric and 
trace-free with respect to the entire set of indices {i, 01,02, •••,az}. On the other hand, Eqs. 
(8.5.18), (8.5.19) are symmetric with respect to the set of indices {ai, 02, ai} but they are 
anti-symmetric with respect to any pair of indices consisting of the index i and either one off 
the set {oi, 02, a^}. This separation of the matching equation (8.5.14) in the symmetric 
and antisymmetric components is due to the fact that this equation has one free (vector-type) 
index i and, hence, can be split in fully symmetric and antisymmetric parts in accordance 
with the Clebsch-Gordan decomposition of the vector spherical harmonics into irreducible 
representations (see [62,83,137] for more details). 

Let us now discuss a freedom in choosing angular speed of rotation of spatial axes of the 
local coordinates. From the discussion in previous sections one knows that there are two 
possible types of the local coordinates - either dynamically or kinematically non-rotating. 
Dynamically non-rotating local coordinates are defined by the condition that the external 
dipole moment Cj = in Eq. (6.2.29). On the other hand, the kinematically non-rotating 
local coordinates are realized if one chooses function jFj = in the coordinate transformation, 
Eq. (7.2.12), for spatial axes. If we do not fix the choice of Cj, then, Eq. (8.5.18) can be 
rewritten as 



(1 + i)W''{x,) - (7 + iHU'^ix,) - \vlQ' 



C, + 0{e^). (8.5.20) 



The first term in Eq. (8.5.20) describes the Lense-Thirring gravitomagnetic precession, the 
second term describes the de Sitter precession, and the third term describes the Thomas 
precession [20,57,76]. We emphasize that in the scalar-tensor theory both the Lense-Thirring 
and the de Sitter precessions depend on the PPN parameter 7 while the Thomas precession 
does not. The reason is that the Thomas precession is a special relativistic effect [76] and, as 
such can not depend at all on a particular choice of a specific gravitational theory. If function 
Ci = 0, spatial axes of the local coordinates are kinematically rotating. On the other hand, 
the choice jFj = gives kinematically non-rotating local coordinate frame which spatial axes 
are fixed with respect to distant quasars with neglected proper motions [81]. 

Functions Bl enter the PPN time coordinate transformation and depend on the gauge- 
dependent multipole moments Zl which can be chosen arbitrary. There are two of the most 
preferable options: 

1) One simplifies the time transformation, Eq. (7.2.11), from the local coordinates to the 
global ones as much as possible. In this case the moments Zl have to be chosen such 
that functions Bl in the time transformation assume the most simple form 

B^,k> = 2vTaf + 0{e') (8.5.21) 

S<^L> = 0(e2), (/>2). (8.5.22) 

Here the choice of Bij in Eq. (8.5.21) eliminates all terms explicitly depending on the 
velocity of the origin of the local coordinates Vg. Substituting Eqs. (8.5.21) and (8.5.22) 
into Eqs. (8.5.16) and (8.5.17) yields 



-2(1 + 7) f/<*'^>(a;,) - v<'U''^>{x,) + O(e^) , (/ > 1). (8.5.23) 



See Eq. (7.2.13). 



This makes the metric tensor and the scalar field in the local coordinates to be deter- 
mined by four sets of the independent external multipoles Pl, Ql, C^, and Z/,. However, 
the multipole moments are not physically significant and describe gauge-dependent 
coordinate effects. 

2) One removes from the metric tensor in the local coordinates all physically insignificant 
multipole moments Z^. The cost of this choice is a more complicated form of the time 
transformation, Eq. (7.2.11), from the local to global coordinates that involves now the 
infinite number of coefficients Bl. Assuming that all Zi^ = for / > 1 one obtains from 
Eqs. (8.5.16) and (8.5.17) 
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(8.5.25) 



At the present step of the matching procedure we prefer to keep the multipole moments Zl 
unspecified to preserve some freedom in making the gauge transformations. 



8.5.4 Matching of goQ{t,x) and gai3{u,w) to the order of 0{e^) 

Matching of the metric tensor at the post-Newtonian order of O(e^) allows us to infer the 
post-Newtonian equations of motion of the origin of the local coordinate system as well as 
the post-Newtonian corrections to the external multipole moments Ql and the remaining 
part of the time transformation formula (7.2.13) between the local and global coordinates. 
Expansion of the matching equation (7.2.4) for the metric tensor component, (700(^)^)1 to 
the post-Newtonian order, generalizes Eq. (8.5.1) 



N{u,w) + e^L{u,w) = N{t,x) + 2^{t,x) 



(8.5.26) 



L{t, x) - 3^^{t, x) + 22)(t, x) - 2<B(t, x)N{t, x) 



+ 4w|<B(t, x) + 2vlN{t, x) + 2viLi{t, x) + 2vi^i{t, x) 



- -vlml{t,x) + 2vlviy^At,x) 
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where the gravitational potentials in the right side of this equations are determined by Eqs. 
(4.2.2) - (4.2.5) and those in the left side are given by Eqs. (6.2.5), (6.2.6), (6.2.36), (6.2.39) 
and (6.2.40). Solution of Eq. (8.5.26) is done in several steps. 

First of all, one substitutes components of the matrix of the PPN coordinate transformation, 
given by Eqs. (7.3.7) - (7.3.10), to Eq. (8.5.26). Then, one matches separately the internal 
gravitational potentials referred to the central body B and those referred to the external 
bodies. The internal gravitational potentials have to be transformed from the local to global 
coordinates by making use of Eqs. (8.3.10) - (8.3.18). One notices that transformation equa- 
tion (8.3.17) for the Newtonian gravitational potential of the central body B can be written 
explicitly in terms of the functions coming about the matching procedure at lowest orders. 
Taking definitions of functions and from section 7.2 and elaborating them by making 
use of the results of previous steps of the matching procedure, one obtains explicit form of 



the relativistic correction, U(b){u,w), from Eq. (8.3.17) which describes transformation of 
the Newtonian potential from the local to global coordinates. It reads 



UiB^u, w) = W^\t, x) - -fU{xs) - alRl + (8.5.27) 

-\alxf{t,x)-vluUt,x) + 0{e'). 

Employing this formula in the matching procedure of the internal Newtonian potential of 
the central body B along with transformation equations for other internal potentials, one 
gets a remarkable result - both the Newtonian and the post-Newtonian terms depending 
on the internal structure of the body B cancel out and, hence, completely vanish from the 
matching equation (8.5.26). This effacing-internal structure effect can be explained in terms 
of the laws of conservation of intrinsic linear and angular momenta of the body B which are 
valid in the scalar-tensor theory of gravity as well as in general relativity [20,65]. In other 
words, the presence of a scalar field in the scalar-tensor theory of gravity does not result in 
the net self-force or self-torque exerted on the body which existence would violate classic 
Newton's laws. The internal-structure effacing principle for spherically-symmetric extended 
bodies was extrapolated to the 2.5 post-Newtonian approximation by one of the authors of 
the present paper [107,138] by applying the Fock-Papapetrou [48,49] method for derivation 
of relativistic equations of motion for a binary system with the conservative and gravitational 
radiation-reaction forces. Validity of the effacing principle was also confirmed by Damour 
[79] who derived equations of motion of two Schwarzschild black holes in the 2.5 post- 
Newtonian approximations by the Einstein-Infeld-Hoffmann technique [139] supplemented 
by the method of analytic continuation of the energy-momentum tensor of point-like particles 
understood in terms of the generalized functions (distributions) [137]. 

At the third step, one equates in Eq. (8.5.26) gravitational potentials generated by all external 
bodies (Moon, Sun, etc.) but the central body B (Earth). This step requires to know how 
function, X) jiQ l{u)w^ , is transformed from the local to global coordinates within the post- 
Newtonian accuracy. General formula of transformation of this function is given by Eq. 
(8.3.22). Substituting to this formula the explicit form of functions and displayed in 
Eqs. (7.2.6), (7.2.11) - (7.2.14), one gets 
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Matching equation (8.5.26) requires to calculate function x,oo(^5 generated by the external 
bodies other than the body B (Earth) and entering goo{t, x) component of the metric tensor 
in the global coordinates as shown in Eqs. (4.2.3) and (8.2.3). Contrary to other potentials, 
like U(t, x) , Wit, x) , $fc(t, x), which are solutions of the homogeneous Laplace equation in 
the vicinity of the body B, function x(t, x) is subject to the Poisson equation [48] 



V^x{t,x) = -2U{t,x) 



(8.5.29) 



After solving this equation and expanding its solution into STF harmonics (see appendix A) 
one obtains 
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Differentiating the left and right sides of Eq. (8.5.30) two times with respect to the global 
coordinate time, t, yields 



oo oo -j^ 



X 
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(8.5.31) 



Finally, one expands all functions in both sides of Eq. (8.5.26) in the Taylor series with 
respect to the distance, i?^, from the central body B, and reduce all similar terms. One 
finds that those terms which do not depend on Kg (that is, functions of time, t, only) form 
an ordinary differential equation of a first order for function B{t) in the coordinate time 
transformation given by Eq. (7.2.11). This differential equation reads 
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(8.5.32) 



Terms which are linear with respect to i?^ give us the post-Newtonian equation of transla- 
tional motion of the origin of the local coordinates, x^(t), in global (barycentric) coordinates. 
It generalizes equation (8.5.3) derived solely in the Newtonian approximation. Barycentric 
acceleration of the origin of the local coordinates with respect to the barycenter of the N-body 
system is 



F^'^Qk + <l,(^s) - + Q^|Y- 2q) 



(8.5.33) 
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This equation contains the external dipole moment, Qi, which must be calculated with the 
post-Newtonian accuracy in order to complete derivation of the post-Newtonian equation of 
translational motion of the origin of the local coordinates. A simple choice of = does 
not allow us to keep the origin of the local coordinates at the center of mass of the central 
body B (Earth) for arbitrary long interval of time. This is because the body is, in general, 
interacting with tidal gravitational field of external bodies (Moon, Sun, etc.) and does not 
move along a geodesic world line while the choice of = makes the origin of the local 
coordinates moving along a geodesic world line [45,69,94]. Thus, function Qi must be defined 
in such a way that the body's center of mass and the origin of the local coordinates coincide 
at any instant of time. This problem is equivalent to solution of the problem of motion of 
the body's center of mass with respect to the origin of the local coordinates and will be 
discussed in the next section. 

Terms which are quadratic, cubic, etc., with respect to the distance, i?^, determine the post- 
Newtonian corrections to the external STF multipole moments Ql = Qf + f^Qif ? where the 
Newtonian-like term QI is shown in equation (8.5.4). The post-Newtonian corrections are 

gf (t) =X<i> + Zl + $,<l>(^b) - lx,tt<L>{xB) + (lY - 2q'^Ul{xb) (8.5.34) 

+2(1 + 7)L^<^"^-i>(a^,) - 2(1 + ^)v';u''^\xb) 

+ (/ _ 27 - 2)v<''tf''^~'>{xB) + (1 + ^HULix^) 

J_^,k^,<Hu,L-i>k^^^^ + (2 - 2/3 - h)UixB)UL{xB) 

-{f -1 + 2 + 27)a<''f/'^-^>(a;s) - IF^<''U^^-^>\xb) , (/ > 2) 
where we used notations 

X<,,> = 3a<'ai> , (8.5.35) 

X<i> = 0, (/>3). (8.5.36) 

These equations finalize description of the STF multipole moments entering external solution 
of the metric tensor in local coordinates w"' = {cu, w) and the parametrized post-Newtonian 
transformation between the local and global coordinates. 



8.6 Final Form of the PPN Coordinate Transformation from Local to Global Coordinates 



For the sake of convenience we summarize the final form of the parametrized post-Newtonian 
coordinate transformation from local to global coordinates which is given by two equations: 
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Here functions A and B are solutions of the ordinary differential equations 
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while the other functions are defined as follows 
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These equations will be used in subsequent sections for derivation of equations of motion of 
extended bodies. They are also convenient for comparison with the relativistic transforma- 
tions derived by other researchers [1,15,43,44,45,61,69,70]. 



9 Translational Equations of Motion of Extended Bodies 



9. 1 Introduction 



In the Newtonian theory of gravity definitions of mass and the center of mass of an extended 
body, which is a member of a many-body system, are quite simple and straightforward 
concepts. Because of the simplicity, they have been directly extrapolated without any change 
to the relativistic theory of gravity by Fock [48] by making use of the invariant density p* 



which obeys the Newtonian-hke equation of continuity, Eq. (3.3.19). The invariant density is 
defined as mass of baryons per unit of the proper volume and after integration over volume 
of a body gives its total baryon mass which is constant [20]. The invariant density is used 
to introduce newtonian definitions of the body's center of mass, and its linear momentum in 
the body's local coordinates. The baryon mass, center of mass and the linear momentum of 
body B are given by integrals [48,127] 



These definitions were used by Fock [48], Papapetrou [49], Petrova [140], Brumberg [131], 
and some other researchers (see, e.g., [128,129,141,142] and references therein) for derivation 
of the post-Newtonian equations of translational motion of extended, spherically-symmetric 
bodies. Physical intuition tells us that equations of motion of such bodies, in principle, have to 
depend only on masses of the bodies which are supposed to be the only parameters character- 
izing the magnitude of their internal gravitational field. This is indeed true in the Newtonian 
theory. It was found, however, that the post-Newtonian equations of motion of spherically- 
symmetric bodies, unlike the Newtonian theory, depend not only on the baryon masses of the 
bodies, Eq. (9.1.1), but also on some other characteristics which are the internal kinetic and 
gravitational energy, elastic energy and moments of inertia [48,127,128,129,131,141,142,152]. 
Appearance of such terms significantly complicates interpretation of the post-Newtonian 
equations of motion and is unsatisfactory from the physical point of view. For instance, due 
to the dependence of the post-Newtonian equations of motion on other parameters rather 
than masses of the bodies, it is possible to argue that motion of a binary system consisting 
of ordinary stars is different from that of a binary system consisting of black holes having 
the same value of the mass as the stars. 

This point of view is incompatible with the Newtonian equations of motion of two black 
holes as it was shown by Demiansky and Grishchuk [125] who proved that motion of black 
holes in the Newtonian approximation obeys to the same laws of gravitational physics as 
for ordinary spherically-symmetric stars. In addition, the Einstein-Infeld-Hoffmann method 
of derivation of the post-Newtonian equations of motion operates with vacuum Einstein's 
equations and does not admit appearance of any terms in the equations of motion of spher- 
ical bodies which would depend on the internal structure of the bodies. Thus, one has to 
expect that the Newtonian definitions of mass, the dipole moment, and the center of mass 
given in Eqs. (9.1.1) - (9.1.3) are not quite appropriate for calculation of the post-Newtonian 
equations of motion of extended bodies. Indeed, our study of relativistic translational mo- 
tion of two spherically-symmetric compact stars in a binary system [107,108] revealed that if 
the Newtonian mass and the center of mass are replaced with their corresponding relativis- 
tic definitions all terms in the equations of motion depending on the internal structure of 
the bodies are effectively eliminated via renormalization of masses. Damour [79] called this 
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We skip in this section label B for all quantities referred to body B as it does not cause mis- 
interpretation. We shall label the bodies with indices A,B,C, etc. every time when it may cause 
confusion. 



property "the effacing principle" and confirmed its validity for spherically-symmetric and 
compact relativistic bodies using the matched asymptotic expansion technique. 

Newtonian theory predicts that if celestial bodies are not spherically-symmetric and rotate, 
their equations of motion must depend on additional parameters which are the mass mul- 
tipole moments of the bodies. It is natural to expect that the post-Newtonian equations of 
motion of such bodies will contain both the mass and current multipole moments given, for 
example, by Eqs. (6.3.1) and (6.3.8) respectively. However, it is not evident whether some 
other parameters have to appear in the relativistic equations of motion in addition to these 
two sets of the internal multipoles. Scrutiny analysis of this problem in general relativity elu- 
cidated that the post-Newtonian equations of motion of extended bodies do contain only the 
mass and current multipoles [69,70,71] and does not depend on any other internal character- 
istic of the bodies. PPN formalism operates with a class of alternative theories of gravity and 
does not obey to the "effacing principle" even for spherically-symmetric bodies because of 
the violation of the strong principle of equivalence [20,38,143,144]. This violation makes two 
masses - inertial and gravitational - be different and, hence, two parameters appear in the 
post-Newtonian equations of motion of spherical bodies as contrasted to general relativistic 
case. It is interesting to answer the question how many parameters have to be introduced to 
characterize motion of extended, non-spherically symmetric and rotating bodies in the PPN 
formalism, thus, extending general relativistic results of Damour, Soffel, and Xu [69,70,71]. 

Like in general relativity, solution of this problem in the PPN formalism could not be achieved 
in the framework of a standard a la Fock-Papapetrou post-Newtonian approach which basi- 
cally operates with a single (global) coordinate system. The global coordinates can not be 
used to define multipole moments of each body in N-body system in a physically consistent 
way. A local coordinate system should be constructed around each of the bodies in order 
to strongly suppress coordinate-dependent contributions to definitions of the multipole mo- 
ments caused by the Lorentz contraction and presence of the background gravitational field 
caused by external bodies. Consistent relativistic concept of the local coordinate system was 
developed in [45,61,69,106]. The concept of a local frame of reference associated with a mov- 
ing body has been also discussed in papers [41,42,43,44]. However, the authors of these papers 
always assumed that the center of mass of the body moves along a geodesic world line in the 
background space-time. This assumption is not valid for a non-spherical and rotating body 
as it follows from the Mathisson-Papapetrou-Dixon equations [49,75] and other arguments 
present in [45,69]. For this reason, the method of construction of the local coordinates in the 
vicinity of a self-gravitating body suggested in [42,43] is not general enough. Thorne [62], 
Blanchet and Damour [63] invented a post-Newtonian definition of the STF mass multipole 
moments which were fruitfully employed by Damour, Soffel and Xu [69,70] for derivation of 
the post-Newtonian translational equations of motion of self-gravitating extended bodies in 
general relativity. These authors had also proved that the only parameters present in the 
equations are Tolman masses of the bodies and their STF multipole moments referred to the 
local frame of reference of each body. 

In this section we derive the parametrized post-Newtonian equations of translational motion 
of extended bodies in the scalar-tensor theory of gravity and prove that these equations 
depend only on the inertial and gravitational masses of the bodies and a set of STF active 
multipole moments defined in section 6.3. In case of spherically-symmetric bodies the only 
parameters present in the equations are the inertial and gravitational masses which are 
different due to the Nordtvedt effect [143,144]. 



Post-Newtonian Thorne-Blanchet-Damour definition of tfie conserved mass and center of 
mass of a single isolated body are given by Eqs. (5.5.9) and (5.5.10). One might think that 
if the body is, in fact, a member of N-body system its gravitational interaction with other 
bodies of the system would violate these conservation laws [94]. The law of conservation of 
mass is indeed violated if the body is not spherically-symmetric (see below). However, the 
law of conservation of the center of mass and linear momentum of the body can be retained. 

There are three types of multipole moments in the scalar-tensor theory - active, conformal 
and scalar - which can be used for definition of body's mass and its center of mass. These 
moments were introduced in section 6.3 and are defined by equations (5.2.14), (5.3.6) and 
(5.4.6). By direct calculation we shall demonstrate that the active and scalar dipole moments 
of the body are not efficient in derivation of the translational equations of motion. This is 
because if one uses either active or scalar dipole and derive equations of motion of the body 
in its local coordinates the equations, we shall obtain, will contain a significant number of 
terms which can be treated as self-accelerations and they can not be removed by simple 
translation of the origin of the local coordinate system to another point. Self- accelerated 
terms in the equations of motion violate Newton's third law and are unacceptable. On the 
other hand, when we use the conformal dipole moment for defining the body's center of 
mass and, then, derive equations of motion of the body, we find that self-acceleration terms 
do not appear, and the equations have remarkably simple structure of the post-Newtonian 
force which is a function of the active multipole moments of the body coupled with external 
multipoles. 

Let us discuss derivation and specific of the translational equations of motion in more detail 
starting from the explicit form of the local (macroscopic) equations of motion of the body's 
matter in local coordinates. 



9.2 Macroscopic post-Newtonian Equations of Motion 



The macroscopic post-Newtonian equations of motion of matter consist of three groups: (1) 
the equation of continuity, (2) the thermodynamic equation relating the elastic energy 11 
and the stress tensor Ha/s, and (3) the Euler equation. 

The equation of continuity in the local coordinates {u, w) has the most simple form for the 
invariant density p* and reads 

fl^^.o. ,9.2.1) 

This equation is exact, that is takes into account all post-Newtonian corrections, as follows 
from the definition of the invariant density p* and Eq. (3.3.19). 

The thermodynamic equation relating the internal elastic energy 11 and the stress tensor 
TTais is required only in a linearized approximation where the stress tensor is completely 
characterized by its spatial components vTjj. Hence, one has from Eq. (3.2.5) the following 
differential equation 

where the operator of convective time derivative is d/du = d/du + u^d/dw^. 



The Euler equation follows from the spatial part of the law of conservation of the energy- 
momentum tensor T*'^^ = 0. It yields 
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where gravitational potentials A^, L, Lj are the metric tensor components in the local 
coordinates and they have been defined in section 6.2. 



9.3 Definitions of Mass, the Center of Mass and the Linear Momentum of an Extended 
Body in the N-hody System 



There are two algebraically independent definitions of the post-Newtonian mass in the scalar- 
tensor theory - the active mass and the conformal mass which are derived from Eqs. (6.3.1) 
and (6.3.5) respectively for index Z = 0. As discussed in section 6.3 one must retain contribu- 
tion of the gravitational potential of external bodies in the definition of the STF multipole 
moments of the body's gravitational field. It will allow us to cancel out in equations of 
motion all terms depending on the internal structure of the central body B which are not 
incorporated to the definition of the STF multipoles. Absence of such terms in equations of 
motion extends validity of the effacing principle from general theory of relativity [69,79,107] 
to the scalar-tensor theory of gravity at least in the first post-Newtonian approximation. The 
question about whether to keep the contribution of the gravitational potentials of external 
bodies in the definition of the STF multipolar decomposition of the body's gravitational 
field was discussed by Thorne and Hartle [94] but they did not come up with a definite an- 
swer. Our approach is based on direct calculation of equations of motion and we have tried 
various possible definitions of the center of mass and the STF multipoles. After tedious and 
cumbersome calculations we came to the conclusion that the equations of motion have the 
most simple form and the minimal set of parameters if we take the conformal definition of 
mass and the center of mass for each body and include the gravitational potential of external 
bodies to the definition of the body's multipole moments 

Con/orma/ multipoles are given by Eq. (6.3.5) and for index I = can be reduced to simpler 
form such that the conformal mass of the body located near the origin of the local coordinates 
is 

M = M-e'}^[Y + {l- 7)P] M + -^Q^^j + ^(^') • 

In what follows, we will need definition of the active mass of the body as well. It is extracted 
from Eq. (6.3.1) when index / = 



In fact, this is the Navier-Stokes equation because the stress tensor is taken in its the most 
general form. 
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is general relativistic definition of tlie post-Newtonian mass of the body [20] and 
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J(2) = / p*w^Sw , 

JVb 

is the second-order rotational moment of inertia of the body. 

It is not difficult to derive a relationship between the active and conformal post-Newtonian 
masses by making use of Eqs. (9.3.1) - (9.3.3). It reads 
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where 77 = 4/? — 7 — 3 in Eq. (9.4.1) is called the Nordtvedt parameter [20]. Numerical value 
of this parameter is known with the precision better than 0.02% from the lunar laser ranging 
experiment [145] which lasts already for more than 30 years. 

One can see that in the scalar-tensor theory of gravity the conformal mass of the body differs 
from its active mass. If the body is completely isolated the difference can be only due to the 
Nordtvedt eSect, that is for 7^ 0, and the time-dependence of the body's rotational moment 
of inertia (for example, radial oscillations, etc.). In the case when the presence of external 
bodies can not be ignored, one has to account for coupling of the external gravitational field 
multipoles, Ql and Pl, with the internal multipole moments of the body. 

It is important to realize that in general case the general relativistic post-Newtonian mass 
of an individual body is not conserved. Indeed, taking a time derivative of Eq. (9.3.3) and 
making use of the macroscopic equations of motion of the body's matter given in Section 
9.2, one gets [129] 
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where overdot means the time derivative with respect to the local coordinate time u. This 
equation reveals that the general relativistic mass of the body is constant, if and only if, 
the mass multipole moments Tl of the body do not depend on time and/or there is no 
external tidal fields, that is = 0. However, one can notice that the conformal and active 
post-Newtonian masses are not constant in the presence of the tidal field even if the body's 
multipole moments Tl are constant. This is because the external multipole moments Ql 
enter definitions of these masses, Eqs. (9.3.1) and (9.3.2), explicitly, and in general case of 
N-body problem they depend on time. 



Direct calculation of the equations of motion elucidates that definition of the conformal 
mass dipole moment given by Eq. (6.3.5) for / = 1 gives the most optimal choice of the 
post-Newtonian center of mass for each body. This is because after differentiation with 
respect to time only the conformal dipole moment leads to the law of conservation of the 
body's linear momentum when one neglects the influence of other external bodies, while 
the post-Newtonian scalar or active dipole moments do not bear such property. Thus, the 
post-Newtonian center of mass of the body, JT"* = X*, is derived from Eq. (6.3.5) for / = 1 
and reads 
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denotes a new STF object [1]. We call attention of the reader to the fact that for / 
scalar function Af = T^'^\ where X^^^ has been defined in Eq. (9.3.4). 



the 



It is worth noting that the post-Newtonian definitions of mass and of the center of mass of the 
body depend not only on the internal distribution of matter's density, velocity, and stresses 
inside the body but also on terms describing the coupling of body's gravitational field with 
that of external masses. As we shall show later, inclusion of these coupling terms in definitions 
(9.3.1), (9.3.2) and (9.3.7) is absolutely necessary in order to simplify translational equations 
of motion as far as possible and bring them to the form which can be reduced to the Einstein- 
Infeld-Hoffmann equations of motion in the limiting case of spherically symmetric bodies. 
In this sense, the question about whether the coupling of internal and external gravitational 
fields should be included in the definitions of the post-Newtonian mass and the center of 
mass, which was a matter of concern for Thorne and Hartle [94], can be considered as having 
been resolved. 



The post-Newtonian linear momentum of the body, V\ is defined as the first time derivative 
of the dipole moment given by Eq. (9.3.7), that is V' = where dot indicates a derivative 
with respect to the local coordinate time u. After taking the derivative one obtains 
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where function 
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Until now the point XB{t) represented location of the origin of the local coordinate system 
in the global coordinates taken at the time t. In general, the origin of the local coordinates 
does not coincide with the center of mass of the body which can move with respect to the 
local coordinates with some velocity and acceleration. In order to be able to keep the center 
of mass of the body at the origin of the local coordinates one must prove that for any instant 
of time the dipole moment defined by Eq. (9.3.7) and its time derivative (that is, the linear 
momentum of the body) given by Eq. (9.3.9) can be made equal to zero. This requirement 
can be achieved, if and only if, the second time derivative of the dipole moment with respect 
to the local coordinate time u is identically equal to zero, that is 

V' = Q (9.3.11) 



It is remarkable that this equation can be satisfied after making an appropriate choice of 
the external dipole moment Qi that characterizes a locally measurable acceleration of the 
origin of the local coordinates with respect to another local coordinate frame whose origin 
moves along a geodesic world line in the background space-time. This statement has been 
proven in [45] in the Newtonian approximation and, then, extended up to the first general 
relativistic post-Newtonian approximation in [69]. In the present paper we shall derive the 
consequences of Eq. (9.3.11) in the first post-Newtonian approximation of the scalar-tensor 
theory of gravity characterized by two PPN parameters, (3 and 7. 



9.4 Translational Equation of Motion in Local Coordinates 



Translational equation of motion of the body in the local coordinates, = (cm,iu), is 
derived by making use of the definition of its conformal linear momentum, V\ displayed in 
Eq. (9.3.9). Differentiating Eq. (9.3.9) one time with respect to the local coordinate time u, 
operating with the macroscopic equations of motion, Eqs. (9.2.1) - (9.2.3), and integrating 
by parts to re-arrange some terms, one obtains 
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where we have shown exphcitly all terms proportional to Qi, rj = 4/3 — 7 — 3, and the 
post-Newtonian correction AP* is given by 
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It is worth noticing that mass Ai is the active mass, and the STF muhipole moments 
(/ > 1) of the body, which appear in the right side of Eqs. (9.4.1) and (9.4.2), are the active 
mass muhipole moments depending on time u. Function Qi in Eq. (9.4.1) has not been yet 
restricted and can be chosen arbitrary. Its choice determines a world line of the origin of the 
local coordinates. If one chooses Qi = 0, then the origin of the local frame, Xg, is moving 
along a geodesic world line defined in the global frame by Eq. (8.5.33) and the center of mass 
of the body under consideration is moving in accordance with the law of motion, Eq. (9.4.1), 
with respect to this geodesic. For practical applications, however, it is more convenient to 
chose the origin of the local frame to be always located at the center of mass of the body. This 
can be accomplished by imposing condition (9.3.11), that is "P* = = 0. If this condition 
is satisfied it allows us to chose JT* = = exactly, which leads to similar conditions for 
the active dipole moment, X* = X* = X* = 0, in the Newtonian approximation. However, in 
the post-Newtonian approximation the active dipole moment of the body, X* 7^ 0, even if the 
conformal multipole moment of the body, JT* = 0, because these two moments are defined 
via different equations (see Section 6.3). 

Fixing the origin of the local frame at the body's center of mass {J"^ = 0) and noticing 
that Ql = Pl for any / > 2 makes the post-Newtonian function AP* = 0. Then, one can 
find solution of Eq. (9.4.1) for function Qi which was considered until now as arbitrary 
variable. We remind that the physical meaning of Qi is acceleration of the origin of the 
local coordinates with respect to a geodesic world line in the background space-time defined 
by the gravitational field of external bodies. Equating P*, separating acceleration Qi of the 
body's center of mass from all other terms in Eq. (9.4.1) and taking all terms with Qi to the 
right side yields the following post-Newtonian equation 



M,,Q, = Fjv + (f;^ + AF;^) + O(e^) , (9.4.3) 
where the conformal anisotropic tensor of mass 
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and the (tidal) gravitational forces 
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Eqs. (9.4.3) - (9.4.7) describe the law of translational motion of the body in the local 
coordinates in the presence of external bodies which create a force dragging motion of the 
body's center of mass from geodesic world line. Newtonian, F^, and the post-Newtonian, 
Fpjv, tidal forces are caused by gravitational coupling of the body's internal (active) multipole 
moments, 2l and Sl, with the external multipole moments, Ql and Cl- The post-Newtonian 
tidal force, F^^y, is reduced in the limit of 7 = 1 to general relativistic expression derived 
previously by Damour, Soffel and Xu [69]. 

It is worthwhile to emphasize that summation with respect to index / in Eq. (9.4.5) begins 
from / = 1. The point is that we have defined the center of mass of the body B in terms of 
the con/orma/ dipole moment J^'^ by the condition J'^ = 0. However, the force depends on 
the active multipole moments of the body but the active dipole moment X* 7^ J'^ and, hence 
X* 7^ 0. For this reason, one has to take into account the contribution to the force F^y coming 
out of the non-zero active dipole of the body. This contribution has a post-Newtonian order 
of magnitude and can be written down in more explicit form as 
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where the active dipole moment X* is 
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We remind that the conformal dipole moment of the body, ^7* = 0. Hence, Eq. (9.4.8) is, in fact. 
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If one takes into account explicit relationships between the multipole moments of the 
scalar field and gravitational potential of the external bodies, then Eq. (9.4.9) can be slightly 
simplified 
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It is clear that dipole moment of the central body can contribute to the equations of motion 
of the body only in the scalar-tensor theory of gravity because in general relativity /? = 7 = 1. 
It is important to notice the presence of the Nordtvedt parameter 77 = 4/3 — 7 — 3 in the active 
dipole moment of the body. It is likely that we observe a general feature of the relativistic 
equations of motion in the scalar-tensor theory. It looks like each active multipole moment 
of the body under consideration has a contribution being proportional to the Nordtvedt 
parameter 77. This, for example, leads to inequality between inertial and gravitational masses 
of the body in case of the monopole mass moment. Presence of the Nordtvedt parameter- 
dependent term(s) in dipole moment is insignificant for the bodies whose shape is close to 
spherically-symmetric. However, it may play a role in motion of sub-systems, like that of 
Earth and Moon, which possess large deviations from spherical symmetry. Experimental 
study of this problem would be desirable. 

Gravitational force AF^^ is also essential only in the scalar-tensor theory of gravity. It is 
proportional to the dipole moment of the scalar field Pi and its time derivatives. Had the 
scalar field absent the dipole moment Pi of the scalar field could not exist and the force AF^^ 
would be zero. The dipole moment of the scalar field Pi couples with the self-gravitational 
energy ~ / p*U^^^ of the body as well as with the energy of external gravitational field and 
kinetic energy of the body's internal motion . In the next section we shall show that it 
is this coupling that is responsible for inequality of inertial and gravitational masses of the 
body and should be treated as a violation of the strong principle of equivalence. 



a difference between the active and conformal dipole moments of the body, X* — J'^ 
^'^ Observe the term T^^^ 



9.5 Equation of Translational Motion in Global Coordinates 



Equation of translational motion of the body's center of mass in the global coordinates 
= (ct, a:;) can be obtained from the equation of motion, Eq. (8.5.33), of the origin of the 
local coordinates, Xb, if the acceleration Qi is subject to obey to the local equation of motion. 
(9.4.3). The acceleration Qi depends on scalar multipoles Pl so that after substitution Eq. 
(9.4.3) to Eq. (8.5.33) one has to make use of Eq. (8.4.1), that is Pl = U^l{xb)- Moreover, one 
has to use Eq. (8.5.3) to replace in forces AF^v, F^^y, and AF^jy all terms depending explicitly 
on Qi with a linear combination of the barycentric acceleration and the gradient U^ii^Xg) of 
the gravitational potential of external bodies, Qi = U^ilxg) — a^. Taking all terms depending 
on the gravitational gradient U^ilxg) to the right side of Eq. (8.5.33) and those depending 
on the barycentric acceleration to the left side, bring Eq. (8.5.33) to the following form 
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where the external potentials (/{xg), U\xb), '^{xg), and x{xb) are defined in Eq. (8.2.5) 
and are taken on the world line of the center of mass of the body B. The external potentials 
can be expanded in multipolar series so that the translational equation of motion (9.5.1) will 
depend only on the active multipole moments of the bodies. We do not present this general 
result in the present paper but consider a more simple case of spherically-symmetric and 
rotating bodies in section 11. 



One has to notice that inertial mass M-g oi the body B is its conformal mass. It is not 
equal to gravitational mass M.g of the body in the right side of this equation which is its 
active mass. Difference between the two masses is given by Eq. (9.3.5) and it causes violation 
of the strong principle of equivalence for massive extended bodies. Existence of the possible 
difference between the inertial and gravitational masses in alternative theories of gravity was 
pointed out by Dicke [36,37,38] and Nordtvedt [16,17,145]. In our calculations the inertial- 
gravitational mass difference originates from Eq. (9.4.7) which has terms proportional to 
scalar dipole. Pi = U^i{xg), and they contribute to the gravitational (active) mass only. 



Forces and F^^ are given by Eqs. (9.4.5) and (9.4.6). Terms in the first and second curled 
brackets of Eq. (9.5.1) being proportional to mass Aig are the post-Newtonian corrections 
to the Newtonian force acting on the body B considered as a monopole massive particle . 
The group of terms in third and forth curled brackets in Eq. (9.5.1) represents the post- 
Newtonian correction to the Newtonian tidal force F^ and takes into account higher order 
multipoles of the body B. In particular, these terms contain time-dependent functions Q 
and Y which define the unit of time and length in the local coordinates. This correction also 
contains the matrix of relativistic precession F*'^ given in Eq. (8.6.8). Eq. (9.5.1) describes a 
generic case of translational equation of motion of extended bodies having arbitrary shape 
and rotation (all multipoles). We derive rotational equation of motion of the body B in the 
next section. 



10 Rotational Equations of Motion of Extended Bodies 

Rotational equations of motion for each body define orientation of the body's angular 
momentum (spin) at each instant of time with respect to the local frame of reference 
w° = (cm, w) which axes are not dynamically-rotating, that is Fermi- Walker transported 
[76] in accordance with Eq. (8.6.8) describing orientation of the axes of the local coordinates 
with respect to the global coordinates at each instant of time. We shall work out relativistic 
equations of the rotational motion by making use of the method proposed in [71]. 

10.1 Post- Newtonian Definition of the Angular Momentum of the Body 

First of all, one needs to introduce definition of the angular momentum of an extended body 
from the N-body system. In principle, one had to use the same principle of the multipolar 
expansion of the metric tensor applied to its goi{u, w) component calculated with taking into 
account all terms of the next post-Newtonian approximation. This procedure was applied 
by Damour and Iyer [64,91] for the post-Newtonian definition of angular momentum of an 
isolated system. It is not known yet how to apply the Damour-Iyer procedure to a single 
body from the N-body system because of its complexity. Therefore, we shall use the approach 
proposed in [71] to bypass this difficulty. 

Let us introduce a bare spin of a single body B by making use of the following post-Newtonian 
definition 

= -^j EijkW^Q'^^d^w , (10.1.1) 

where 

eo^- = [-g)<p (^T°^- + , (10.1.2) 

is a linear combination of tensor of energy-momentum of matter T^^ and the pseudo-tensor 
of gravitational field t^" . We have defined the bare spin of the body B by Eq. (10.1.1) because 
it corresponds to the conserved spin of an isolated system (see Eq. 5.5.6 for more detail). We 
also assume that the center of mass of the body B is chosen such that its conformal dipole 
moment is zero. 

Integration in Eq. (10.1.1) is formally performed in local coordinates over entire space. How- 
ever, tensor of energy-momentum T^^ includes the matter of the body B only and depends 



on the complete metric tensor g^^ in the local coordinates. We assume that the pseudo-tensor 
t^^ depends only on the internal part of the local metric tensor for the body B. Integration 
by parts allows us to reduce the bare spin of the body B in Eq. (10.1.1) to the following 
expression 
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where integration is over the volume of the body B, and potential W^^^' is defined by Eq. 
(9.3.10). We shall use Eq. (10.1.3) to derive rotational equations of motion of body's angular 
momentum. 



10.2 Equations of Rotational Motion in Local Coordinates 



Rotational equations of motion for body's spin are derived by differentiation of Eq. (10.1.3) 
with respect to the local coordinate time u. After taking the time derivative and making use 
of the macroscopic equations of motion in local coordinates given in section 9.2, one makes 
several transformations of the integrand to reduce similar terms and to simplify final result. 
After tedious but straightforward calculations done in the spirit of paper [71], one obtains 
equations of the rotational motion of body B in its own local coordinate frame 




AT' - -^AS' 
du 



(10.2.1) 



where T* is a general- relativistic torque for 7 = 1, and AT* is its post-Newtonian correction 
due to the presence of the scalar field, while AiS* can be considered as a supplementary 
post-Newtonian contribution to the bare spin 5*. The torque and other terms in the right 
side of Eq. (10.2.1) read as follows: 
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is additional set of multipole moments which has been used aheady in definition of the 
multipole moments in section 5. 

General relativistic torque T* depends on the multipole moments Zi which define the residual 
gauge freedom. They can be used to simplify the post-Newtonian correction to the torque, 
AT*. This correction is, in fact, exactly equivalent to AT* = Sij^a^g (X* — J"^), where X* 
and JT* are active and conformal dipole moments of the body B respectively. The difference 
between the two dipole moments taken under condition that = is given by Eq. (9.4.11) 
and has been reproduced in Eq. (10.2.3). We have taken into account the external monopole 
moments Q and Y defining the units of measurement of the local time and spatial coordinates 
respectively. Their contribution is to the rotational equations of motion is extremely small 
and can be omitted. 



We re-define the spin of the body as 

so that equations of rotational motion acquire their final form 

d 



du 
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(10.2.7) 



These equations should be compared with analogous equations derived by Klioner and Soffel 
[1] First of all, we notice that our definition of the multipole moments 
differs by a numerical factor (/ + 1)// from that, C'Kii°ii°>-Sofrei^ ^^^^ Klioner and Soffel, 
that is 

-2(1+7)'^^ 



j^Kopeikin — Vlasov 



^Klioner— Soffel 



10.2.8 



Comparison of spins, 5^ - our notations, and, 5"* - Klioner-Soffel's notation, shows that 
they are equal if sign minus in front of three last terms in Eq. (9.45) of Klioner-Soffel's 
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These equations are numbered (9.42) - (9.47) in [1]. 



paper [1], is replaced with sign plus. The general relativistic torque in the present paper 
coincides naturally with that derived in [1]. The biggest difference occurs between the post- 
Newtonian correction AT* to the torque in this paper and a corresponding quantity given 
in two equations, (9.43) and (9.46), in the paper [1]. First three terms in our Eq. (10.2.3) 
completely coincides with equation (9.43) derived by Klioner and Soffel, thus, confirming the 
presence of the Nordtvedt effect for rotational motion of the bodies. However, we obtained 
different terms in the second and third lines of our Eq. (10.2.3) as contrasted with equation 
(9.46) of Klioner and Soffel's paper [1]. We suppose that the difference may come out of 
slightly different gauge conditions used in this paper and in [1]. Additional origin of the 
difference is that we used definition of the center of mass of the body which is not reduced 
to that used by Klioner and Soffel. 



11 Motion of spherically-symmetric and rigidly-rotating bodies 

11.1 Definition of spherically- symmetric and rigidly rotating body 

It is well understood that the notion of a spherically-symmetric and rigidly rotating body 
is not invariant but coordinate-dependent [66,76]. According to special theory of relativity, 
coordinate grid of a moving coordinate frame is linearly deformed and magnitude of this 
deformation depends on velocity of the frame with respect to a reference frame being at 
rest. This consideration assumes that if one considers a spherically-symmetric body in a 
static frame it will be not spherically-symmetric in a moving frame. Deformation of the 
body's shape can be calculated by applying the Lorentz transformation to the equation 
describing the shape of the body in the static frame [146] . The Lorentz deformation is solely 
coordinate effect which does not lead to appearance of physical stresses (tensions) inside the 
moving body. Nonetheless, the Lorentz deformation of the body's shape has to be taken into 
account for correct calculation of observed physical effects associated with motion of the 
body. Poincare and Lorentz were first who took into account special relativistic deformation 
of a moving electron for calculation of the electromagnetic radiation-reaction force exerted 
on the electron due to the emission of electromagnetic radiation [147]. 

In general relativity gravitational field causes deformation of coordinate grid of a static 
frame with respect to the grid of the same frame taken in the absence of gravitational 
field . Hence, this distortion of the coordinate grid represent pure mathematical effect 
and does not cause physical deformation of the body being at rest with respect to this 
frame. However, gravity-caused deformations of the local coordinate's grid must be taken into 
account in calculation of translational equations of motion of the body with respect to the 
global coordinate frame. It is worth mentioning that one has to distinguish the mathematical 
deformations of the frame from the physical (tidal) deformations of the body itself [72] . This 
can be achieved if a precise relativistic theory of reference frames is employed. The PPN 
formalism brings new complications due to the presence of additional fields which can cause 
both coordinate and physical deformations of the body's volume. 

Post-Newtonian definition of the multipole moments of the gravitational field is coordinate- 



This represents pure mathematical comparison. It is not physically possible to turn off or to 
screen gravitational field. 



dependent. Therefore, the exphcit structure of the multipolar expansion and the number 
of terms present in this expansion crucially depend on the choice of coordinates. Transfor- 
mation from local coordinates to global ones will change mathematical description of the 
multipole moments and we have to be careful in finding the most precise formulation of 
the notion of spherically-symmetric body to avoid introduction of non-physical multipole 
moments. Any misunderstanding of this concept will lead to inconsistencies in calculation of 
equations of motion for N-body problem in the first [107] and higher-order post-Newtonian 
approximations [148] and/or appearance of spurious coordinate-dependent terms having no 
physical meaning. 

We assume that for each body of N-body system the geometrical center of the body's spher- 
ical symmetry is located at the center of mass of the body that coincides with the origin of 
the local coordinates associated with this body. We assume that all functions characterizing 
internal structure of the body have spherically-symmetric distribution in the local coordi- 
nates. These functions are: the invariant density p*, the internal energy 11, and the stress 
tensor TTjj. Spherical symmetry in the local coordinates means that these functions depend 
only on the local radial coordinate r = \w\: 



Moreover, we assume that the internal distribution of matter does not depend on the local 
coordinate time u that excludes radial oscillations of the body from consideration. Radial 
oscillations can be easily included in our version of the PPN formalism but we postpone 
treatment of this problem for future work . 

Spherically-symmetric distribution of matter must generate a spherically-symmetric gravi- 
tational field. Therefore, the multipolar expansion of the Newtonian gravitational potential 
of the body must have in the local coordinates only a monopole term 



where the baryon (Newtonian) mass Ai^^B is defined in Eq. (9.1.1). Strictly speaking, this 
monopole expansion will be violated at some order of approximation because the external 
tidal force of the background gravitational field acts on the body and deforms its spherically- 
symmetric distribution of matter. This tidal distortion is proportional to some numerical 
coefficient which characterizes elastic properties of the body under consideration. Assum- 
ing that the body is made of matter with sufficiently low elasticity one can reduce the tidal 
deformation of the body to a negligibly small value, at least in the first post-Newtonian 
approximation . 

We shall consider the case of rigidly rotating bodies for which the internal velocity of matter 
(as defined in the local coordinates) is a vector product of the angular velocity Q^, referred 
to the local frame, and the radius-vector w*, that is 



" Notice that radial oscillations are irrelevant for consideration of this problem in general relativity 
due to the Birkhoff 's theorem [66,76] 
'^^ Love's number /c2 [68] . 

Analytic estimate of the magnitude of the tidal deformation comparatively with the magnitude 
of the post-Newtonian forces has been done in [79,107]. 



p*{u,w) = p*{r) , 



U{u,w) = n(r) , 



(11.1.1) 




(11.1.2) 



v^=ey,niw^ (11.1.3) 

u' = 'kfy-ninlw<^'> , (11.1.4) 

in the local frame of the body B. Again, one should remind that rotation causes rotational 
deformation of the body and distorts spherical symmetry. However, the rotational deforma- 
tion is proportional to the same Love's number, /c2, [68] and by assuming that the body is 
rigid enough and rotates sufficiently slow, one can make the rotational deformation to be 
negligibly small. In what follows we shall use this assumption and neglect the rotational 
deformation. 

Spherical symmetry of each body assumes that one can use the following (pure geometric) 
properties which are valid for any function /(r), depending on radial coordinate r only [62]: 



f{r)w^-^--^-^d^W = ^ 5(a.a....5a._,a,, /^^ fivy'd'w , (11.1.5) 

f{r)w'^'^-'^'+'d^w = 0, (11.1.6) 

Vb 

where 5(aia2---'5a2i_ia2i) is the fully symmetric linear combination of the Kronecker delta sym- 
bols [62]. In particular, for any / > 1 one has 

/ f{r)w<'''^-''>d^w = . (11.1.7) 

JVb 



One will also need several other equations for performing integration over sphere in the local 
coordinates of body B. They are as follows: 



A<.L>B^N> p*w<'^>w<''>d'w=\ (2/ + 1)!! B ' ^ ^ ^^^^^^^ 



Vb 
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[ , (/ > 2) 

where we used Eq. (11.1.4), A^l^ and -B<l> are arbitrary STF tensors, and 

Jf) = / p*r^^d^w, (11.1.10) 

JVb 



is 2/-th order rotational moment of inertia of the body B 
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Eqs. (11.1.5) - (11.1.10) will be used for calculation of multipolar expansions of various 
gravitational potentials entering translational equations of motion of the bodies. 



Notice that all odd rotational moments I^'^'-'^^^ = 0. 



11.2 Coordinate Transformation of Multipole Moments 



Multipolar expansion of the Newtonian potential in the global coordinates, = {ct,x), 
introduces multipole moments of a body defined as integrals over the body's volume taken 
on hypersurface of constant time t, that is 

1^^= [ p*{t,x)R'^R'^...Rld^x , (11.2.1) 

JVb 

where i?^ = — x^, and is the origin of the local coordinates coinciding with the center 
of mass of the body. We have postulated that the density and other structure-dependent 
functions inside the body have a spherical-symmetric distribution in the local coordinates 
= (cM,tu), so that according to Eq. (11.1.7) the following relationship must held for any 
/ > 1 

/ p*(r')tf'<^>rfV = , (/ > 1) (11.2.2) 

JVb 

where the integration is over a hypersurface of constant local coordinate time u. However, 
Eq. (11.2.2) does not assume that the multipole moments of the body B defined in the global 
coordinates, I^, equal zero, in fact 7^ for / > 1. One can calculate directly from Eq. 
(11.2.2) after making use of transformation formula, Eq. (8.3.6), from the local to global 
coordinates 



w 
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where R'l = x'* 



™i pi 



X* — X* , v'^ 



dx'^/dt, = dx^/dt, and 



F,, = -e,,kJ''' , (11.2.4) 
Dij=^5ijU{xB) , (11.2.5) 

D^Jk = I (^ai6'' + a';6'^ - a'J^"^ , (11.2.6) 
and function J-''' is defined in Eq. (8.5.20). 

Eq. (11.2.3) must be used for transforming integrals shown in Eq. (11.2.2) from the local 
to global coordinates. This is because the integration in Eq. (11.2.2) is performed over the 
hypersurface of constant (local) time coordinate u while similar integrals in the multipolar 
decomposition of the Newtonian gravitational potential in the global coordinates are defined 
on the hypersurface of constant (global) time coordinate t. Transformation of space coordi- 
nates from the space-like hypersurface of constant time u to that of time t depends on space 
coordinates, x, of the point at which matching of the local and global coordinates is done so 
that the two space coordinates: x' - the point of integration, and, x - the matching point, 
appear in Eq. (11.2.3) simultaneously and they both belong to the hypersurface of constant 
time t. 



^ This kind of transformation of the multipole moments turns out to be important for calculation 
of the 3-d post-Newtonian equations of motion [148]. 



Substitution of Eq. (11.2.3) into Eq. (11.2.2) yields [138] 



/ p*{u,w')w'<'^>d'w' = I<'^> + ( -viv<n^^~'>^ - IF^<H'^~'>^ (11.2.7) 
Jvb V2 

Taking into account Eq. (11.2.2) one concludes that only the dipole, I^, and the quadrupole, 
I^, moments differ from zero. More specifically, one has 



K = [e^^kv^nl + la^) + O [e') , (11.2.8) 

lf^> = -^T(f)^;<X> + O {e') , (11.2.9) 



l<L> = 0(e^] , (/ > 3) (11.2.10) 



The same expressions for the multipole moments can be obtained in a different way by 
making use of multipolar expansions of the Newtonian potential of body B in the local 
and global coordinates and their subsequent comparison with the help of transformation 
formula shown in Eq. (8.5.27). We have checked that both derivations are self-consistent 
and yield identical expressions for multipole moments given in Eqs. (11.2.8) - (11.2.10). 
Transformation of multipole moments of Earth's gravitational field from the global to local 
coordinates were used in our papers [60,61] in order to derive the post-Newtonian equation 
of motion of Earth's artificial satellites (valid also for Moon) in the geocentric frame with 
taking into account relativistic corrections due to the presence of Earth's quadrupole field. 



11.3 Multipolar Decomposition of Gravitational Potentials in Global Coordinates 

In order to derive equations of motion of the bodies in global coordinates one will need to 
know the multipolar decomposition of gravitational potentials f/fl(t, x), U^g-^it, x), $(fl)(t, x), 
and X(B)(t, x) in these coordinates. The potentials under discussion are defined in Eqs. (4.2.8) 
- (4.2.14). For the Newtonian potential one has 

where the dipole, P, and quadrupole, I^*-'^, moments are given by Eqs. (11.2.8) and (11.2.9). 



Vector-potential UJ^g^ {t, x) is decomposed as follows 

uu, ii;) + o (.^) . (11,3,2) 

where we have used the fact that inside the body, = + + O (e^), and the internal 
velocity, i^*, is defined in Eq. (11.1.3). 

Superpotential X(B){i,'x) has the following multipolar decomposition 



1 cx'^^ 

X,,,{t,x) = -GM.sRs--^ + 0{e') . (11.3.3) 



Potential $(B)(t, x) consists of a linear combination of four functions as shown in Eq. (4.2.6). 
For each of these functions one has 



J \X X \ XL R IL R J Vr 



p*{t,x')v'^£x' _GM,bvI , G_ 

Vb i^' ^ ^' I J^B Rb ■'Vb 

2G eijkR'B<^llB^ G n<'n^>R%RiI^B^ 



O e 



R% 5 R% 



JVb \x — x'\ Kb JVb 

°^ ^ 1)' )i?<i> 5W 1 ^ ^ . 2 



, V ±j J-i*A-^B ^^BA ^ ^ \ \ n r 



^ / p-(t.a|-)n(t,x-)fe- G ^ ,.MnMA,H-o(.') , (11,3,6) 

J^s |a; — X'\ Kb JVb ^ ' 

^rit,x) = G[ =^ f p^r)d'^ + o (e^) . (11.3.7) 

JVb \x — x'\ Rb JVb ^ ' 

This concludes the set of equations describing the multipolar decomposition of the gravita- 
tional potentials in the global coordinates. 



11. Translational Equations of Motion 



Both conditions of spherical symmetry and rigid rotation, Eqs. (11.1.1) and (11.1.3), allow 
us to simplify equation (9.5.1) of translational motion of body B drastically. For example, 
Eq. (11.1.7) assumes that all multipole moments of the body 1^ = (e^) for all / > 1, and 
= O (e^) for all / > 2. Therefore, calculation of the tidal Newtonian force for body B, 
taken from Eq. (9.5.1), yields 
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which has the post-Newtonian, (~ e^), order of magnitude. 



The post-Newtonian gravitomagnetic tidal force from Eq. (9.5.1) in case of a spherically- 
symmetric body is reduced to 
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B is the mass, and 



is the spin of the body B. 



5; = ^x(f)fi^, 



;ii.4.4) 



All other terms in Eq. (9.5.1) depending on are equal to zero. Hence, Eq. (9.5.1) is 
drastically simplified for spherically-symmetric bodies and reads 



Msg's = MsViixs) - F 



iSikvl - f >^ - 2(7 + p)SikU (xs) 



U.kiXs) 



;il.4.5) 
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where the conformal, Ais, and active, A^s, masses od the body B are related to each other 
via Eq. (9.3.5), that is 



MB = MB + e^ 
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and the gravitational potential 

V{x) = U{x)+e^ 



;il.4.7) 



The tidal force F^y is given by Eq. (11.4.1) and Cij is shown in Eq. (11.4.3). 



We are to calculate all terms in the right side of Eq. (11.4.5) explicitly in terms of body's 
mass, rotational moment of inertia, and spin. Among them the most complicated is the first 
one, that is Vi{xg). By making use of Eqs. (11.3.1) - (11.3.7) we obtain 
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where 



Mc = Mc- t 
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is the active mass of the body C, is the general relativistic mass of the body C defined 
by Eq. (9.3.3) where for the sake of simplicity we assumed F = 0, and 77 = 4/3 — 7 — 3 is the 
Nordtvedt parameter. 



After calculating derivatives from potentials V{t, x), U{t, x) and substituting them into Eq. 
(11.4.5) one obtains the following expression for acceleration of the center of mass of the 
body B: 
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where is the Newtonian force and F^^^, F^, F^ are the post-Newtonian relativistic cor- 
rections. Gravitational forces in the right side of this equation are given by the following 
expressions 
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T3<iL> A/r p<L> 



BC iJ^^C 



where we have defined -R' 



BC 



Rf 



I, and spin 5^ of body B relates to 



the angular speed of the body's rotation fi^ via Eq. (11.4.4). 



Eq. (11.4.10) elucidates that inertial mass Mg of the body B is simply its general relativistic 
mass given by Eqs. (9.3.3). This mass is conserved (constant) for spherically-symmetric 
bodies as follows from Eq. (9.3.6). The Nordtvedt gravitational mass Tig of the body B 
depends on the gravitational defect of mass multiplied with the Nordtvedt parameter rj = 



47-/9-3 

= Ms - -e^T] I p*UBd^w . (11.4.16) 
2 Jvb 

The Newtonian gravitational force F^, Eq. (11.4.11), depends in the scalar-tensor theory only 
on the Nordtvedt gravitational masses of the bodies. Will [20] distinguishes the "active" and 
"passive" gravitational masses which depend on the entire set of the PPN parameters. In 
our approach used in the present paper only two PPN parameters, /3 and 7, exist. In this 
case the "active" and "passive" gravitational masses coincide and reduce to one and the 
same expression given by Eq. (11.4.16). Inertial and gravitational masses of the body are 
not equal in the scalar-tensor theory of gravity [38,143,144]. This inequality violates the 
strong principle of equivalence for massive bodies. This violation can be also explained from 
the point of view of interaction of the gravitational field of the body under consideration 
with the scalar field generated by external bodies. This interaction leads to a local force 
which brings about a non-zero value of the time derivative of the body's linear momentum. 
Indeed, assuming that the body under consideration has finite size, is non-rotating, and 
spherically-symmetric, one obtains from Eq. (9.4.1) 



r = MQi (1 + ^e^T] 1^ p*U^^'d'w^ - (11.4.17) 
2 Jvb 

where Pi is a gradient of the scalar field of the external bodies. If one keeps the body's 
center of mass at the origin of the local coordinate system (P* = 0), then, the body's center 
of mass will experience an anomalous acceleration Qi 7^ 0. This anomalous acceleration is 
due to the interaction of the gravitational energy ("gravitational charge") of the body under 
consideration with the gradient of the external scalar field. The coupling constant of this 
interaction is the dimensionless Nordtvedt parameter 77. 

The post-Newtonian forces, Eqs. (11.4.12)-(11.4.15), depend only on the general relativistic 
masses of the bodies which coincide with the Newtonian definition of mass, Eq. (9.1.1) in 
the approximation under consideration. The post-Newtonian force (11.4.12) is known as the 
(Lorentz-Droste) Einstein-Infeld-Hoffmann (EIH) force [149] presently used as a basis of JPL 
ephemerides [21]. It was derived in general relativity by Lorentz and Droste [150] and later 
by Einstein, Infeld and Hoffmann [139], Petrova [140], and Fock [48]. In the Brans-Dicke 
theory this force was derived by Estabrook [151] in case of /3 = 1,7 7^ 1 and by Dallas [152] 
in the case /5 7^ 1,7 7^ 1 (see also [153]). These derivations assumed that the bodies have 
negligible ratio of their radii to the characteristic distance between the bodies (a point-like 
body approximation) as well as that they are non-rotating and move along geodesic world 
lines. 

Corrections to the EIH force are given by Eqs. (11.4.13) -(11.4.15). The force given by 
Eq. (11.4.13) describes the relativistic post-Newtonian correction to the EIH force due to 
the coupling of the body's spin with orbital angular momentum and rotational spins of other 
bodies. It depends on the PPN parameter 7 only. If one takes 7 = 1 in Eq. (11.4.13), the 
force -Fj is reduced exactly to its general relativistic expression obtained earlier by other 
researchers [70,127,154,155]. Our Eq. (11.4.13) for the PPN force coincides with that 
derived by Klioner and Soffel [1]. 

The force (11.4.14) describes general relativistic correction to the EIH force due to the finite 



size of the bodies. This correction is proportional to the forth-order rotational moments of 
inertia of the bodies, X^^-* , while all terms, which are proportional to the second-order body's 
rotational moments of inertia, X^'^\ cancelled mutually out. Nordtvedt [156] considered the 
problem of translational motion of extended bodies in the general class of scalar-tensor 
theories of gravity. He came to the conclusion that covariant formulation of the variational 
principle requires the second-order moment of inertia of extended body to be coupled with the 
Ricci tensor of the background gravitational field generated by the external bodies. Hence, 
such coupling must disappear in general theory of relativity by virtue of the vanishing of 
the background Ricci tensor in vacuum. However, body's moments of inertia of higher order 
couple with the full Riemann tensor and its derivatives, and for this reason they can present 
in general relativistic equations of motion of spherically-symmetric bodies as demonstrated 
in Eq. (11.4.14). 

Nordtvedt's calculation [156] of the equations of motion of extended spherically-symmetric 
bodies agrees with our derivation of these equations based on the implementation of the 
matched asymptotic expansion technique and separate solution of the internal and exter- 
nal problems for gravitational field of the N-body system. However, vanishing of all terms 
depending on the second order rotational moment of inertia in general relativity is in dis- 
agreement with calculations of Brumberg [127], Spyrou [128], Dallas [152] and Vincent [153] 
who came to the conclusion that the general relativistic Lagrangian for the system of N 
spherically-symmetric bodies must depend on the second-order moments of inertia of these 
bodies, X^'^\ Brumberg's expression for the force due to the finite size of the bodies is [127] 
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We have analyzed and pinned down the origin of the disagreement between Nordtvedt's 
[156] and Brumberg's calculations. Nordtvedt [156] used covariant approach while Brumberg 
[127] followed Fock's method [48] and operated with the coordinate-dependent definitions 
of the multipole moments. In fact, Brumberg defined multipole moments of the bodies in 
the global (barycentric) frame of the N-body system. Spherical symmetry of the bodies was 
also defined by Brumberg in the global frame. Such definition of the spherical symmetry 
does not comply with the relativistic law of transformation between local and global frames. 
Hence, the bodies can sustain their spherically-symmetric shape in the global frame if and 
only if there are internal stresses in the body's matter which compensate for the Lorentz and 
gravitational contractions of the body's shapes [138] . Existence of such coordinate-dependent 
internal stresses inside the body under consideration is unnatural. Furthermore, assumption 
about spherical symmetry of the bodies in the global barycentric frame brings about the 
force -F* . which is a pure coordinate effect. 

Brumberg X 



In order to prove that the force -^Brumborg has no physical origin we have considered transla- 
tional equation of motion for body B defined in the global frame as follows 

where the post-Newtonian corrections (not shown explicitly) include the relativistic point- 
like effects and Brumberg's force (11.4.18). Barycentric velocity of the body's matter is 
obtained by differentiation with respect to time of both sides of Eq. (11.2.3) and can be 
decomposed to the sum of the barycentric velocity of the body's center of mass, Vg, and the 
internal velocity's field in local coordinates 



v' = v'g{t) + u\u,w) + e^Au\u,w) , (11.4.20) 

where Au^ is the relativistic correction to the local velocity u^. This correction is a quadratic 
function of the local coordinates of the body. Hence, subsequent calculation of the time 
derivative of and calculation of the integral in the left side of Eq. (11.4.19) bring about 
terms which actually depend on the moment of inertia of body B. This moment of inertia 
is reduced to the rotational moment of inertia in case of a spherically-symmetric body. 
Calculation of the integral from the Newtonian potential in the right side of Eq. (11.4.19) has 
to be done by splitting the potential in two parts - internal and external (see Eq. (8.2.1)), and 
applying Eqs. (11.3.1) and (11.2.8), (11.2.9) for calculation of the integrals from the external 
potential. This again gives a number of terms depending on the rotational moment of inertia 

(2) 

Ib of the body B. Summing up all these terms one obtains exactly the same expression 
as in Eq. (11.4.18) but with opposite sign. It means that these terms cancel out with the 
force -^Brumbcrg- ^hls Completely agrees with our calculations of the force F^^j^ which does 
not depend on the body's rotational moments of inertia of the second order. We conclude 
that the origin of the coordinate-dependent force . is directly associated with an 
inappropriate choice of the body's center of mass and the property of its spherical symmetry 
which must be defined with respect to the local coordinate frame co-moving with the body 
under consideration. Brumberg attempted to make more physical calculation of the force 
F^Q^ in his another book [131] but he did not arrive to any definite conclusion regarding 
whether the force -Fjc^ depends on the body's rotational moments of inertia of the second 
order or not. Our calculations resolve the problem and demonstrate that force . does 

j~ Hrumberg 

not exist and that general relativistic correction due to the finite size of the moving bodies 
is proportional to the forth and higher order rotational moments of inertia of the bodies. 
These corrections are extremely small for the bodies comprising the solar system and can 
be neglected in treatment of the solar system gravitational experiments. However, finite size 
effects can become important during the final stage of coalescence of binary neutron stars 
so that they should be included in the precise calculation of templates of the gravitational 
waveforms. 

The force SF^^^ describes relativistic correction due to the finite size of the bodies in the 
scalar-tensor theory of gravity. This force is proportional to the parameter (3 — 1 only and, 
in contrast to general relativity, depends on the second order rotational moments of inertia, 
X^i\ This dependence was noticed by Nordtvedt [156,157] who has found that in the case 
of weakly self-gravitating bodies the finite-size effects are proportional to ?7 = 4/3 — 7 — 3. 



In the sense that i^BrumbcrE physical impact on the motion of the bodies having finite size. 



This is in disagreement with our calculations of the force SF^^^ but one can easily reconcile 
the two formulations. The matter is that Nordtvedt [157] worked in harmonic coordinate 
system defined by the condition da{\/—gg°'^) = while we worked in the quasi-harmonic 
coordinates defined by the Nutku condition da{4'y/—gg°'^) = (see Eq. (3.4.1) in the present 
paper). This leads to two different forms of the transformation between spatial global and 
local coordinates. In harmonic coordinates this transformation reads [1,157] 
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while in the quasi-harmonic coordinates, used in the present paper, we have 
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The two transformations have different dependence on 7 in terms proportional to the accel- 
eration so that 



harmonic = + {l ' l)e' " \<RI) ■ {IIA.2?>) 



It is due to this difference the parameter 7 had appeared in Nordtvedt's calculations of the 
finite size effects and made the parameter 77 characterizing the magnitude of the finite-size 
effects in the harmonic coordinates. However, dependence of the magnitude of the finite-size 
effects on the parameter 7 in Nordtvedt's calculations is a pure coordinate effect which has 
no physical meaning. Parameter 7 can be eliminated from the force SF^^^ if one works in 
the quasi-harmonic coordinates defined by the Nutku condition (3.4.1) of this paper. 



11.5 Rotational Equations of Motion 



Derivation of rotational equations of motion for spherically-symmetric bodies requires cal- 
culation of the multipole moments TZl oi the body under consideration. One has 

7^L = O (e^) . (11.5.1) 

All other multipole moments have been calculated in the previous section. Performing cal- 
culations of the torques and body's spin given by Eqs. (10.2.2)-(10.2.4) one obtains 
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Consequently, the rotational equation of motion for the body's spin is 
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where Ck is angular velocity of rotation of the local coordinate frame with respect to that 
which axes are subject to the Fermi- Walker transport. 

Our Eq. (11.5.5) has one extra term comparatively with the corresponding equation (9.75) 
from the paper [1] by Klioner and Soffel. This term depends on the forth-order rotational 
moment of inertia, x'^b \ of the body B and has pure general relativistic origin. This term 
was not taken into account by Klioner and Soffel because they neglected finite size of the 
rotating body. Contribution of the the forth-order rotational moment of inertia, X^f'' to the 
rotational torque is negligibly small for the theory of Earth's rotation. However, it may 
become significant during last several orbits of a coalescing binary neutron star. It would be 
interesting to study the impact of this term on the form of gravitational waves emitted such 
binaries. 



A Solution of the Laplace Equation for Scalar, Vector and Tensor Fields 



In this appendix we find solutions of the Laplace equations for scalar F(t, a?), vector -Fj(t, x), 
and tensor Fij{t,x) fields. These equations are 

AF{t,x) = , AFi{t,x) = , AFij{t,x) = 0. (A.l) 

The procedure of finding solutions of equations (A.l) is based on the approach developed in 
[137] (see also [62], [83] and references therein). 

Basic spherical functions are 

¥^"^{9, 0) = -^e''^'^PJ^{cosd) , (-/ < m < /), (A.2) 



where Pp^{cos6) are the associated Legendre polynomials. According to the general theorem 
any arbitrary function f{6, (p) such that its square is integrable over the sphere, can be 
expanded in the convergent series 
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Making use of transformation from spherical to Cartesian coordinates one can obtain [62] 
a one-to-one mapping between the spherical harmonics and the symmetric trace-free (STF) 
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where N^^ 



are products of components of unit vector n = r/r. Tensors 3^<^ 



with —l<m<l form a basis in {21 + l)-dimensional space of symmetric and trace-free 
tensors with / indices, that is any STF tensor of rank / can be represented as 
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Hence, equation (A. 3) can be recast to the following form 

oo 



(A.5) 



(A.6) 



Spherical functions are the eigenf unctions of the orbital angular momentum operator 
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that is a consequence of definition of canonical basis in the sub-space in which the irreducible 
representation with weight I is realized. 

Thus, equations (A.6) and (A.7) reveal that any scalar function F{t, x) that is solution of 
the Laplace equation (A.l) is given by 
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where v4<i> and B^l^ are STF multipole moments depending on time t only. 



Vector and tensor spherical harmonics are obtained from the direct product of two irreducible 
representations of the rotation group with weights I' and /" which can be expanded into 
irreducible representations with weights \l' — l"\ < I < \l' + Canonical orthonormal basis 
in the sub-space of the vector spherical functions in which the irreducible representation has 
weight /, is the set of 3(2/ + 1) functions 
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where /' can take either of three values 1 — 1,1,1 + 1, notation m"m' \ Im) stands for the 
Clebsch-Gordan coefficients [137], and three vectors 



represent canonical basis of the main matrix representation. 
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Canonical basis in the three-dimensional space of constant tensors of second rank is made 
of 9 tensors such that five of them, 
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are symmetric trace-free tensors, three tensors, 
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are fully antisymmetric, and one, 
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m' = ~l' m"=—l 



is the unit tensor. We shall abandon the antisymmetric part of the basis as we are interested 
only in the symmetric tensors. Hence, the canonical basis in the sub-space of such tensors of 
the second rank with irreducible representation with weight / is formed from 6(2/ + 1) tensor 
harmonics 
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where / — 2 </'</ + 2, and 2/ + 1 scalar spherical harmonics 
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Finally, solutions of the Laplace equations for vector, Fj, and tensor, Fjj, functions (A.l) are 
given as follows 
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Fij (t, 3j) 
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where the symbol Sym(ij) denotes symmetrization, and Cl-, Dl,---, Tl represent STF mul- 
tipole moments depending on time t only. 



B The ChristofFel Symbols and the Riemann Tensor 



B.l The Christoffel Symbols 



In this appendix we give formulas for the Christoffel symbols and the Riemann tensor which 
elucidates the physical meaning of external multipole moments in the expression for the local 
metric tensor (6.2.36) - (6.2.38). Christoffel symbols are defined by standard expression 

where the metric tensor components are taken from equations (3.3.5) - (3.3.7). Calculation 
results in 
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We neglect in our calculations all terms which are quadratic with respect to the angular 

(1) 

speed of rotation f2* and linear velocity V* in hoi- In addition, we note many terms in Eqs. 

(1) 

(B.1)-(B.7) equal zero due to the specific structure of hoi- Thus, after simplification the 
expressions for the Christoffel symbols are reduced to 
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These expressions have been used in the present paper. 



B.2 The Riemann Tensor of External Gravitational Field in Loeal Coordinates 



Components of the Riemann tensor computed by making use of the external metric tensor 
only, are 
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Here, all multipole moments Pl are caused by the presence of scalar field. 



C Comparison with the Khoner-Soffel Approach to Reference Frames in the 
PPN Formalism 

Klioner and Soffel [1] have worked out an independent approach aimed to amend the standard 
PPN formalism [18,19,20] with a mathematical procedure allowing us to construct the local 
PPN coordinates in the vicinity of moving massive bodies (Earth, Moon, etc.) comprising 
an N-body system. Klioner-Soffel approach is a straightforward extension of the standard 
PPN formalism and, thus, is not based on the field equations of a particular parametrized 
theory of gravity in order to preserve that generality which the original PPN formalism 
was trying to reach in the global PPN coordinates. Our point of view is that such general 
PPN formulation is not possible without having recourse to the field equations of a particular 
theory of gravity (or a set of theories). This is because the goal of gravitational experiments is 
to explain physics of gravitational field and exclude the field equations of untenable theories 
of gravity. Formal procedure for extension of the PPN formalism, proposed by Klioner and 
Soffel, is not able to solve this task. Even in the case of two parameters, 7 and (3, one can 
not say what physical theory is hiding behind mathematical manipulations present in the 
Klioner-Soffel approach [1]. 

The two alternative points of view on one and the same problem led naturally to a number of 
differences in many equations which have been derived. In this appendix we summarize the 
differences we have found between the Klioner-Soffel formulation [1] of the local reference 
frames in the PPN formalism and that based directly on the scalar-tensor theory of gravity 
used in the present paper. These differences are as follows: 

• Klioner and Soffel restricted their attention to the case of two PPN parameters, j3 and 7 like 
we did. However, Klioner and Soffel did not rely upon any specific theory of gravitation, 
thus, leaving open the question about what kind of field equations is compatible and 
can be used along with the PPN coordinate transformations and the equations of motion 
they have obtained. Calculations in the present paper have been done in the framework 
of a general class of scalar-tensor theories of gravity with one scalar field which is also 
parametrized by the PPN parameters 7 and (3. In our formulation we know exactly the 
field equations which are compatible with our coordinate transformations and equations of 
motion. Hence, measuring these parameters restricts domain of applicability of a particular 
physical theory of gravity. 

• Klioner and Soffel worked in harmonic coordinates defined by the condition da{^/—gg°'^) = 
while we worked in the quasi-harmonic coordinates defined by the Nutku condition 
doi{,4'\/~99"^) = . Spatial transformation from local to global coordinates given by 
Eqs. (8.6.2) in our approach does not contain the PPN parameter 7 in the quadratic terms 
with respect to i?^ while that of Klioner and Soffel does. 

• Klioner and Soffel postulated (guessed) a specific form of the metric tensor in the local 
frame with one free function - potential ^ entering their Eqs. (3.33) and (4.25). Then, 
they postulated specific rules for matching of the local and global frames and used them in 
order to determine the post-Newtonian coordinate transformation along with the a priory 
unknown potential ^ in the local frame. We have used the field equations in order to find 
the metric tensor in the local frame, so that it gets fully determined, and matched it with 
the metric tensor in the global frame by making use of the same procedure as that used 
in general theory of relativity [14,15]. Hence, our matching procedure is not voluntary but 



See Eq. (3.4.1) in the present paper. 



is in a complete agreement with the field equations. 

• Klioner and Soff'el noticed two distinguished properties of the metric tensor in the local 
frame valid in general relativity: 

A. Gravitational field of external bodies is represented only in the form of a relativistic tidal 
potential being at least of the second order in the local spatial coordinates. It coincides 
with the usual Newtonian tidal potential in the Newtonian limit of the PPN formalism; 

B. The internal gravitational field of a body (or sub-system of bodies) coincides with the 
gravitational field of a corresponding isolated source provided that the tidal infiuence of 
the external matter is neglected. 

Klioner-Soffel approach can not keep the two properties of the metric tensor in the local 
frame simultaneously while in our approach to the problem under consideration these two 
properties are retained similar to general relativity. 

• Klioner and Soffel do not distinguish active, conformal, and scalar multipole moments 
characterizing different properties of the gravitational field of an isolated body. In fact, 
they operate only with active multipole moments. We use two types of the multipole mo- 
ments in our calculations and show that the conformal dipole moment has more advanced 
mathematical properties than the active one for definition of body's center of mass. It 
allows us to simplify translational equations of motion of the bodies drastically. 

• Klioner and Soffel derived equations of motion of a body in its local frame and discovered 
that in their approach these equations depend explicitly not only on the set of multipoles 
being inherent to general theory of relativity but contain one more family of the multipole 
moments which they denoted as A/l . In our approach the equations of motion of the 
body in the local frame do not explicitly depend on such multipole moments except for the 
second time derivative of the monopole term, J\f = X^'^\ entering Eq. (9.4.7). The reader 
is invited to compare equations (9.37)-(9.40) from the paper [1] with equations (9.4.1) - 
(9.4.7) of the present paper to decide which type of equations is more economic. 

• Klioner and Soffel defined the center of mass of a body by making use of active dipole mo- 
ment of the body (see their Eq. (9.24)). They have found that if the Nordtvedt parameter 
?7 7^ and/or 7 7^ 0, the second time derivative of this dipole does not vanish for the case of 
one isolated body unless specific physical conditions are met (secular stationarity). Thus, 
Klioner and Soffel admit the existence of a self-accelerated motion of the body in their for- 
malism which violates the third (action-counteraction) Newton's law. On the other hand, 
we defined the center of mass of each body by making use of conformal dipole moment of 
the body (see our Eq. (9.3.7). The second time derivative of the con/orma/ dipole moment 
vanishes perfectly in the case of a single body and no self-accelerated terms in motion of 
the body comes about in our formalism. 

• Equations of rotational motions in Klioner-Soffel formalism contain torques which depend 
on the Nordtvedt parameter 77 = 4 — 3/3 — 1. We confirm this observation. However, we 
noticed that Klioner and Soffel did not take into account finite-size effects of the bodies and 
they worked in a different gauge. This led to appearance of several terms in our equations 
which are not present in Klioner-Soffel's analysis [1]. These terms can be important in 
some particular situations, for example, in rotational equations of motion of coalescing 
binaries. 



See our Eq. (9.3.8) for its precise definition. 
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Fig. 1. The picture illustrates an astronomical N-body system and coordinate charts associated 
with it. Global coordinates, x" = {ct,x^), cover the entire space, have origin at the barycenter 
of the system, and are Minkowskian at infinity. Each body has its own local coordinate chart, 
w" = {cu, w^), having origin at the center of mass of the body under consideration. Local coordinates 
are not asymptotically Minkowskian far away from the body and do not cover the entire space. 
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Barycenter of the subsystem 




Body #2 



Barycenter of the N-body system 



Body #3 




Barycenter of the subsystem 




Body #5 
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Fig. 2. The picture illustrates a hierarchy of coordinate charts existing in N-body astronomical 
system (solar system) consisting of several sub-systems: Mi, M2,...,Mp. One global coordinate 
chart, x°' = {ct, x^), covers the entire space and has origin at the barycenter of the N-body system. 
Each sub-system consists of several gravitationally-bounded bodies (a planet and its moons) and 
has its own local coordinate chart, w"' = {cu,w''), having origin at the center of mass of the 
sub-system. At the same time, each body from the sub-system possesses its own local coordinate 
chart, = (cs, ^*). The hierarchy can have as many levels as necessary for adequate description of 
motion of the bodies. 



Body's World Tube 




Fig. 3. Two 3-dimensional hypersurfaces of constant time related to the global (G) and local (L) 
coordinates are shown. The two hypersurfaces do not coincide because of relativistic transformation 
of the local coordinate grid with respect to the global one. A body's world tube contains a world line 
of the origin of the local coordinates. The world tube cross-section, ^i, is a part of the hypersurface 
of constant time t of the global coordinates, and that, ^42, is a part of the hypersurface of constant 
time u of the local coordinates. The line CD marks intersection of the two hypersurfaces. All 
integrals depending on density, velocity and other internal characteristics of the body are performed 
in the global coordinates over the cross-section Ai whereas in the local coordinates the integration 
is over the cross-section A2- 




' = const. + zt 



t = const. 



Fig. 4. Lie Transfer 

Matching the local and global coordinates requires to calculate integrals of various functions 
depending on the internal structure of body over two different hypersurfaces of constant 
time, u and t, as shown in Fig. 3. Relationship between the integrals taken on the two 
hypersurfaces is set up by making use of the Lie transfer from one hypersurfaces to another. 
The integral curves used in the Lie transfer are the world lines of the four-velocity of body's 
matter. Only one of such lines, MN, is shown in the figure. Dashed line is the world line of 
the origin of the local coordinates which coincides with the center of mass of the body under 
consideration. The center of mass is not associated with four-velocity of body's matter. 




Fig. 5. World lines of the origins of two local coordinates are shown. The line (A) indicates the world 
line of the origin of the local coordinates A falling freely in the background space-time described by 
the external metric tensor g^^ . The line (B) depicts motion of the origin of the local coordinates 
B always located at the center of mass of the body under consideration. In general, the two world 
lines do not coincide due to the existence of the inertial force, AiQi, in the local coordinate system 
B. This force arises due to the gravitational coupling between mass and spin multipole moments 
of the body with the gravitoelectric and gravitomagnetic tidal gravitational fields [106,135] of the 
background metric produced by other bodies of the N-body system. 



